In this study we explored the mode of genetic inheritance of Varroa destructor, a parasitic mite of honeybees, so far known as a haplodiploid species with a arrhenotokous parthenogenesis reproductive mode. The mite shows contrasting phenomena of highly inbreeding life style (sib-mating), yet maintaining relatively high genetic variation.

The experimental setup consisted by a three-generational pedigree with sib-mating.
Sample size: 30 families, total of 223 individuals 3 different colonies, from OIST apiary
For details of mite collection and pedigree construction, please see the original manuscript.

load libraries

library("tidyverse")
library("plyr")
library("dplyr")
library("ggplot2")
library("scales")
library("ggpubr")
library("gridExtra")
library("grid")
library("GGally")
library("vcfR") # for extracting genotype data from a vcf file
library("data.table")
library("stringr")
library("janitor")
library("gmodels")
library("rstatix")
library("freqtables")
library("broom")
library("DescTools") # for the Goodness-of-Fit test, 3 variables; and for Breslow-Day Test for Homogeneity of the Odds Ratios
library("vcd") # for the woolf test testign log odds for each pair
library("patchwork") # for gathering the plots
#library("fuzzyjoin") # to join tables based on a string in a column
#library("aspi") # Repeated G–tests of Goodness-of-Fit, work only for 2 variables..
#library("RVAideMemoire") # Repeated G–tests of Goodness-of-Fit, work only for 2 variables...
#library("InfoTrad")
#library("ggthemes") # for more colors in the ggplot
knitr::opts_chunk$set(echo = TRUE, warning = FALSE, message = FALSE,
                     fig.width = 8,
                      fig.asp = 0.6,
                      out.width = "100%")
#fig.width = 6,fig.asp = 0.8,out.width = "100%"

Load Variant Call Format (VCF) file.

Extract genotypes for each site and individual. The metadata for all samples can be found in here.

vcf <- read.vcfR("/Users/nuriteliash/Documents/GitHub/varroa-linkage-map/data/vcf_filter/Q40BIALLDP16HDP40mis.5Chr7/Q40BIALLDP16HDP40mis.5Chr7.recode.vcf", verbose = FALSE )
vcf
## ***** Object of Class vcfR *****
## 223 samples
## 7 CHROMs
## 35,169 variants
## Object size: 293.5 Mb
## 0 percent missing data
## *****        *****         *****
# extract the genotype for each site in each individual
gt <- extract.gt(vcf, element = "GT") 
gt <- as.data.frame(t(gt)) %>%
  rownames_to_column("ID")
#clean the ID names
gt$ID <- sub("_[^_]+$", "", gt$ID)

The individual ID name nomenclature is composed of 3 parts, separated by an underscore (_).
The first part is the family ID, the second is the unique name of the individual, and the the third part indicates its generation, sex, and its position in relation to the foundress mite (generation F0).
For example:
Individual ID 240_240a_son, belong to family 240, has a unique name of 240a and is the son of the foundress mite of family 240, that is, its a male of F1 generation.
Individual ID 240_241b_grndat, belong to family 240, has a unique name of 241b and is the granddaughter of the foundress mite of family 240, that is, its a female of F2 generation (and the daughter of 240_240a_son).


From the 223 individuals of 30 families, we excluded non adults mites (for which its hard to determine the sex).
Eventually, we kept 144 individuals (adult males and females in F0, F1 and F2 generations) of 26 families, and all 35,169 biallelic sites.


For example, these are the first 6 sites (13565-25513) of chromosome (NW_019211454.1), in one family (family number 240).

table <-  gt %>% 
  t() %>%
  as.data.frame() %>%
  row_to_names(row_number = 1) %>% 
  dplyr::select(contains(c("son", "dat", "fnd"))) # keep only adults of F0, F1 and F2 

table %>% select(contains(c("240_"))) %>% head()
##                      240_240a_son 240_241c_grnson 240_241b_grndat
## NW_019211454.1_13565          0/1            <NA>            <NA>
## NW_019211454.1_18037          0/0            <NA>            <NA>
## NW_019211454.1_20998         <NA>            <NA>            <NA>
## NW_019211454.1_24935         <NA>             0/0            <NA>
## NW_019211454.1_25470         <NA>            <NA>            <NA>
## NW_019211454.1_25513         <NA>            <NA>            <NA>
##                      240_241a_grndat 240_241_dat 240_240_fnd
## NW_019211454.1_13565             0/1         0/1         0/0
## NW_019211454.1_18037             0/0         0/0         0/0
## NW_019211454.1_20998             0/0         0/0         0/0
## NW_019211454.1_24935             0/0         0/0         0/0
## NW_019211454.1_25470             0/0         0/0         0/0
## NW_019211454.1_25513             0/0        <NA>         0/0

The family members include:

  • F0 generation: foundress female mite (240_240_fnd).
  • F1 generation: adult son (240_240a_son) and adult daughter (240_241_dat) of the foundress F0 mite. Because in varroa mite reproduction is via sib-mating, these brother and sister will also mate, to produce the F2 generation.
  • F2 generation: adult grandson (240_241c_grnson) and two adult granddaughters (240_241a_grndat and 240_241b_grndat), of the foundress F0 mite.

Each individual, was genotyped for each site with one of the three genotypes:

  • Homozygote for the reference allele = 0/0
  • Heterozygote = 0/1
  • Homozygote for the alternate allele = 1/1
  • “NA” = site genotype not determined

For more information about the mapping, variant calling and variant filtration, please see the Snakemake pipeline in here.


In the following code we viewed the F2 offspring genotype of all possible nine crosses between F1 male and females.

  • Control crosses of homozygotic sites: The first four crosses are aiming mainly to detect false sites:
    1. F1 male (0/0) x female (0/0)
    2. F1 male (1/1) x female (1/1)
    3. F1 male (0/0) x female (1/1)
    4. F1 male (1/1) x female (0/0)

Then we cross heterozygotic sites, to explore the mode of genetic inheritance in varroa mite:

  • Homozygotic male, with heterozygotic female:
    1. F1 male (0/0) x female (0/1)
    2. F1 male (1/1) x female (0/1)
  • Heterozygotic male, with homozygotic female:
    1. F1 male (0/1) x female (0/0)
    2. F1 male (0/1) x female (1/1)
  • Heterozygotic male, with heterozygotic female:
    1. F1 male (0/1) x female (0/1)

For each crossing we have expected genotype-ratio of the F2, based on the following assumptions under arrhenotokous haplodiploid system:

  • Fertilized egg develops into diploid female ♀ (2n).
  • Unfertilized egg develops into haploid male ♂ (n).
  • No paternal inheritance in the males.

Control: Homozygotic sites (crosses 1-4)

To detect false sites

(1) F1 male (0/0) x female (0/0)

F2 offspring

we expect all F2 offspring to be homozygotic (0/0) like their parents (F1)

x <- dat %>%
  dplyr::select(c("sample","gt","sex", "exp","prop", "total")) %>%
  mutate(type = "exp") %>%
  dplyr::rename(count = exp)

y <- dat %>%
  dplyr::select(c("sample","gt","sex", "obs","total")) %>%
  mutate(type = "obs") %>%
  dplyr::rename(count = obs) %>%
  mutate(prop = (count/total))
  
df_male <- rbind(x, y) %>%
  dplyr::filter(sex == "male")
      
df_fem <- rbind(x, y) %>%
  dplyr::filter(sex == "female")

p_male_pooled <- ggplot(df_male, aes(fill=gt, y=prop, x=type)) + 
    geom_bar(position="fill", stat="identity", ) +
    xlab("F2 males") + 
    ylab("Proportion of F2 genotype") +
    labs(title = "F2 Males of F1 male (0/0) x female (0/0)", 
         subtitle = paste0("Pooled sites = ", df_male %>% dplyr::filter(gt == "0/1") %>% dplyr::filter(type == "obs") %>% summarise(sum(total)) %>% as.numeric())) + 
   scale_x_discrete("type", labels = c("Expected","Observed")) +
labs(fill = "Genotype") + 
  theme_classic() +
  theme(axis.text.x = element_text(size=15),
        axis.title.x = element_blank()) 

p_female_pooled <- ggplot(df_fem, aes(fill=gt, y=prop, x=type)) + 
    geom_bar(position="fill", stat="identity", ) +
    xlab("F2 females") + 
    ylab("Proportion of F2 genotype") +
    labs(title = "F2 females of F1 male (0/0) x female (0/0)", 
         subtitle = paste0("Pooled sites = ", df_fem %>% dplyr::filter(gt == "0/1") %>% dplyr::filter(type == "obs") %>% summarise(sum(total)) %>% as.numeric())) + 
    scale_x_discrete("type", labels = c("Expected","Observed")) +
labs(fill = "Genotype") + 
  theme_classic() +
  theme(axis.text.x = element_text(size=15),
        axis.title.x = element_blank()) 

grid.arrange(p_female_pooled, p_male_pooled, ncol = 2)

# plot ratio for each individual: 
# males:
ggplot(df_male, aes(fill=gt, y=prop, x=type)) + 
    geom_bar(position="fill", stat="identity", ) +
    xlab("F2 siblings") + 
    ylab("Proportion of F2 genotype") +
    labs(title = "F2 Males of F1 male (0/0) x female (0/0)") + 
   # geom_text(aes(label = n), position = position_stack(vjust = 0.5)) +
    geom_text(data = filter(df_male, gt == "1/1"), aes(y=prop, x=type, label=total), vjust = 0) +
    facet_wrap(~ sample) +
  labs(fill = "Genotype") + 
  theme_classic() +
  theme(axis.text.x = element_text(size=15),
        axis.title.x = element_blank()) 

# and females:
ggplot(df_fem, aes(fill=gt, y=prop, x=type)) + 
    geom_bar(position="fill", stat="identity", ) +
    xlab("F2 siblings") + 
    ylab("Proportion of F2 genotype") +
    labs(title = "F2 Females of F1 male (0/0) x female (0/0)") + 
   # geom_text(aes(label = n), position = position_stack(vjust = 0.5)) +
    geom_text(data = filter(df_fem, gt == "1/1"), aes(y=prop, x=type, label=total), vjust = 0) +
    facet_wrap(~ sample) +
  labs(fill = "Genotype") + 
  theme_classic() +
  theme(axis.text.x = element_text(size=15),
        axis.title.x = element_blank()) 

F0 foundress

we expect the F0 foundress to be homozygotic (0/0), like her son

# plot the pooled data for all F0 foundress with male (0/0) x female (0/0) offspring
ggplot(observed, aes(fill=gt, y=prop, x = type)) + 
    geom_bar(position="fill", stat="identity") +
    ylab("Proportion of F0 genotype") +
    labs(title = "F0 foundress of F1 male (0/0) x female (0/0)", 
         subtitle = paste0("Pooled sites = ", observed %>% dplyr::filter(gt == "0/1") %>% summarise(sum(total)) %>% as.numeric())) + 
    theme_classic() +
    theme(axis.text.x = element_blank(), axis.title.x = element_blank()) + 
    labs(fill = "Genotype")

# plot for each individual: 
ggplot(observed, aes(fill=gt, y=prop, x=type)) + 
    geom_bar(position="fill", stat="identity", ) +
    xlab("Observed") + 
    ylab("Proportion of F0 genotype") +
    labs(title = "F0 foundress of F1 male (0/0) x female (0/0)") + 
    geom_text(data = filter(observed, gt == "1/1"), aes(y=prop, x=type, label=total), vjust = 0) +
    facet_wrap(~ sample) +
  labs(fill = "Genotype") + 
  theme_classic() +
  theme(axis.text.x =  element_blank(), axis.title.x = element_blank()) 

(2) F1 male (1/1) x female (1/1)

F2 offspring

we expect all F2 offspring to be homozygotic (1/1) like their parents (F1)

# set the families (include only families with at least one adult F2)
family = grep("grndat|grnson",gt$ID, value=TRUE) %>%
  str_extract("[^_]+")  %>%
  unique()

# define a list to put all the data frames in
obs <- list()

for (fam in family) {
 
obs[[fam]] <- table %>%
  dplyr::select(starts_with(fam)) %>%
  dplyr::filter_at(vars(matches("_son")), all_vars(. == "1/1")) %>%
  dplyr::filter_at(vars(matches("_dat")), all_vars(. == "1/1")) %>% 
  dplyr::select(contains("grn")) %>%
  tidyr::pivot_longer(everything())  %>% 
  dplyr::rename(sample = name, gt = value) %>%
  dplyr::count(sample, gt, .drop = FALSE) %>% 
  dplyr::filter(gt %in% c("0/0", "1/1", "0/1")) %>%
  mutate(n = as.numeric(n)) %>%
  group_by(sample) %>%
  mutate(total = as.numeric(sum(n))) %>%
  dplyr::rename(obs = n) %>%
  mutate(sex = case_when(
    grepl("son", sample) ~ "male",
    grepl("dat", sample) ~ "female"))
}

# bind all families together, to a final data frame containing all observed counts
observed <- do.call("rbind", obs)

# state the expected site proportion of each genotype and sex under a perfect mendelian segregation of F2
fem_exp <- data.frame(sex = "female", gt = c("0/0", "0/1", "1/1") , prop = c(0,0,1))
male_exp <- data.frame(sex = "male",gt = c("0/0", "0/1", "1/1"),prop = c(0,0,1))
mendel <- bind_rows(fem_exp,male_exp) %>% mutate(sex = as.character(sex))%>% mutate(gt = as.character(gt)) 
 # mutate(type = "exp")

# using on the total genotypes count of sites per sample, calculate the expected counts, and add them to the observed counts into a final dat table
dat <- observed %>%
  #mutate(sample = as.character(sample)) %>%
  mutate(gt = as.character(gt)) %>%
  left_join(mendel, by = c("gt", "sex")) %>% 
  mutate(exp = total*prop)
x <- dat %>%
  dplyr::select(c("sample","gt","sex", "exp","prop", "total")) %>%
  mutate(type = "exp") %>%
  dplyr::rename(count = exp)

y <- dat %>%
  dplyr::select(c("sample","gt","sex", "obs","total")) %>%
  mutate(type = "obs") %>%
  dplyr::rename(count = obs) %>%
  mutate(prop = (count/total))
  
df_male <- rbind(x, y) %>%
  dplyr::filter(sex == "male")
      
df_fem <- rbind(x, y) %>%
  dplyr::filter(sex == "female")

p_male_pooled <- ggplot(df_male, aes(fill=gt, y=prop, x=type)) + 
    geom_bar(position="fill", stat="identity", ) +
    xlab("F2 males") + 
    ylab("Proportion of F2 genotype") +
    labs(title = "F2 Males of F1 male (1/1) x female (1/1)", 
         subtitle = paste0("Pooled sites = ", df_male %>% dplyr::filter(gt == "0/0") %>% dplyr::filter(type == "obs") %>% summarise(sum(total)) %>% as.numeric())) + 
   scale_x_discrete("type", labels = c("Expected","Observed")) +
labs(fill = "Genotype") + 
  theme_classic() +
  theme(axis.text.x = element_text(size=15),
        axis.title.x = element_blank()) 

p_female_pooled <- ggplot(df_fem, aes(fill=gt, y=prop, x=type)) + 
    geom_bar(position="fill", stat="identity", ) +
    xlab("F2 females") + 
    ylab("Proportion of F2 genotype") +
    labs(title = "F2 females of F1 male (1/1) x female (1/1)", 
         subtitle = paste0("Pooled sites = ", df_fem %>% dplyr::filter(gt == "0/0") %>% dplyr::filter(type == "obs") %>% summarise(sum(total)) %>% as.numeric())) + 
    scale_x_discrete("type", labels = c("Expected","Observed")) +
labs(fill = "Genotype") + 
  theme_classic() +
  theme(axis.text.x = element_text(size=15),
        axis.title.x = element_blank()) 

grid.arrange(p_female_pooled, p_male_pooled, ncol = 2)

ggplot(df_male, aes(fill=gt, y=prop, x=type)) + 
    geom_bar(position="fill", stat="identity", ) +
    xlab("F2 siblings") + 
    ylab("Proportion of F2 genotype") +
    labs(title = "F2 Males of F1 male (1/1) x female (1/1)") + 
   # geom_text(aes(label = n), position = position_stack(vjust = 0.5)) +
    geom_text(data = filter(df_male, gt == "1/1"), aes(y=prop, x=type, label=total), vjust = 0) +
    facet_wrap(~ sample) +
  labs(fill = "Genotype") + 
  theme_classic() +
  theme(axis.text.x = element_text(size=15),
        axis.title.x = element_blank()) 

ggplot(df_fem, aes(fill=gt, y=prop, x=type)) + 
    geom_bar(position="fill", stat="identity", ) +
    xlab("F2 siblings") + 
    ylab("Proportion of F2 genotype") +
    labs(title = "F2 Females of F1 male (1/1) x female (1/1)") + 
   # geom_text(aes(label = n), position = position_stack(vjust = 0.5)) +
    geom_text(data = filter(df_fem, gt == "1/1"), aes(y=prop, x=type, label=total), vjust = 0) +
    facet_wrap(~ sample) +
  labs(fill = "Genotype") + 
  theme_classic() +
  theme(axis.text.x = element_text(size=15),
        axis.title.x = element_blank()) 

F0 foundress

we expect the F0 foundress to be homozygotic (1/1), like her son

# plot the pooled data for all F0 foundress with male (1/1) x female (1/1) offspring
ggplot(observed, aes(fill=gt, y=prop, x = type)) + 
    geom_bar(position="fill", stat="identity") +
    ylab("Proportion of F0 genotype") +
    labs(title = "F0 foundress of F1 male (1/1) x female (1/1)", 
         subtitle = paste0("Pooled sites = ", observed %>% dplyr::filter(gt == "0/1") %>% summarise(sum(total)) %>% as.numeric())) + 
    theme_classic() +
    theme(axis.text.x = element_blank(), axis.title.x = element_blank()) + 
    labs(fill = "Genotype")

# plot for each individual: 
ggplot(observed, aes(fill=gt, y=prop, x=type)) + 
    geom_bar(position="fill", stat="identity", ) +
    xlab("Observed") + 
    ylab("Proportion of F0 genotype") +
    labs(title = "F0 foundress of F1 male (1/1) x female (1/1)") + 
    geom_text(data = filter(observed, gt == "1/1"), aes(y=prop, x=type, label=total), vjust = 0) +
    facet_wrap(~ sample) +
  labs(fill = "Genotype") + 
  theme_classic() +
  theme(axis.text.x =  element_blank(), axis.title.x = element_blank()) 

(3) F1 male (0/0) x female (1/1)

F2 offpsring

We expect zero sites for this cross, since there is no paternal inheritance to the males.
Indeed, there are only 103 sites in the F2 pooled females, and 71 for the F2 pooled males:

# set the families (include only families with at least one adult F2)
family = grep("grndat|grnson",gt$ID, value=TRUE) %>%
  str_extract("[^_]+")  %>%
  unique()

# define a list to put all the data frames in
obs <- list()

for (fam in family) {
 
obs[[fam]] <- table %>%
  dplyr::select(starts_with(fam)) %>%
  dplyr::filter_at(vars(matches("_son")), all_vars(. == "0/0")) %>%
  dplyr::filter_at(vars(matches("_dat")), all_vars(. == "1/1")) %>% 
  dplyr::select(contains("grn")) %>%
  tidyr::pivot_longer(everything())  %>% 
  dplyr::rename(sample = name, gt = value) %>%
  dplyr::count(sample, gt, .drop = FALSE) %>% 
  dplyr::filter(gt %in% c("0/0", "1/1", "0/1")) %>%
  mutate(n = as.numeric(n)) %>%
  group_by(sample) %>%
  mutate(total = as.numeric(sum(n))) %>%
  dplyr::rename(obs = n) %>%
  mutate(sex = case_when(
    grepl("son", sample) ~ "male",
    grepl("dat", sample) ~ "female"))
}

# bind all families together, to a final data frame containing all observed counts
observed <- do.call("rbind", obs)

# state the expected site proportion of each genotype and sex under a perfect mendelian segregation of F2
fem_exp <- data.frame(sex = "female", gt = c("0/0", "0/1", "1/1") , prop = c(0,0,0))
male_exp <- data.frame(sex = "male",gt = c("0/0", "0/1", "1/1"),prop = c(0,0,0))
mendel <- bind_rows(fem_exp,male_exp) %>% mutate(sex = as.character(sex))%>% mutate(gt = as.character(gt)) 
 # mutate(type = "exp")

# using on the total genotypes count of sites per sample, calculate the expected counts, and add them to the observed counts into a final dat table
dat <- observed %>%
  #mutate(sample = as.character(sample)) %>%
  mutate(gt = as.character(gt)) %>%
  left_join(mendel, by = c("gt", "sex")) %>% 
  mutate(exp = total*prop)
x <- dat %>%
  dplyr::select(c("sample","gt","sex", "exp","prop", "total")) %>%
  mutate(type = "exp") %>%
  dplyr::rename(count = exp)

y <- dat %>%
  dplyr::select(c("sample","gt","sex", "obs","total")) %>%
  mutate(type = "obs") %>%
  dplyr::rename(count = obs) %>%
  mutate(prop = (count/total))
  
df_male <- rbind(x, y) %>%
  dplyr::filter(sex == "male")
      
df_fem <- rbind(x, y) %>%
  dplyr::filter(sex == "female")

p_male_pooled <- ggplot(df_male, aes(fill=gt, y=prop, x=type)) + 
    geom_bar(position="fill", stat="identity", ) +
    xlab("F2 males") + 
    ylab("Proportion of F2 genotype") +
    labs(title = "F2 Males of F1 male (0/0) x female (1/1)", 
         subtitle = paste0("Pooled sites = ", df_male %>% dplyr::filter(gt == "0/0") %>% dplyr::filter(type == "obs") %>% summarise(sum(total)) %>% as.numeric())) + 
   scale_x_discrete("type", labels = c("Expected","Observed")) +
labs(fill = "Genotype") + 
  theme_classic() +
  theme(axis.text.x = element_text(size=15),
        axis.title.x = element_blank()) 

p_female_pooled <- ggplot(df_fem, aes(fill=gt, y=prop, x=type)) + 
    geom_bar(position="fill", stat="identity", ) +
    xlab("F2 females") + 
    ylab("Proportion of F2 genotype") +
    labs(title = "F2 females of F1 male (0/0) x female (1/1)", 
         subtitle = paste0("Pooled sites = ", df_fem %>% dplyr::filter(gt == "0/0") %>% dplyr::filter(type == "obs") %>% summarise(sum(total)) %>% as.numeric())) + 
    scale_x_discrete("type", labels = c("Expected","Observed")) +
labs(fill = "Genotype") + 
  theme_classic() +
  theme(axis.text.x = element_text(size=15),
        axis.title.x = element_blank()) 

grid.arrange(p_female_pooled, p_male_pooled, ncol = 2)

ggplot(df_male, aes(fill=gt, y=prop, x=type)) + 
    geom_bar(position="fill", stat="identity", ) +
    xlab("F2 siblings") + 
    ylab("Proportion of F2 genotype") +
    labs(title = "F2 Males of F1 male (0/0) x female (1/1)") + 
   # geom_text(aes(label = n), position = position_stack(vjust = 0.5)) +
    geom_text(data = filter(df_male, gt == "1/1"), aes(y=prop, x=type, label=total), vjust = 0) +
    facet_wrap(~ sample) +
  labs(fill = "Genotype") + 
  theme_classic() +
  theme(axis.text.x = element_text(size=15),
        axis.title.x = element_blank()) 

ggplot(df_fem, aes(fill=gt, y=prop, x=type)) + 
    geom_bar(position="fill", stat="identity", ) +
    xlab("F2 siblings") + 
    ylab("Proportion of F2 genotype") +
    labs(title = "F2 Females of F1 male (0/0) x female (1/1)") + 
   # geom_text(aes(label = n), position = position_stack(vjust = 0.5)) +
    geom_text(data = filter(df_fem, gt == "1/1"), aes(y=prop, x=type, label=total), vjust = 0) +
    facet_wrap(~ sample) +
  labs(fill = "Genotype") + 
  theme_classic() +
  theme(axis.text.x = element_text(size=15),
        axis.title.x = element_blank()) 

F0 foundress

we expect ?????…

# plot the pooled data for all F0 foundress with male (0/0) x female (1/1) offspring
ggplot(observed, aes(fill=gt, y=prop, x = type)) + 
    geom_bar(position="fill", stat="identity") +
    ylab("Proportion of F0 genotype") +
    labs(title = "F0 foundress of F1 male (0/0) x female (1/1)", 
         subtitle = paste0("Pooled sites = ", observed %>% dplyr::filter(gt == "0/1") %>% summarise(sum(total)) %>% as.numeric())) + 
    theme_classic() +
    theme(axis.text.x = element_blank(), axis.title.x = element_blank()) + 
    labs(fill = "Genotype")

# plot for each individual: 
ggplot(observed, aes(fill=gt, y=prop, x=type)) + 
    geom_bar(position="fill", stat="identity", ) +
    xlab("Observed") + 
    ylab("Proportion of F0 genotype") +
    labs(title = "F0 foundress of F1 male (0/0) x female (1/1)") + 
    geom_text(data = filter(observed, gt == "1/1"), aes(y=prop, x=type, label=total), vjust = 0) +
    facet_wrap(~ sample) +
  labs(fill = "Genotype") + 
  theme_classic() +
  theme(axis.text.x =  element_blank(), axis.title.x = element_blank()) 

(4) F1 male (1/1) x female (0/0)

F2 offspring

We expect zero sites for this cross, since there is no paternal inheritance to the males.
Indeed, there are only 98 sites in the F2 pooled females, and 75 for the F2 pooled males:

# set the families (include only families with at least one adult F2)
family = grep("grndat|grnson",gt$ID, value=TRUE) %>%
  str_extract("[^_]+")  %>%
  unique()

# define a list to put all the data frames in
obs <- list()

for (fam in family) {
 
obs[[fam]] <- table %>%
  dplyr::select(starts_with(fam)) %>%
  dplyr::filter_at(vars(matches("_son")), all_vars(. == "1/1")) %>%
  dplyr::filter_at(vars(matches("_dat")), all_vars(. == "0/0")) %>% 
  dplyr::select(contains("grn")) %>%
  tidyr::pivot_longer(everything())  %>% 
  dplyr::rename(sample = name, gt = value) %>%
  dplyr::count(sample, gt, .drop = FALSE) %>% 
  dplyr::filter(gt %in% c("0/0", "1/1", "0/1")) %>%
  mutate(n = as.numeric(n)) %>%
  group_by(sample) %>%
  mutate(total = as.numeric(sum(n))) %>%
  dplyr::rename(obs = n) %>%
  mutate(sex = case_when(
    grepl("son", sample) ~ "male",
    grepl("dat", sample) ~ "female"))
}

# bind all families together, to a final data frame containing all observed counts
observed <- do.call("rbind", obs)

# state the expected site proportion of each genotype and sex under a perfect mendelian segregation of F2
fem_exp <- data.frame(sex = "female", gt = c("0/0", "0/1", "1/1") , prop = c(0,0,0))
male_exp <- data.frame(sex = "male",gt = c("0/0", "0/1", "1/1"),prop = c(0,0,0))
mendel <- bind_rows(fem_exp,male_exp) %>% mutate(sex = as.character(sex))%>% mutate(gt = as.character(gt)) 
 # mutate(type = "exp")

# using on the total genotypes count of sites per sample, calculate the expected counts, and add them to the observed counts into a final dat table
dat <- observed %>%
  #mutate(sample = as.character(sample)) %>%
  mutate(gt = as.character(gt)) %>%
  left_join(mendel, by = c("gt", "sex")) %>% 
  mutate(exp = total*prop)
x <- dat %>%
  dplyr::select(c("sample","gt","sex", "exp","prop", "total")) %>%
  mutate(type = "exp") %>%
  dplyr::rename(count = exp)

y <- dat %>%
  dplyr::select(c("sample","gt","sex", "obs","total")) %>%
  mutate(type = "obs") %>%
  dplyr::rename(count = obs) %>%
  mutate(prop = (count/total))
  
df_male <- rbind(x, y) %>%
  dplyr::filter(sex == "male")
      
df_fem <- rbind(x, y) %>%
  dplyr::filter(sex == "female")

p_male_pooled <- ggplot(df_male, aes(fill=gt, y=prop, x=type)) + 
    geom_bar(position="fill", stat="identity", ) +
    xlab("F2 males") + 
    ylab("Proportion of F2 genotype") +
    labs(title = "F2 Males of F1 male (1/1) x female (0/0)", 
         subtitle = paste0("Pooled sites = ", df_male %>% dplyr::filter(gt == "0/0") %>% dplyr::filter(type == "obs") %>% summarise(sum(total)) %>% as.numeric())) + 
   scale_x_discrete("type", labels = c("Expected","Observed")) +
labs(fill = "Genotype") + 
  theme_classic() +
  theme(axis.text.x = element_text(size=15),
        axis.title.x = element_blank()) 

p_female_pooled <- ggplot(df_fem, aes(fill=gt, y=prop, x=type)) + 
    geom_bar(position="fill", stat="identity", ) +
    xlab("F2 females") + 
    ylab("Proportion of F2 genotype") +
    labs(title = "F2 females of F1 male (1/1) x female (0/0)", 
         subtitle = paste0("Pooled sites = ", df_fem %>% dplyr::filter(gt == "0/0") %>% dplyr::filter(type == "obs") %>% summarise(sum(total)) %>% as.numeric()))+ 
   scale_x_discrete("type", labels = c("Expected","Observed")) +
labs(fill = "Genotype") + 
  theme_classic() +
  theme(axis.text.x = element_text(size=15),
        axis.title.x = element_blank()) 

grid.arrange(p_female_pooled, p_male_pooled, ncol = 2)

ggplot(df_male, aes(fill=gt, y=prop, x=type)) + 
    geom_bar(position="fill", stat="identity", ) +
    xlab("F2 siblings") + 
    ylab("Proportion of F2 genotype") +
    labs(title = "F2 Males of F1 male (1/1) x female (0/0)") + 
   # geom_text(aes(label = n), position = position_stack(vjust = 0.5)) +
    geom_text(data = filter(df_male, gt == "1/1"), aes(y=prop, x=type, label=total), vjust = 0) +
    facet_wrap(~ sample) +
  labs(fill = "Genotype") + 
  theme_classic() +
  theme(axis.text.x = element_text(size=15),
        axis.title.x = element_blank()) 

ggplot(df_fem, aes(fill=gt, y=prop, x=type)) + 
    geom_bar(position="fill", stat="identity", ) +
    xlab("F2 siblings") + 
    ylab("Proportion of F2 genotype") +
    labs(title = "F2 Females of F1 male (1/1) x female (0/0)") + 
   # geom_text(aes(label = n), position = position_stack(vjust = 0.5)) +
    geom_text(data = filter(df_fem, gt == "1/1"), aes(y=prop, x=type, label=total), vjust = 0) +
    facet_wrap(~ sample) +
  labs(fill = "Genotype") + 
  theme_classic() +
  theme(axis.text.x = element_text(size=15),
        axis.title.x = element_blank()) 

F0 foundress

we expect ?????…

# plot the pooled data for all F0 foundress with male (1/1) x female (0/0) offspring
ggplot(observed, aes(fill=gt, y=prop, x = type)) + 
    geom_bar(position="fill", stat="identity") +
    ylab("Proportion of F0 genotype") +
    labs(title = "F0 foundress of F1 male (1/1) x female (0/0)", 
         subtitle = paste0("Pooled sites = ", observed %>% dplyr::filter(gt == "0/1") %>% summarise(sum(total)) %>% as.numeric())) + 
    theme_classic() +
    theme(axis.text.x = element_blank(), axis.title.x = element_blank()) + 
    labs(fill = "Genotype")

# plot for each individual: 
ggplot(observed, aes(fill=gt, y=prop, x=type)) + 
    geom_bar(position="fill", stat="identity", ) +
    xlab("Observed") + 
    ylab("Proportion of F0 genotype") +
    labs(title = "F0 foundress of F1 male (1/1) x female (0/0)") + 
    geom_text(data = filter(observed, gt == "1/1"), aes(y=prop, x=type, label=total), vjust = 0) +
    facet_wrap(~ sample) +
  labs(fill = "Genotype") + 
  theme_classic() +
  theme(axis.text.x =  element_blank(), axis.title.x = element_blank()) 


Based the ‘Control crossing’ results, I made a vector with all ‘False sites’. These sites will be excluded from the following analyses (5-9).

List of ‘False sites’:

# based on the former control crossing - make a vector of the false sites:


# deduct these sites form the 'table' file

Crossing homozygotic male, with heterozygotic female

(5) F1 male (0/0) x female (0/1)

F2 offspring

We expect F2 females’ sites to be evenly distributed between 0/1 and 0/0 genotypes, and F2 males to be evenly distributed between 0/0 and 1/1 (as they are hemizygote).

# set the families (include only families with at least one adult F2)
family = grep("grndat|grnson",gt$ID, value=TRUE) %>%
  str_extract("[^_]+")  %>%
  unique()

# define a list to put all the data frames in
obs <- list()

for (fam in family) {
 
obs[[fam]] <- table %>%
  dplyr::select(starts_with(fam)) %>%
  dplyr::filter_at(vars(matches("_son")), all_vars(. == "0/0")) %>%
  dplyr::filter_at(vars(matches("_dat")), all_vars(. == "0/1")) %>% 
  dplyr::select(contains("grn")) %>%
  tidyr::pivot_longer(everything())  %>% 
  dplyr::rename(sample = name, gt = value) %>%
  dplyr::count(sample, gt, .drop = FALSE) %>% 
  dplyr::filter(gt %in% c("0/0", "1/1", "0/1")) %>%
  mutate(n = as.numeric(n)) %>%
  group_by(sample) %>%
  mutate(total = as.numeric(sum(n))) %>%
  dplyr::rename(obs = n) %>%
  mutate(sex = case_when(
    grepl("son", sample) ~ "male",
    grepl("dat", sample) ~ "female"))
}

# bind all families together, to a final data frame containing all observed counts
observed <- do.call("rbind", obs)

# state the expected site proportion of each genotype and sex under a perfect mendelian segregation of F2
fem_exp <- data.frame(sex = "female", gt = c("0/0", "0/1", "1/1") , prop = c(0.5, 0.5,0))
male_exp <- data.frame(sex = "male",gt = c("0/0", "0/1", "1/1"),prop = c(0.5, 0,0.5))
mendel <- bind_rows(fem_exp,male_exp) %>% mutate(sex = as.character(sex))%>% mutate(gt = as.character(gt)) 
 # mutate(type = "exp")

# using on the total genotypes count of sites per sample, calculate the expected counts, and add them to the observed counts into a final dat table
dat <- observed %>%
  #mutate(sample = as.character(sample)) %>%
  mutate(gt = as.character(gt)) %>%
  left_join(mendel, by = c("gt", "sex")) %>% 
  mutate(exp = total*prop)
x <- dat %>%
  dplyr::select(c("sample","gt","sex", "exp","prop", "total")) %>%
  mutate(type = "exp") %>%
  dplyr::rename(count = exp)

y <- dat %>%
  dplyr::select(c("sample","gt","sex", "obs","total")) %>%
  mutate(type = "obs") %>%
  dplyr::rename(count = obs) %>%
  mutate(prop = (count/total))
  
df_male <- rbind(x, y) %>%
  dplyr::filter(sex == "male")
      
df_fem <- rbind(x, y) %>%
  dplyr::filter(sex == "female")

p_male_pooled <- ggplot(df_male, aes(fill=gt, y=prop, x=type)) + 
    geom_bar(position="fill", stat="identity", ) +
    xlab("F2 males") + 
    ylab("Proportion of F2 genotype") +
    labs(title = "F2 Males of F1 male (0/0) x female (0/1)", 
         subtitle = paste0("Pooled sites = ", df_male %>% dplyr::filter(gt == "0/0") %>% dplyr::filter(type == "obs") %>% summarise(sum(total)) %>% as.numeric())) + 
   scale_x_discrete("type", labels = c("Expected","Observed")) +
labs(fill = "Genotype") + 
  theme_classic() +
  theme(axis.text.x = element_text(size=15),
        axis.title.x = element_blank()) 

p_female_pooled <- ggplot(df_fem, aes(fill=gt, y=prop, x=type)) + 
    geom_bar(position="fill", stat="identity", ) +
    xlab("F2 females") + 
    ylab("Proportion of F2 genotype") +
    labs(title = "F2 females of F1 male (0/0) x female (0/1)", 
         subtitle = paste0("Pooled sites = ", df_fem %>% dplyr::filter(gt == "0/0") %>% dplyr::filter(type == "obs") %>% summarise(sum(total)) %>% as.numeric())) + 
   scale_x_discrete("type", labels = c("Expected","Observed")) +
labs(fill = "Genotype") + 
  theme_classic() +
  theme(axis.text.x = element_text(size=15),
        axis.title.x = element_blank()) 

grid.arrange(p_female_pooled, p_male_pooled, ncol = 2)

ggplot(df_male, aes(fill=gt, y=prop, x=type)) + 
    geom_bar(position="fill", stat="identity", ) +
    xlab("F2 siblings") + 
    ylab("Proportion of F2 genotype") +
    labs(title = "F2 Males of F1 male (0/0) x female (0/1)") + 
   # geom_text(aes(label = n), position = position_stack(vjust = 0.5)) +
    geom_text(data = filter(df_male, gt == "1/1"), aes(y=prop, x=type, label=total), vjust = 0) +
    facet_wrap(~ sample) +
  labs(fill = "Genotype") +
  theme_classic()

ggplot(df_fem, aes(fill=gt, y=prop, x=type)) + 
    geom_bar(position="fill", stat="identity", ) +
    xlab("F2 siblings") + 
    ylab("Proportion of F2 genotype") +
    labs(title = "F2 Females of F1 male (0/0) x female (0/1)") + 
   # geom_text(aes(label = n), position = position_stack(vjust = 0.5)) +
    geom_text(data = filter(df_fem, gt == "1/1"), aes(y=prop, x=type, label=total), vjust = 0) +
    facet_wrap(~ sample) +
  labs(fill = "Genotype") +
  theme_classic()

F0 foundress

we expect ?????…

# plot the pooled data for all F0 foundress with male (0/0) x female (0/1) offspring
ggplot(observed, aes(fill=gt, y=prop, x = type)) + 
    geom_bar(position="fill", stat="identity") +
    ylab("Proportion of F0 genotype") +
    labs(title = "F0 foundress of F1 male (0/0) x female (0/1)", 
         subtitle = paste0("Pooled sites = ", observed %>% dplyr::filter(gt == "0/1") %>% summarise(sum(total)) %>% as.numeric())) + 
    theme_classic() +
    theme(axis.text.x = element_blank(), axis.title.x = element_blank()) + 
    labs(fill = "Genotype")

# plot for each individual: 
ggplot(observed, aes(fill=gt, y=prop, x=type)) + 
    geom_bar(position="fill", stat="identity", ) +
    xlab("Observed") + 
    ylab("Proportion of F0 genotype") +
    labs(title = "F0 foundress of F1 male (0/0) x female (0/1)") + 
    geom_text(data = filter(observed, gt == "1/1"), aes(y=prop, x=type, label=total), vjust = 0) +
    facet_wrap(~ sample) +
  labs(fill = "Genotype") + 
  theme_classic() +
  theme(axis.text.x =  element_blank(), axis.title.x = element_blank()) 


Based the ‘Control crossing’ results, I made a vector with all ‘False sites’. These sites will be excluded from the following analyses (5-9).

(6) F1 male (1/1) x female (0/1)

F2 offpsring

We expect F2 females’ sites to be evenly distributed between 0/1 and 1/1 genotypes, and F2 males to be evenly distributed between 0/0 and 1/1 (as they are hemizygote).

# set the families (include only families with at least one adult F2)
family = grep("grndat|grnson",gt$ID, value=TRUE) %>%
  str_extract("[^_]+")  %>%
  unique()

# define a list to put all the data frames in
obs <- list()

for (fam in family) {
 
obs[[fam]] <- table %>%
  dplyr::select(starts_with(fam)) %>%
  dplyr::filter_at(vars(matches("_son")), all_vars(. == "1/1")) %>%
  dplyr::filter_at(vars(matches("_dat")), all_vars(. == "0/1")) %>% 
  dplyr::select(contains("grn")) %>%
  tidyr::pivot_longer(everything())  %>% 
  dplyr::rename(sample = name, gt = value) %>%
  dplyr::count(sample, gt, .drop = FALSE) %>% 
  dplyr::filter(gt %in% c("0/0", "1/1", "0/1")) %>%
  mutate(n = as.numeric(n)) %>%
  group_by(sample) %>%
  mutate(total = as.numeric(sum(n))) %>%
  dplyr::rename(obs = n) %>%
  mutate(sex = case_when(
    grepl("son", sample) ~ "male",
    grepl("dat", sample) ~ "female"))
}

# bind all families together, to a final data frame containing all observed counts
observed <- do.call("rbind", obs)

# state the expected site proportion of each genotype and sex under a perfect mendelian segregation of F2
fem_exp <- data.frame(sex = "female", gt = c("0/0", "0/1", "1/1") , prop = c(0,0.5,1))
male_exp <- data.frame(sex = "male",gt = c("0/0", "0/1", "1/1"),prop = c(0.5,0,0.5))
mendel <- bind_rows(fem_exp,male_exp) %>% mutate(sex = as.character(sex))%>% mutate(gt = as.character(gt)) 
 # mutate(type = "exp")

# using on the total genotypes count of sites per sample, calculate the expected counts, and add them to the observed counts into a final dat table
dat <- observed %>%
  #mutate(sample = as.character(sample)) %>%
  mutate(gt = as.character(gt)) %>%
  left_join(mendel, by = c("gt", "sex")) %>% 
  mutate(exp = total*prop)
x <- dat %>%
  dplyr::select(c("sample","gt","sex", "exp","prop", "total")) %>%
  mutate(type = "exp") %>%
  dplyr::rename(count = exp)

y <- dat %>%
  dplyr::select(c("sample","gt","sex", "obs","total")) %>%
  mutate(type = "obs") %>%
  dplyr::rename(count = obs) %>%
  mutate(prop = (count/total))
  
df_male <- rbind(x, y) %>%
  dplyr::filter(sex == "male")
      
df_fem <- rbind(x, y) %>%
  dplyr::filter(sex == "female")

p_male_pooled <- ggplot(df_male, aes(fill=gt, y=prop, x=type)) + 
    geom_bar(position="fill", stat="identity", ) +
    xlab("F2 males") + 
    ylab("Proportion of F2 genotype") +
    labs(title = "F2 Males of F1 male (1/1) x female (0/1)", 
         subtitle = paste0("Pooled sites = ", df_male %>% dplyr::filter(gt == "0/0") %>% dplyr::filter(type == "obs") %>% summarise(sum(total)) %>% as.numeric())) + 
   scale_x_discrete("type", labels = c("Expected","Observed")) +
labs(fill = "Genotype") + 
  theme_classic() +
  theme(axis.text.x = element_text(size=15),
        axis.title.x = element_blank()) 

p_female_pooled <- ggplot(df_fem, aes(fill=gt, y=prop, x=type)) + 
    geom_bar(position="fill", stat="identity", ) +
    xlab("F2 females") + 
    ylab("Proportion of F2 genotype") +
    labs(title = "F2 females of F1 male (1/1) x female (0/1)", 
         subtitle = paste0("Pooled sites = ", df_fem %>% dplyr::filter(gt == "0/0") %>% dplyr::filter(type == "obs") %>% summarise(sum(total)) %>% as.numeric())) + 
   scale_x_discrete("type", labels = c("Expected","Observed")) +
labs(fill = "Genotype") + 
  theme_classic() +
  theme(axis.text.x = element_text(size=15),
        axis.title.x = element_blank()) 

grid.arrange(p_female_pooled, p_male_pooled, ncol = 2)

ggplot(df_male, aes(fill=gt, y=prop, x=type)) + 
    geom_bar(position="fill", stat="identity", ) +
    xlab("F2 siblings") + 
    ylab("Proportion of F2 genotype") +
    labs(title = "F2 Males of F1 male (1/1) x female (0/1)") + 
   # geom_text(aes(label = n), position = position_stack(vjust = 0.5)) +
    geom_text(data = filter(df_male, gt == "1/1"), aes(y=prop, x=type, label=total), vjust = 0) +
    facet_wrap(~ sample) +
  labs(fill = "Genotype") +
  theme_classic()

ggplot(df_fem, aes(fill=gt, y=prop, x=type)) + 
    geom_bar(position="fill", stat="identity", ) +
    xlab("F2 siblings") + 
    ylab("Proportion of F2 genotype") +
    labs(title = "F2 Females of F1 male (1/1) x female (0/1)") + 
   # geom_text(aes(label = n), position = position_stack(vjust = 0.5)) +
    geom_text(data = filter(df_fem, gt == "1/1"), aes(y=prop, x=type, label=total), vjust = 0) +
    facet_wrap(~ sample) +
  labs(fill = "Genotype") +
  theme_classic()

F0 foundress

we expect ?????…

# plot the pooled data for all F0 foundress with male (1/1) x female (0/1) offspring
ggplot(observed, aes(fill=gt, y=prop, x = type)) + 
    geom_bar(position="fill", stat="identity") +
    ylab("Proportion of F0 genotype") +
    labs(title = "F0 foundress of F1 male (1/1) x female (0/1)", 
         subtitle = paste0("Pooled sites = ", observed %>% dplyr::filter(gt == "0/1") %>% summarise(sum(total)) %>% as.numeric())) + 
    theme_classic() +
    theme(axis.text.x = element_blank(), axis.title.x = element_blank()) + 
    labs(fill = "Genotype")

# plot for each individual: 
ggplot(observed, aes(fill=gt, y=prop, x=type)) + 
    geom_bar(position="fill", stat="identity", ) +
    xlab("Observed") + 
    ylab("Proportion of F0 genotype") +
    labs(title = "F0 foundress of F1 male (1/1) x female (0/1)") + 
    geom_text(data = filter(observed, gt == "1/1"), aes(y=prop, x=type, label=total), vjust = 0) +
    facet_wrap(~ sample) +
  labs(fill = "Genotype") + 
  theme_classic() +
  theme(axis.text.x =  element_blank(), axis.title.x = element_blank()) 


Based the ‘Control crossing’ results, I made a vector with all ‘False sites’. These sites will be excluded from the following analyses (5-9).

Crossing heterozygotic male, with homozygotic female

The former crosses of heterozygotic females (5 and 6) show that F2 males can be heterozygotic and carry two alleles, like their mother.
But are these sites real? and if they are, can they be transmitted to their daughters?
Looking at the F2 offspring of crossing heterozygotic F1 males with females, we tested 3 possible predictions:

  1. Heterozygotic sites are false: they are actually hemizygotic (0)
  2. Heterozygotic sites are false: they are actually hemizygotic (1)
  3. Heterozygotic sites are real: Males can carry and inherit two alleles (0/1)

(7) F1 male (0/1) x female (0/0)

F2 offspring

# set the families (include only families with at least one adult F2)
family = grep("grndat|grnson",gt$ID, value=TRUE) %>%
  str_extract("[^_]+")  %>%
  unique()

# define a list to put all the data frames in
obs <- list()

for (fam in family) {
 
obs[[fam]] <- table %>%
  dplyr::select(starts_with(fam)) %>%
  dplyr::filter_at(vars(matches("_son")), all_vars(. == "0/1")) %>%
  dplyr::filter_at(vars(matches("_dat")), all_vars(. == "0/0")) %>% 
  dplyr::select(contains("grn")) %>%
  tidyr::pivot_longer(everything())  %>% 
  dplyr::rename(sample = name, gt = value) %>%
  dplyr::count(sample, gt, .drop = FALSE) %>% 
  dplyr::filter(gt %in% c("0/0", "1/1", "0/1")) %>%
  mutate(n = as.numeric(n)) %>%
  group_by(sample) %>%
  mutate(total = as.numeric(sum(n))) %>%
  dplyr::rename(obs = n) %>%
  mutate(sex = case_when(
    grepl("son", sample) ~ "male",
    grepl("dat", sample) ~ "female"))
}

# bind all families together, to a final data frame containing all observed counts
observed <- do.call("rbind", obs) 

# state the expected site proportion of each genotype and sex under the different hypothesis:
# Heterozygotic sites are real:Males can carry and inherit two alleles
fem_sitesReal <- data.frame(type = "sitesReal", sex = "female", gt = c("0/0", "0/1", "1/1") , prop = c(0.5, 0.5, 0))

# Heterozygotic sites are false:Males can carry and inherit only one allele
# the heterozygotic sites are actually (0)
fem_sitesFalse_0 <- data.frame(type = "sitesFalse_0",sex = "female", gt = c("0/0", "0/1", "1/1") , prop = c(1,0,0))
# the heterozygotic sites are actually (1)
fem_sitesFalse_1 <- data.frame(type = "sitesFalse_1", sex = "female", gt = c("0/0", "0/1", "1/1") , prop = c(0,1,0))

# males F2 genotype is expected to stay the same , as it is produced from unfertilized egg, hence have no paternal inheritance
male_sitesReal <- data.frame(type = "sitesReal",sex = "male",gt = c("0/0", "0/1", "1/1"),prop = c(1,0,0))
male_sitesFalse_0 <- data.frame(type = "sitesFalse_0",sex = "male",gt = c("0/0", "0/1", "1/1"),prop = c(1,0,0))
male_sitesFalse_1 <- data.frame(type = "sitesFalse_1", sex = "male",gt = c("0/0", "0/1", "1/1"),prop = c(1,0,0))

expected_prop <- bind_rows(fem_sitesReal,fem_sitesFalse_0,fem_sitesFalse_1,male_sitesReal,male_sitesFalse_0,male_sitesFalse_1) %>% mutate(sex = as.character(sex))%>% mutate(gt = as.character(gt)) 
 # mutate(type = "exp")

# using on the total genotypes count of sites per sample, calculate the expected counts, and add them to the observed counts into a final dat table
expected <- observed %>%
  mutate(gt = as.character(gt)) %>%
  full_join(expected_prop, by = c("gt", "sex")) %>% 
  mutate(count = total*prop)

x <- observed %>%
  dplyr::select(c("sample","sex","gt","total", "obs")) %>%
  mutate(type = "Observed") %>%
  dplyr::rename(count = obs) %>%
  mutate(prop = (count/total))

dat <-  bind_rows(x, expected) %>% 
  dplyr::select(-obs)
# subset for males and females, and reorder the "types" levels:
dat_fem <- dat %>% dplyr::filter(sex == "female")
level_order_fem <- factor(dat_fem$type, level = c("sitesFalse_0", "sitesFalse_1", "sitesReal", "Observed"))

dat_male <- dat %>% dplyr::filter(sex == "male")
level_order_male <- factor(dat_male$type, level = c("sitesFalse_0", "sitesFalse_1", "sitesReal", "Observed"))

p_male_pooled <- ggplot(dat_male, aes(fill=gt, y=prop, x=level_order_male)) + 
    geom_bar(position="fill", stat="identity", width = 0.8) +
    ylab("Proportion of F2 genotype") +
    labs(title = paste0("F2 Male, pooled sites = ", dat_male %>% dplyr::filter(type == "Observed") %>% summarise(sum(count)) %>% as.numeric()))  + 
  labs(fill = "Genotype:") + 
  theme_classic() +
  theme(axis.text.x = element_text(size=12, angle = 45,hjust = 1),
        axis.title.x = element_blank())

p_female_pooled <- ggplot(dat_fem, aes(fill=gt, y=prop, x=level_order_fem)) + 
    geom_bar(position="fill", stat="identity", width = 0.8) +
    ylab("Proportion of F2 genotype") +
    labs(title = paste0("F2 Female, pooled sites = ", dat_fem %>% dplyr::filter(type == "Observed") %>% summarise(sum(count)) %>% as.numeric())) + 
  labs(fill = "Genotype:") +
  theme_classic() +
  theme(axis.text.x = element_text(size=12, angle = 45,hjust = 1),
        axis.title.x = element_blank())

combined <- p_female_pooled + p_male_pooled & theme(legend.position = "bottom")
combined + plot_layout(guides = "collect") + plot_annotation(
  title = 'F2 genotype. F1 Cross: male (0/1) x female (0/0)')

# plot each individual:
ggplot(dat_male, aes(fill=gt, y=prop, x=level_order_male)) + 
    geom_bar(position="fill", stat="identity", ) +
    ylab("Proportion of F2 genotype") +
    labs(title = "F2 male genotype. F1 Cross: male (0/1) x female (0/0)", fill = "Genotype:") +
    geom_text(data = dat_male %>% filter(gt == "1/1") %>% filter(type == "Observed"), aes(x = 4, y = 0.8, label = total), inherit.aes = FALSE) +
    facet_wrap(~ sample) +
    theme_classic()+ 
    theme(legend.position = "bottom", axis.title.x = element_blank(), axis.text.x = element_text(size=8, angle = 45,hjust = 1)) 

ggplot(dat_fem, aes(fill=gt, y=prop, x=level_order_fem)) + 
    geom_bar(position="fill", stat="identity", ) +
    ylab("Proportion of F2 genotype") +
    labs(title = "F2 female genotype. F1 Cross: male (0/1) x female (0/0)", fill = "Genotype:") + 
    geom_text(data = dat_fem %>% filter(gt == "1/1") %>% filter(type == "Observed"), aes(x = 4, y = 0.8, label = total), inherit.aes = FALSE) +
    facet_wrap(~ sample) +
    theme_classic()+ 
    theme(legend.position = "bottom", axis.title.x = element_blank(), axis.text.x = element_text(size=8, angle = 45,hjust = 1)) 

F0 foundress

we expect ?????…

# plot the pooled data for all F0 foundress with male (0/1) x female (0/0) offspring
ggplot(observed, aes(fill=gt, y=prop, x = type)) + 
    geom_bar(position="fill", stat="identity") +
    ylab("Proportion of F0 genotype") +
    labs(title = "F0 foundress of F1 male (0/1) x female (0/0)", 
         subtitle = paste0("Pooled sites = ", observed %>% dplyr::filter(gt == "0/1") %>% summarise(sum(total)) %>% as.numeric())) + 
    theme_classic() +
    theme(axis.text.x = element_blank(), axis.title.x = element_blank()) + 
    labs(fill = "Genotype")

# plot for each individual: 
ggplot(observed, aes(fill=gt, y=prop, x=type)) + 
    geom_bar(position="fill", stat="identity", ) +
    xlab("Observed") + 
    ylab("Proportion of F0 genotype") +
    labs(title = "F0 foundress of F1 male (0/1) x female (0/0)") + 
    geom_text(data = filter(observed, gt == "1/1"), aes(y=prop, x=type, label=total), vjust = 0) +
    facet_wrap(~ sample) +
  labs(fill = "Genotype") + 
  theme_classic() +
  theme(axis.text.x =  element_blank(), axis.title.x = element_blank()) 


Based the ‘Control crossing’ results, I made a vector with all ‘False sites’. These sites will be excluded from the following analyses (5-9).

(8) F1 male (0/1) x female (1/1)

F2 offspring

# set the families (include only families with at least one adult F2)
family = grep("grndat|grnson",gt$ID, value=TRUE) %>%
  str_extract("[^_]+")  %>%
  unique()

# define a list to put all the data frames in
obs <- list()

for (fam in family) {
 
obs[[fam]] <- table %>%
  dplyr::select(starts_with(fam)) %>%
  dplyr::filter_at(vars(matches("_son")), all_vars(. == "0/1")) %>%
  dplyr::filter_at(vars(matches("_dat")), all_vars(. == "1/1")) %>% 
  dplyr::select(contains("grn")) %>%
  tidyr::pivot_longer(everything())  %>% 
  dplyr::rename(sample = name, gt = value) %>%
  dplyr::count(sample, gt, .drop = FALSE) %>% 
  dplyr::filter(gt %in% c("0/0", "1/1", "0/1")) %>%
  mutate(n = as.numeric(n)) %>%
  group_by(sample) %>%
  mutate(total = as.numeric(sum(n))) %>%
  dplyr::rename(obs = n) %>%
  mutate(sex = case_when(
    grepl("son", sample) ~ "male",
    grepl("dat", sample) ~ "female"))
}

# bind all families together, to a final data frame containing all observed counts
observed <- do.call("rbind", obs) 

# state the expected site proportion of each genotype and sex under the different hypothesis:
# Heterozygotic sites are real:Males can carry and inherit two alleles
fem_sitesReal <- data.frame(type = "sitesReal", sex = "female", gt = c("0/0", "0/1", "1/1") , prop = c(0,0.5,0.5))

# Heterozygotic sites are false:Males can carry and inherit only one allele
# the heterozygotic sites are actually (0)
fem_sitesFalse_0 <- data.frame(type = "sitesFalse_0",sex = "female", gt = c("0/0", "0/1", "1/1") , prop = c(0,1,0))
# the heterozygotic sites are actually (1)
fem_sitesFalse_1 <- data.frame(type = "sitesFalse_1", sex = "female", gt = c("0/0", "0/1", "1/1") , prop = c(0,0,1))

# males F2 genotype is expected to stay the same , as it is produced from unfertilized egg, hence have no paternal inheritance
male_sitesReal <- data.frame(type = "sitesReal",sex = "male",gt = c("0/0", "0/1", "1/1"),prop = c(0,0,1))
male_sitesFalse_0 <- data.frame(type = "sitesFalse_0",sex = "male",gt = c("0/0", "0/1", "1/1"),prop = c(0,0,1))
male_sitesFalse_1 <- data.frame(type = "sitesFalse_1", sex = "male",gt = c("0/0", "0/1", "1/1"),prop = c(0,0,1))

expected_prop <- bind_rows(fem_sitesReal,fem_sitesFalse_0,fem_sitesFalse_1,male_sitesReal,male_sitesFalse_0,male_sitesFalse_1) %>% mutate(sex = as.character(sex))%>% mutate(gt = as.character(gt)) 
 # mutate(type = "exp")

# using on the total genotypes count of sites per sample, calculate the expected counts, and add them to the observed counts into a final dat table
expected <- observed %>%
  mutate(gt = as.character(gt)) %>%
  full_join(expected_prop, by = c("gt", "sex")) %>% 
  mutate(count = total*prop)

x <- observed %>%
  dplyr::select(c("sample","sex","gt","total", "obs")) %>%
  mutate(type = "Observed") %>%
  dplyr::rename(count = obs) %>%
  mutate(prop = (count/total))

dat <-  bind_rows(x, expected) %>% 
  dplyr::select(-obs)
# subset for males and females, and reorder the "types" levels:
dat_fem <- dat %>% dplyr::filter(sex == "female")
level_order_fem <- factor(dat_fem$type, level = c("sitesFalse_0", "sitesFalse_1", "sitesReal", "Observed"))

dat_male <- dat %>% dplyr::filter(sex == "male")
level_order_male <- factor(dat_male$type, level = c("sitesFalse_0", "sitesFalse_1", "sitesReal", "Observed"))

p_male_pooled <- ggplot(dat_male, aes(fill=gt, y=prop, x=level_order_male)) + 
    geom_bar(position="fill", stat="identity", width = 0.8) +
    ylab("Proportion of F2 genotype") +
    labs(title = paste0("F2 Male, pooled sites = ", dat_male %>% dplyr::filter(type == "Observed") %>% summarise(sum(count)) %>% as.numeric()))  + 
  labs(fill = "Genotype:") + 
  theme_classic() +
  theme(axis.text.x = element_text(size=12, angle = 45,hjust = 1),
        axis.title.x = element_blank())

p_female_pooled <- ggplot(dat_fem, aes(fill=gt, y=prop, x=level_order_fem)) + 
    geom_bar(position="fill", stat="identity", width = 0.8) +
    ylab("Proportion of F2 genotype") +
    labs(title = paste0("F2 Female, pooled sites = ", dat_fem %>% dplyr::filter(type == "Observed") %>% summarise(sum(count)) %>% as.numeric())) + 
  labs(fill = "Genotype:") +
  theme_classic() +
  theme(axis.text.x = element_text(size=12, angle = 45,hjust = 1),
        axis.title.x = element_blank())

combined <- p_female_pooled + p_male_pooled & theme(legend.position = "bottom")
combined + plot_layout(guides = "collect") + plot_annotation(
  title = 'F2 genotype. F1 Cross: male (0/1) x female (1/1)')

# plot each individual:
ggplot(dat_male, aes(fill=gt, y=prop, x=level_order_male)) + 
    geom_bar(position="fill", stat="identity", ) +
    ylab("Proportion of F2 genotype") +
    labs(title = "F2 male genotype. F1 Cross: male (0/1) x female (1/1)", fill = "Genotype:") +
    geom_text(data = dat_male %>% filter(gt == "1/1") %>% filter(type == "Observed"), aes(x = 4, y = 0.8, label = total), inherit.aes = FALSE) +
    facet_wrap(~ sample) +
    theme_classic()+ 
    theme(legend.position = "bottom", axis.title.x = element_blank(), axis.text.x = element_text(size=8, angle = 45,hjust = 1)) 

ggplot(dat_fem, aes(fill=gt, y=prop, x=level_order_fem)) + 
    geom_bar(position="fill", stat="identity", ) +
    ylab("Proportion of F2 genotype") +
    labs(title = "F2 female genotype. F1 Cross: male (0/1) x female (1/1)", fill = "Genotype:") + 
    geom_text(data = dat_fem %>% filter(gt == "1/1") %>% filter(type == "Observed"), aes(x = 4, y = 0.8, label = total), inherit.aes = FALSE) +
    facet_wrap(~ sample) +
    theme_classic()+ 
    theme(legend.position = "bottom", axis.title.x = element_blank(), axis.text.x = element_text(size=8, angle = 45,hjust = 1)) 

F0 foundress

we expect ?????…

# plot the pooled data for all F0 foundress with male (0/1) x female (1/1) offspring
ggplot(observed, aes(fill=gt, y=prop, x = type)) + 
    geom_bar(position="fill", stat="identity") +
    ylab("Proportion of F0 genotype") +
    labs(title = "F0 foundress of F1 male (0/1) x female (1/1)", 
         subtitle = paste0("Pooled sites = ", observed %>% dplyr::filter(gt == "0/1") %>% summarise(sum(total)) %>% as.numeric())) + 
    theme_classic() +
    theme(axis.text.x = element_blank(), axis.title.x = element_blank()) + 
    labs(fill = "Genotype")

# plot for each individual: 
ggplot(observed, aes(fill=gt, y=prop, x=type)) + 
    geom_bar(position="fill", stat="identity", ) +
    xlab("Observed") + 
    ylab("Proportion of F0 genotype") +
    labs(title = "F0 foundress of F1 male (0/1) x female (1/1)") + 
    geom_text(data = filter(observed, gt == "1/1"), aes(y=prop, x=type, label=total), vjust = 0) +
    facet_wrap(~ sample) +
  labs(fill = "Genotype") + 
  theme_classic() +
  theme(axis.text.x =  element_blank(), axis.title.x = element_blank()) 


Based the ‘Control crossing’ results, I made a vector with all ‘False sites’. These sites will be excluded from the following analyses (5-9).

(9) F1 male (0/1) x female (0/1)

F2 offspring

# set the families (include only families with at least one adult F2)
family = grep("grndat|grnson",gt$ID, value=TRUE) %>%
  str_extract("[^_]+")  %>%
  unique()

# define a list to put all the data frames in
obs <- list()

for (fam in family) {
 
obs[[fam]] <- table %>%
  dplyr::select(starts_with(fam)) %>%
  dplyr::filter_at(vars(matches("_son")), all_vars(. == "0/1")) %>%
  dplyr::filter_at(vars(matches("_dat")), all_vars(. == "0/1")) %>% 
  dplyr::select(contains("grn")) %>%
  tidyr::pivot_longer(everything())  %>% 
  dplyr::rename(sample = name, gt = value) %>%
  dplyr::count(sample, gt, .drop = FALSE) %>% 
  dplyr::filter(gt %in% c("0/0", "1/1", "0/1")) %>%
  mutate(n = as.numeric(n)) %>%
  group_by(sample) %>%
  mutate(total = as.numeric(sum(n))) %>%
  dplyr::rename(obs = n) %>%
  mutate(sex = case_when(
    grepl("son", sample) ~ "male",
    grepl("dat", sample) ~ "female"))
}

# bind all families together, to a final data frame containing all observed counts
observed <- do.call("rbind", obs) 

# state the expected site proportion of each genotype and sex under the different hypothesis:
# Heterozygotic sites are real:Males can carry and inherit two alleles
fem_sitesReal <- data.frame(type = "sitesReal", sex = "female", gt = c("0/0", "0/1", "1/1") , prop = c(0.25, 0.5, 0.25))

# Heterozygotic sites are false:Males can carry and inherit only one allele
# the heterozygotic sites are actually (0)
fem_sitesFalse_0 <- data.frame(type = "sitesFalse_0",sex = "female", gt = c("0/0", "0/1", "1/1") , prop = c(0.5, 0.5, 0))
# the heterozygotic sites are actually (1)
fem_sitesFalse_1 <- data.frame(type = "sitesFalse_1", sex = "female", gt = c("0/0", "0/1", "1/1") , prop = c(0, 0.5, 0.5))

# males F2 genotype is expected to stay the same , as it is produced from unfertilized egg, hence have no paternal inheritance
male_sitesReal <- data.frame(type = "sitesReal",sex = "male",gt = c("0/0", "0/1", "1/1"),prop = c(0.5, 0, 0.5))
male_sitesFalse_0 <- data.frame(type = "sitesFalse_0",sex = "male",gt = c("0/0", "0/1", "1/1"),prop = c(0.5, 0, 0.5))
male_sitesFalse_1 <- data.frame(type = "sitesFalse_1", sex = "male",gt = c("0/0", "0/1", "1/1"),prop = c(0.5, 0, 0.5))

expected_prop <- bind_rows(fem_sitesReal,fem_sitesFalse_0,fem_sitesFalse_1,male_sitesReal,male_sitesFalse_0,male_sitesFalse_1) %>% mutate(sex = as.character(sex))%>% mutate(gt = as.character(gt)) 
 # mutate(type = "exp")

# using on the total genotypes count of sites per sample, calculate the expected counts, and add them to the observed counts into a final dat table
expected <- observed %>%
  mutate(gt = as.character(gt)) %>%
  full_join(expected_prop, by = c("gt", "sex")) %>% 
  mutate(count = total*prop)

x <- observed %>%
  dplyr::select(c("sample","sex","gt","total", "obs")) %>%
  mutate(type = "Observed") %>%
  dplyr::rename(count = obs) %>%
  mutate(prop = (count/total))

dat <-  bind_rows(x, expected) %>% 
  dplyr::select(-obs)
# subset for males and females, and reorder the "types" levels:
dat_fem <- dat %>% dplyr::filter(sex == "female")
level_order_fem <- factor(dat_fem$type, level = c("sitesFalse_0", "sitesFalse_1", "sitesReal", "Observed"))

dat_male <- dat %>% dplyr::filter(sex == "male")
level_order_male <- factor(dat_male$type, level = c("sitesFalse_0", "sitesFalse_1", "sitesReal", "Observed"))

p_male_pooled <- ggplot(dat_male, aes(fill=gt, y=prop, x=level_order_male)) + 
    geom_bar(position="fill", stat="identity", width = 0.8) +
    ylab("Proportion of F2 genotype") +
    labs(title = paste0("F2 Male, pooled sites = ", dat_male %>% dplyr::filter(type == "Observed") %>% summarise(sum(count)) %>% as.numeric()))  + 
  labs(fill = "Genotype:") + 
  theme_classic() +
  theme(axis.text.x = element_text(size=12, angle = 45,hjust = 1),
        axis.title.x = element_blank())

p_female_pooled <- ggplot(dat_fem, aes(fill=gt, y=prop, x=level_order_fem)) + 
    geom_bar(position="fill", stat="identity", width = 0.8) +
    ylab("Proportion of F2 genotype") +
    labs(title = paste0("F2 Female, pooled sites = ", dat_fem %>% dplyr::filter(type == "Observed") %>% summarise(sum(count)) %>% as.numeric())) + 
  labs(fill = "Genotype:") +
  theme_classic() +
  theme(axis.text.x = element_text(size=12, angle = 45,hjust = 1),
        axis.title.x = element_blank())

combined <- p_female_pooled + p_male_pooled & theme(legend.position = "bottom")
combined + plot_layout(guides = "collect") + plot_annotation(
  title = 'F2 genotype. F1 Cross: male (0/1) x female (0/1)')

# plot each individual:
ggplot(dat_male, aes(fill=gt, y=prop, x=level_order_male)) + 
    geom_bar(position="fill", stat="identity", ) +
    ylab("Proportion of F2 genotype") +
    labs(title = "F2 male genotype. F1 Cross: male (0/1) x female (0/1)", fill = "Genotype:") +
    geom_text(data = dat_male %>% filter(gt == "1/1") %>% filter(type == "Observed"), aes(x = 4, y = 0.8, label = total), inherit.aes = FALSE) +
    facet_wrap(~ sample) +
    theme_classic()+ 
    theme(legend.position = "bottom", axis.title.x = element_blank(), axis.text.x = element_text(size=8, angle = 45,hjust = 1)) 

ggplot(dat_fem, aes(fill=gt, y=prop, x=level_order_fem)) + 
    geom_bar(position="fill", stat="identity", ) +
    ylab("Proportion of F2 genotype") +
    labs(title = "F2 female genotype. F1 Cross: male (0/1) x female (0/1)", fill = "Genotype:") + 
    geom_text(data = dat_fem %>% filter(gt == "1/1") %>% filter(type == "Observed"), aes(x = 4, y = 0.8, label = total), inherit.aes = FALSE) +
    facet_wrap(~ sample) +
    theme_classic()+ 
    theme(legend.position = "bottom", axis.title.x = element_blank(), axis.text.x = element_text(size=8, angle = 45,hjust = 1)) 

F0 foundress

we expect ?????…

# plot the pooled data for all F0 foundress with male (0/1) x female (1/1) offspring
ggplot(observed, aes(fill=gt, y=prop, x = type)) + 
    geom_bar(position="fill", stat="identity") +
    ylab("Proportion of F0 genotype") +
    labs(title = "F0 foundress of F1 male (0/1) x female (0/1)", 
         subtitle = paste0("Pooled sites = ", observed %>% dplyr::filter(gt == "1/1") %>% summarise(sum(total)) %>% as.numeric())) + 
    theme_classic() +
    theme(axis.text.x = element_blank(), axis.title.x = element_blank()) + 
    labs(fill = "Genotype")

# plot for each individual: 
ggplot(observed, aes(fill=gt, y=prop, x=type)) + 
    geom_bar(position="fill", stat="identity", ) +
    xlab("Observed") + 
    ylab("Proportion of F0 genotype") +
    labs(title = "F0 foundress of F1 male (0/1) x female (0/1)") + 
    geom_text(data = filter(observed, gt == "1/1"), aes(y=prop, x=type, label=total), vjust = 0) +
    facet_wrap(~ sample) +
  labs(fill = "Genotype") + 
  theme_classic() +
  theme(axis.text.x =  element_blank(), axis.title.x = element_blank()) 


Statistics

We did 2 types of a-parametric tests for independence, testing the predicted counts against the observed counts:

  • Cochran–Mantel–Haenszel Test (CMH) for Repeated Tests of Independence
  • Replicated G-test of goodness of fit

Cochran–Mantel–Haenszel Test for Repeated Tests of Independence

CMH Test for females

CMH Test for males

Replicated G-test of goodness of fit

רק לדוגמא

for one family, 284

for one family, 338

for one family, 110

all genotypes for the 3 families

old code for false sites

# set the families (include only families with at least one adult F2)
family = grep("grndat|grnson",gt$ID, value=TRUE) %>%
  str_extract("[^_]+")  %>%
  unique()

# define a list to put all the data frames in
list_grndat <- list()

for (fam in family) {
 
list_grndat[[fam]] <- table %>%
  dplyr::select(starts_with(fam)) %>%
  dplyr::filter_at(vars(matches("_son")), all_vars(. == "0/0")) %>%
  dplyr::filter_at(vars(matches("_dat")), all_vars(. == "0/0")) %>% 
  dplyr::select(contains("grndat")) %>%
  rownames_to_column("site") 
}

#%>% dplyr::filter(across(everything(), ~ !str_detect(., "0/0")))%>% na.omit()
  
list_grndat$`240` %>% view()

# for male F2
list_grnson <- list()

for (fam in family) {
 
list_grnson[[fam]] <- table %>%
  dplyr::select(starts_with(fam)) %>%
  dplyr::filter_at(vars(matches("_son")), all_vars(. == "0/0")) %>%
  dplyr::filter_at(vars(matches("_dat")), all_vars(. == "0/0")) %>% 
  dplyr::select(contains("grnson")) %>%
  rownames_to_column("site") 
}
#%>% dplyr::filter(across(everything(), ~ !str_detect(., "0/0"))) %>% na.omit()
list_grnson$`600`
##                         site 600_601a_grnson
## 1     NW_019211454.1_4912263             0/0
## 2     NW_019211454.1_4920532             0/0
## 3     NW_019211454.1_4921004             0/0
## 4     NW_019211454.1_4952477             0/0
## 5     NW_019211454.1_4966874            <NA>
## 6     NW_019211454.1_4975120            <NA>
## 7     NW_019211454.1_4978792             0/0
## 8     NW_019211454.1_4980011            <NA>
## 9     NW_019211454.1_5013515             0/0
## 10    NW_019211454.1_5021396             0/0
## 11    NW_019211454.1_5095835             0/0
## 12    NW_019211454.1_5122236            <NA>
## 13    NW_019211454.1_5166671            <NA>
## 14    NW_019211454.1_5175616            <NA>
## 15    NW_019211454.1_5184381             0/0
## 16    NW_019211454.1_5214324             0/0
## 17    NW_019211454.1_5229776            <NA>
## 18    NW_019211454.1_5257386             0/0
## 19    NW_019211454.1_5261860            <NA>
## 20    NW_019211454.1_5266356             0/0
## 21    NW_019211454.1_5268357             0/0
## 22    NW_019211454.1_5279110             0/0
## 23    NW_019211454.1_5918916             0/0
## 24    NW_019211454.1_6033069             0/0
## 25    NW_019211454.1_6041909             0/0
## 26    NW_019211454.1_6167320            <NA>
## 27    NW_019211454.1_6226168             0/0
## 28    NW_019211454.1_6251092             0/0
## 29    NW_019211454.1_6259775            <NA>
## 30    NW_019211454.1_6366438            <NA>
## 31    NW_019211454.1_6453874            <NA>
## 32    NW_019211454.1_6479391             0/0
## 33    NW_019211454.1_6503569             0/0
## 34    NW_019211454.1_6509668            <NA>
## 35    NW_019211454.1_6512550             0/0
## 36    NW_019211454.1_6512723            <NA>
## 37    NW_019211454.1_6516631            <NA>
## 38    NW_019211454.1_6538442            <NA>
## 39    NW_019211454.1_8203978            <NA>
## 40    NW_019211454.1_9513244            <NA>
## 41    NW_019211454.1_9526241             0/0
## 42    NW_019211454.1_9526415             0/0
## 43   NW_019211454.1_10787954             0/0
## 44   NW_019211454.1_10789989             0/0
## 45   NW_019211454.1_12645630             0/0
## 46   NW_019211454.1_12703153             0/0
## 47   NW_019211454.1_12852621            <NA>
## 48   NW_019211454.1_12860342            <NA>
## 49   NW_019211454.1_12869979            <NA>
## 50   NW_019211454.1_12874184            <NA>
## 51   NW_019211454.1_12889655            <NA>
## 52   NW_019211454.1_12943805            <NA>
## 53   NW_019211454.1_13005638             0/0
## 54   NW_019211454.1_13010879             0/0
## 55   NW_019211454.1_13013741            <NA>
## 56   NW_019211454.1_13044739            <NA>
## 57   NW_019211454.1_13056888            <NA>
## 58   NW_019211454.1_13063554             0/0
## 59   NW_019211454.1_13073686            <NA>
## 60   NW_019211454.1_13100815            <NA>
## 61   NW_019211454.1_13107115             0/0
## 62   NW_019211454.1_13146372            <NA>
## 63   NW_019211454.1_13160204            <NA>
## 64   NW_019211454.1_13177229            <NA>
## 65   NW_019211454.1_13186306            <NA>
## 66   NW_019211454.1_13204249             0/0
## 67   NW_019211454.1_14865144             0/0
## 68   NW_019211454.1_14955656             0/0
## 69   NW_019211454.1_14999663            <NA>
## 70   NW_019211454.1_14999711            <NA>
## 71   NW_019211454.1_15001520             0/0
## 72   NW_019211454.1_15003776            <NA>
## 73   NW_019211454.1_15011010            <NA>
## 74   NW_019211454.1_15037118            <NA>
## 75   NW_019211454.1_15153923            <NA>
## 76   NW_019211454.1_15158713             0/0
## 77   NW_019211454.1_15161035             0/0
## 78   NW_019211454.1_15168967             0/0
## 79   NW_019211454.1_15181015             0/0
## 80   NW_019211454.1_15610546            <NA>
## 81   NW_019211454.1_15648222             0/0
## 82   NW_019211454.1_15648701             0/0
## 83   NW_019211454.1_15701177            <NA>
## 84   NW_019211454.1_15712292             0/0
## 85   NW_019211454.1_16244152             0/0
## 86   NW_019211454.1_16246384             0/0
## 87   NW_019211454.1_16252347            <NA>
## 88   NW_019211454.1_16252834             0/0
## 89   NW_019211454.1_16266943             0/0
## 90   NW_019211454.1_16297722            <NA>
## 91   NW_019211454.1_16312031             0/0
## 92   NW_019211454.1_16312209             0/0
## 93   NW_019211454.1_16340860            <NA>
## 94   NW_019211454.1_16348147            <NA>
## 95   NW_019211454.1_16350959             0/0
## 96   NW_019211454.1_16353460            <NA>
## 97   NW_019211454.1_16357130             0/0
## 98   NW_019211454.1_16368943            <NA>
## 99   NW_019211454.1_16380802            <NA>
## 100  NW_019211454.1_16387023             0/0
## 101  NW_019211454.1_16389486             0/0
## 102  NW_019211454.1_16410682            <NA>
## 103  NW_019211454.1_16419671             0/0
## 104  NW_019211454.1_16461135            <NA>
## 105  NW_019211454.1_16469261             0/0
## 106  NW_019211454.1_16472031             0/0
## 107  NW_019211454.1_16483028             0/0
## 108  NW_019211454.1_16502287             0/0
## 109  NW_019211454.1_16519175            <NA>
## 110  NW_019211454.1_16548057             0/0
## 111  NW_019211454.1_16548075             0/0
## 112  NW_019211454.1_16548390             0/0
## 113  NW_019211454.1_16548759             0/0
## 114  NW_019211454.1_16549009             0/0
## 115  NW_019211454.1_16550594            <NA>
## 116  NW_019211454.1_16555740            <NA>
## 117  NW_019211454.1_16557116             0/0
## 118  NW_019211454.1_16557423             0/0
## 119  NW_019211454.1_16557649             0/0
## 120  NW_019211454.1_16561544            <NA>
## 121  NW_019211454.1_16574668            <NA>
## 122  NW_019211454.1_16580293            <NA>
## 123  NW_019211454.1_16645606             0/0
## 124  NW_019211454.1_16647761            <NA>
## 125  NW_019211454.1_16647867             0/0
## 126  NW_019211454.1_16655076             0/0
## 127  NW_019211454.1_16656654            <NA>
## 128  NW_019211454.1_16668875             0/0
## 129  NW_019211454.1_16679234            <NA>
## 130  NW_019211454.1_16689923             0/0
## 131  NW_019211454.1_16693222             0/0
## 132  NW_019211454.1_16698054             0/0
## 133  NW_019211454.1_16969415             0/0
## 134  NW_019211454.1_16978613             0/0
## 135  NW_019211454.1_16995141             0/0
## 136  NW_019211454.1_17114649            <NA>
## 137  NW_019211454.1_17129653            <NA>
## 138  NW_019211454.1_17146268            <NA>
## 139  NW_019211454.1_17148834            <NA>
## 140  NW_019211454.1_17182202            <NA>
## 141  NW_019211454.1_17193528            <NA>
## 142  NW_019211454.1_17198346             0/0
## 143  NW_019211454.1_17198496            <NA>
## 144  NW_019211454.1_17211037             0/0
## 145  NW_019211454.1_17217149            <NA>
## 146  NW_019211454.1_17219103             0/0
## 147  NW_019211454.1_17223506             0/0
## 148  NW_019211454.1_17243529             0/0
## 149  NW_019211454.1_17246645            <NA>
## 150  NW_019211454.1_17248571             0/0
## 151  NW_019211454.1_17257077             0/0
## 152  NW_019211454.1_17260462             0/0
## 153  NW_019211454.1_17271714             0/0
## 154  NW_019211454.1_17290015            <NA>
## 155  NW_019211454.1_17291980            <NA>
## 156  NW_019211454.1_17301347             0/0
## 157  NW_019211454.1_17314319             0/0
## 158  NW_019211454.1_17320038            <NA>
## 159  NW_019211454.1_17322790             0/0
## 160  NW_019211454.1_17331569             0/0
## 161  NW_019211454.1_17344167            <NA>
## 162  NW_019211454.1_17350818            <NA>
## 163  NW_019211454.1_17352533            <NA>
## 164  NW_019211454.1_17374556             0/0
## 165  NW_019211454.1_17374579             0/0
## 166  NW_019211454.1_17375864             0/0
## 167  NW_019211454.1_17376941            <NA>
## 168  NW_019211454.1_17541597             0/0
## 169  NW_019211454.1_17551873             0/0
## 170  NW_019211454.1_17558887             0/0
## 171  NW_019211454.1_17576766             0/0
## 172  NW_019211454.1_17578196             0/0
## 173  NW_019211454.1_17583341             0/0
## 174  NW_019211454.1_17606030             0/0
## 175  NW_019211454.1_17610400             0/0
## 176  NW_019211454.1_17622107             0/0
## 177  NW_019211454.1_17772729            <NA>
## 178  NW_019211454.1_17803877            <NA>
## 179  NW_019211454.1_17837365             0/0
## 180  NW_019211454.1_17837385             0/0
## 181  NW_019211454.1_17837976             0/0
## 182  NW_019211454.1_17844163             0/0
## 183  NW_019211454.1_17954895            <NA>
## 184  NW_019211454.1_17955731             0/0
## 185  NW_019211454.1_17955924            <NA>
## 186  NW_019211454.1_17965364            <NA>
## 187  NW_019211454.1_17982176             0/0
## 188  NW_019211454.1_17986906             0/0
## 189  NW_019211454.1_17999036            <NA>
## 190  NW_019211454.1_18109983            <NA>
## 191  NW_019211454.1_18113968            <NA>
## 192  NW_019211454.1_18906639            <NA>
## 193  NW_019211454.1_21189838            <NA>
## 194  NW_019211454.1_21205822            <NA>
## 195  NW_019211454.1_23149732            <NA>
## 196  NW_019211454.1_23162575            <NA>
## 197  NW_019211454.1_23188311             0/0
## 198  NW_019211454.1_23192907            <NA>
## 199  NW_019211454.1_23199148            <NA>
## 200  NW_019211454.1_23212140            <NA>
## 201  NW_019211454.1_23218946             0/0
## 202  NW_019211454.1_23251760             0/0
## 203  NW_019211454.1_23254977             0/0
## 204  NW_019211454.1_23264224             0/0
## 205  NW_019211454.1_23268905            <NA>
## 206  NW_019211454.1_23272086            <NA>
## 207  NW_019211454.1_23276013             0/0
## 208  NW_019211454.1_23278009             0/0
## 209  NW_019211454.1_23298719            <NA>
## 210  NW_019211454.1_23299008            <NA>
## 211  NW_019211454.1_23300769            <NA>
## 212  NW_019211454.1_23300881            <NA>
## 213  NW_019211454.1_23301117             0/0
## 214  NW_019211454.1_23333824            <NA>
## 215  NW_019211454.1_23366149             0/0
## 216  NW_019211454.1_23382553             0/0
## 217  NW_019211454.1_23391385             0/0
## 218  NW_019211454.1_23410544            <NA>
## 219  NW_019211454.1_23425501             0/0
## 220  NW_019211454.1_23442161             0/0
## 221  NW_019211454.1_23577877            <NA>
## 222  NW_019211454.1_23595719            <NA>
## 223  NW_019211454.1_23624888            <NA>
## 224  NW_019211454.1_23634167            <NA>
## 225  NW_019211454.1_23637875             0/0
## 226  NW_019211454.1_23644002            <NA>
## 227  NW_019211454.1_24570146            <NA>
## 228  NW_019211454.1_24570473            <NA>
## 229  NW_019211454.1_24578003            <NA>
## 230  NW_019211454.1_24683480             0/0
## 231  NW_019211454.1_24727073             0/0
## 232  NW_019211454.1_24775786             0/0
## 233  NW_019211454.1_24819583             0/0
## 234  NW_019211454.1_27474697            <NA>
## 235  NW_019211454.1_27476484             0/0
## 236  NW_019211454.1_27491198             0/0
## 237  NW_019211454.1_27498234             0/0
## 238  NW_019211454.1_27498752             0/0
## 239  NW_019211454.1_27520650             0/0
## 240  NW_019211454.1_27531194            <NA>
## 241  NW_019211454.1_27540325            <NA>
## 242  NW_019211454.1_27580634             0/0
## 243  NW_019211454.1_27581718             0/0
## 244  NW_019211454.1_27582068             0/0
## 245  NW_019211454.1_27593680             0/0
## 246  NW_019211454.1_27594931            <NA>
## 247  NW_019211454.1_27596087             0/0
## 248  NW_019211454.1_27596126             0/0
## 249  NW_019211454.1_27612215             0/0
## 250  NW_019211454.1_27617579             0/0
## 251  NW_019211454.1_27699731            <NA>
## 252  NW_019211454.1_27727090            <NA>
## 253  NW_019211454.1_27735379            <NA>
## 254  NW_019211454.1_27738749             0/0
## 255  NW_019211454.1_27751954             0/0
## 256  NW_019211454.1_27753097             0/0
## 257  NW_019211454.1_27758233            <NA>
## 258  NW_019211454.1_27758304             0/0
## 259  NW_019211454.1_27758399             0/0
## 260  NW_019211454.1_27790477            <NA>
## 261  NW_019211454.1_27795450            <NA>
## 262  NW_019211454.1_27816715            <NA>
## 263  NW_019211454.1_27892109             0/0
## 264  NW_019211454.1_27892155             0/0
## 265  NW_019211454.1_27959399             0/0
## 266  NW_019211454.1_27961020            <NA>
## 267  NW_019211454.1_27996384             0/0
## 268  NW_019211454.1_28024213             0/0
## 269  NW_019211454.1_28039755            <NA>
## 270  NW_019211454.1_28065845             0/0
## 271  NW_019211454.1_28066125            <NA>
## 272  NW_019211454.1_29231990             0/0
## 273  NW_019211454.1_29232287             0/0
## 274  NW_019211454.1_29261601             0/0
## 275  NW_019211454.1_29919924             0/0
## 276  NW_019211454.1_30101274            <NA>
## 277  NW_019211454.1_30402966             0/0
## 278  NW_019211454.1_30419096            <NA>
## 279  NW_019211454.1_30426006             0/0
## 280  NW_019211454.1_30431348             0/0
## 281  NW_019211454.1_30435231            <NA>
## 282  NW_019211454.1_30447033            <NA>
## 283  NW_019211454.1_30454380            <NA>
## 284  NW_019211454.1_30454704            <NA>
## 285  NW_019211454.1_30457049             0/0
## 286  NW_019211454.1_30470808            <NA>
## 287  NW_019211454.1_30480822             0/0
## 288  NW_019211454.1_30498655            <NA>
## 289  NW_019211454.1_30502656            <NA>
## 290  NW_019211454.1_30505451             0/0
## 291  NW_019211454.1_30519580             0/0
## 292  NW_019211454.1_30520762            <NA>
## 293  NW_019211454.1_30527028             0/0
## 294  NW_019211454.1_30527215             0/0
## 295  NW_019211454.1_30528862            <NA>
## 296  NW_019211454.1_30756276            <NA>
## 297  NW_019211454.1_30781127             0/0
## 298  NW_019211454.1_30802390            <NA>
## 299  NW_019211454.1_30815022             0/0
## 300  NW_019211454.1_30850228             0/0
## 301  NW_019211454.1_31019741             0/0
## 302  NW_019211454.1_31117387             0/0
## 303  NW_019211454.1_32039254             0/0
## 304  NW_019211454.1_32054483             0/0
## 305  NW_019211454.1_32099381            <NA>
## 306  NW_019211454.1_32136801            <NA>
## 307  NW_019211454.1_32140314            <NA>
## 308  NW_019211454.1_32142383             0/0
## 309  NW_019211454.1_32143449             0/0
## 310  NW_019211454.1_32145009             0/0
## 311  NW_019211454.1_32156007            <NA>
## 312  NW_019211454.1_32166326            <NA>
## 313  NW_019211454.1_32171625             0/0
## 314  NW_019211454.1_32175316             0/0
## 315  NW_019211454.1_32186517            <NA>
## 316  NW_019211454.1_32209139            <NA>
## 317  NW_019211454.1_32278692            <NA>
## 318  NW_019211454.1_32293756             0/0
## 319  NW_019211454.1_32297314            <NA>
## 320  NW_019211454.1_32317836             0/0
## 321  NW_019211454.1_32964012             0/0
## 322  NW_019211454.1_33338832             0/0
## 323  NW_019211454.1_33352032             0/0
## 324  NW_019211454.1_33356322             0/0
## 325  NW_019211454.1_33357679            <NA>
## 326  NW_019211454.1_33372762            <NA>
## 327  NW_019211454.1_33421667            <NA>
## 328  NW_019211454.1_33448354            <NA>
## 329  NW_019211454.1_33452617             0/0
## 330  NW_019211454.1_33457401             0/0
## 331  NW_019211454.1_33470836             0/0
## 332  NW_019211454.1_33481705             0/0
## 333  NW_019211454.1_34661141             0/0
## 334  NW_019211454.1_34696147            <NA>
## 335  NW_019211454.1_34710662             0/0
## 336  NW_019211454.1_34735114             0/0
## 337  NW_019211454.1_34750224            <NA>
## 338  NW_019211454.1_34755204            <NA>
## 339  NW_019211454.1_34767703             0/0
## 340  NW_019211454.1_34769592             0/0
## 341  NW_019211454.1_34780333             0/0
## 342  NW_019211454.1_34780346             0/0
## 343  NW_019211454.1_34809110             0/0
## 344  NW_019211454.1_34843865            <NA>
## 345  NW_019211454.1_34871348            <NA>
## 346  NW_019211454.1_34947904             0/0
## 347  NW_019211454.1_35028182            <NA>
## 348  NW_019211454.1_35029455            <NA>
## 349  NW_019211454.1_35038046            <NA>
## 350  NW_019211454.1_35081471            <NA>
## 351  NW_019211454.1_35300566             0/0
## 352  NW_019211454.1_35531494            <NA>
## 353  NW_019211454.1_35537358             0/0
## 354  NW_019211454.1_35635742            <NA>
## 355  NW_019211454.1_35636953             0/0
## 356  NW_019211454.1_35657554             0/0
## 357  NW_019211454.1_36225107             0/0
## 358  NW_019211454.1_36231049             0/0
## 359  NW_019211454.1_36249867             0/0
## 360  NW_019211454.1_36303013            <NA>
## 361  NW_019211454.1_36303334            <NA>
## 362  NW_019211454.1_36305588            <NA>
## 363  NW_019211454.1_36305722             0/0
## 364  NW_019211454.1_36306986            <NA>
## 365  NW_019211454.1_36309014             0/0
## 366  NW_019211454.1_36312828            <NA>
## 367  NW_019211454.1_36322719             0/0
## 368  NW_019211454.1_36325864             0/0
## 369  NW_019211454.1_36327293             0/0
## 370  NW_019211454.1_36400301            <NA>
## 371  NW_019211454.1_36429883             0/0
## 372  NW_019211454.1_36431797             0/0
## 373  NW_019211454.1_36431811             0/0
## 374  NW_019211454.1_36434432             0/0
## 375  NW_019211454.1_36437879            <NA>
## 376  NW_019211454.1_36439773             0/0
## 377  NW_019211454.1_36456925             0/0
## 378  NW_019211454.1_36467442             0/0
## 379  NW_019211454.1_36476993            <NA>
## 380  NW_019211454.1_36485502            <NA>
## 381  NW_019211454.1_36487316            <NA>
## 382  NW_019211454.1_36488320             0/0
## 383  NW_019211454.1_36489299            <NA>
## 384  NW_019211454.1_36490354             0/0
## 385  NW_019211454.1_36492502            <NA>
## 386  NW_019211454.1_36517090             0/0
## 387  NW_019211454.1_36559306            <NA>
## 388  NW_019211454.1_36579567            <NA>
## 389  NW_019211454.1_36581439             0/0
## 390  NW_019211454.1_36592385            <NA>
## 391  NW_019211454.1_36594808            <NA>
## 392  NW_019211454.1_36625536             0/0
## 393  NW_019211454.1_36627365             0/0
## 394  NW_019211454.1_36651633             0/0
## 395  NW_019211454.1_36707318             0/0
## 396  NW_019211454.1_36757145            <NA>
## 397  NW_019211454.1_36761381             0/0
## 398  NW_019211454.1_36779314            <NA>
## 399  NW_019211454.1_36785002             0/0
## 400  NW_019211454.1_36786796             0/0
## 401  NW_019211454.1_36791692             0/1
## 402  NW_019211454.1_36799139             0/0
## 403  NW_019211454.1_36827402             0/0
## 404  NW_019211454.1_36827809             0/0
## 405  NW_019211454.1_36848731            <NA>
## 406  NW_019211454.1_36852527            <NA>
## 407  NW_019211454.1_36854992            <NA>
## 408  NW_019211454.1_36859160             0/0
## 409  NW_019211454.1_36860378            <NA>
## 410  NW_019211454.1_36866653            <NA>
## 411  NW_019211454.1_36870521             0/0
## 412  NW_019211454.1_36873027            <NA>
## 413  NW_019211454.1_36874395             0/0
## 414  NW_019211454.1_36881402            <NA>
## 415  NW_019211454.1_36885217             0/0
## 416  NW_019211454.1_36922811            <NA>
## 417  NW_019211454.1_36935960             0/0
## 418  NW_019211454.1_36941774            <NA>
## 419  NW_019211454.1_36961149            <NA>
## 420  NW_019211454.1_36963866            <NA>
## 421  NW_019211454.1_36970514             0/0
## 422  NW_019211454.1_36974827             0/0
## 423  NW_019211454.1_37003339            <NA>
## 424  NW_019211454.1_37003387            <NA>
## 425  NW_019211454.1_37010589            <NA>
## 426  NW_019211454.1_37011724            <NA>
## 427  NW_019211454.1_37024601             0/0
## 428  NW_019211454.1_37025330            <NA>
## 429  NW_019211454.1_37032799            <NA>
## 430  NW_019211454.1_37039254            <NA>
## 431  NW_019211454.1_37047454             0/0
## 432  NW_019211454.1_37049914             0/0
## 433  NW_019211454.1_37056997            <NA>
## 434  NW_019211454.1_37070881             0/0
## 435  NW_019211454.1_37100466            <NA>
## 436  NW_019211454.1_37128833             0/0
## 437  NW_019211454.1_37130569            <NA>
## 438  NW_019211454.1_37133936             0/0
## 439  NW_019211454.1_37134076             0/0
## 440  NW_019211454.1_37135404             0/0
## 441  NW_019211454.1_37135625             0/0
## 442  NW_019211454.1_37138864            <NA>
## 443  NW_019211454.1_37138887            <NA>
## 444  NW_019211454.1_37143549             0/0
## 445  NW_019211454.1_37753628            <NA>
## 446  NW_019211454.1_37910377            <NA>
## 447  NW_019211454.1_37938433            <NA>
## 448  NW_019211454.1_37940118            <NA>
## 449  NW_019211454.1_37940348             0/0
## 450  NW_019211454.1_37947399             0/0
## 451  NW_019211454.1_37948672            <NA>
## 452  NW_019211454.1_37987970            <NA>
## 453  NW_019211454.1_37987998            <NA>
## 454  NW_019211454.1_37988822            <NA>
## 455  NW_019211454.1_37989288            <NA>
## 456  NW_019211454.1_37992108             0/0
## 457  NW_019211454.1_37993898            <NA>
## 458  NW_019211454.1_37994047             0/0
## 459  NW_019211454.1_37994181            <NA>
## 460  NW_019211454.1_37998762            <NA>
## 461  NW_019211454.1_38000338            <NA>
## 462  NW_019211454.1_38001316            <NA>
## 463  NW_019211454.1_38001809             0/0
## 464  NW_019211454.1_38004005            <NA>
## 465  NW_019211454.1_38008590            <NA>
## 466  NW_019211454.1_38046727            <NA>
## 467  NW_019211454.1_38624087             0/0
## 468  NW_019211454.1_38639387             0/0
## 469  NW_019211454.1_38642178             0/0
## 470  NW_019211454.1_38646340             0/0
## 471  NW_019211454.1_38653750             0/0
## 472  NW_019211454.1_38693611            <NA>
## 473  NW_019211454.1_38694566            <NA>
## 474  NW_019211454.1_38710500             0/0
## 475  NW_019211454.1_38719649             0/0
## 476  NW_019211454.1_38719683             0/0
## 477  NW_019211454.1_38719882            <NA>
## 478  NW_019211454.1_38723980             0/0
## 479  NW_019211454.1_38728405            <NA>
## 480  NW_019211454.1_38731284            <NA>
## 481  NW_019211454.1_38732279             0/0
## 482  NW_019211454.1_38734590             0/0
## 483  NW_019211454.1_38760743             0/0
## 484  NW_019211454.1_38760785             0/0
## 485  NW_019211454.1_38760800             0/0
## 486  NW_019211454.1_38772689             0/0
## 487  NW_019211454.1_39463236             0/0
## 488  NW_019211454.1_39511091            <NA>
## 489  NW_019211454.1_39719447             0/0
## 490  NW_019211454.1_39762977             0/0
## 491  NW_019211454.1_39788462             0/0
## 492  NW_019211454.1_39798196            <NA>
## 493  NW_019211454.1_39808695            <NA>
## 494  NW_019211454.1_39810900            <NA>
## 495  NW_019211454.1_40034273            <NA>
## 496  NW_019211454.1_40070100             0/0
## 497  NW_019211454.1_40070316             0/0
## 498  NW_019211454.1_40089652             0/0
## 499  NW_019211454.1_40090061             0/0
## 500  NW_019211454.1_40114751            <NA>
## 501  NW_019211454.1_40126972             0/0
## 502  NW_019211454.1_40129662             0/0
## 503  NW_019211454.1_40132976             0/0
## 504  NW_019211454.1_40144234            <NA>
## 505  NW_019211454.1_40147921            <NA>
## 506  NW_019211454.1_40147937            <NA>
## 507  NW_019211454.1_40157895             0/0
## 508  NW_019211454.1_40158668             0/0
## 509  NW_019211454.1_40158739             0/0
## 510  NW_019211454.1_40162517             0/0
## 511  NW_019211454.1_40167978             0/0
## 512  NW_019211454.1_40169514             0/0
## 513  NW_019211454.1_40170666             0/0
## 514  NW_019211454.1_40187999             0/0
## 515  NW_019211454.1_40206130             0/0
## 516  NW_019211454.1_40208386             0/0
## 517  NW_019211454.1_40211570             0/0
## 518  NW_019211454.1_40227396            <NA>
## 519  NW_019211454.1_40241170             0/0
## 520  NW_019211454.1_40243531            <NA>
## 521  NW_019211454.1_40250427             0/0
## 522  NW_019211454.1_40261072             0/0
## 523  NW_019211454.1_40266318             0/0
## 524  NW_019211454.1_40276818            <NA>
## 525  NW_019211454.1_40280330             0/0
## 526  NW_019211454.1_40626334             0/0
## 527  NW_019211454.1_40630404            <NA>
## 528  NW_019211454.1_40631887             0/0
## 529  NW_019211454.1_40638403             0/0
## 530  NW_019211454.1_40666283             0/0
## 531  NW_019211454.1_40680609             0/0
## 532  NW_019211454.1_40683288            <NA>
## 533  NW_019211454.1_40693798            <NA>
## 534  NW_019211454.1_40698206             0/0
## 535  NW_019211454.1_40717536             0/0
## 536  NW_019211454.1_40721776             0/0
## 537  NW_019211454.1_40722145            <NA>
## 538  NW_019211454.1_40725411             0/0
## 539  NW_019211454.1_40725787             0/0
## 540  NW_019211454.1_41101184            <NA>
## 541  NW_019211454.1_41155118             0/0
## 542  NW_019211454.1_41157815             0/0
## 543  NW_019211454.1_41158461            <NA>
## 544  NW_019211454.1_41176200             0/0
## 545  NW_019211454.1_41178977             0/0
## 546  NW_019211454.1_41187641            <NA>
## 547  NW_019211454.1_41189701             0/0
## 548  NW_019211454.1_41195384             0/0
## 549  NW_019211454.1_41198484            <NA>
## 550  NW_019211454.1_41200614             0/0
## 551  NW_019211454.1_41200809             0/1
## 552  NW_019211454.1_41201413             0/0
## 553  NW_019211454.1_41213762            <NA>
## 554  NW_019211454.1_41220506             0/0
## 555  NW_019211454.1_41230368             0/0
## 556  NW_019211454.1_41240538             0/0
## 557  NW_019211454.1_41244157             0/0
## 558  NW_019211454.1_41250910            <NA>
## 559  NW_019211454.1_41252294             0/0
## 560  NW_019211454.1_41271694             0/0
## 561  NW_019211454.1_41275567            <NA>
## 562  NW_019211454.1_41276956            <NA>
## 563  NW_019211454.1_41276969             0/0
## 564  NW_019211454.1_41347283            <NA>
## 565  NW_019211454.1_41379837            <NA>
## 566  NW_019211454.1_41389149             0/0
## 567  NW_019211454.1_44066323             0/0
## 568  NW_019211454.1_44075125             0/0
## 569  NW_019211454.1_44103469             0/0
## 570  NW_019211454.1_44105552             0/0
## 571  NW_019211454.1_44112775            <NA>
## 572  NW_019211454.1_44140420             0/0
## 573  NW_019211454.1_46214928            <NA>
## 574  NW_019211454.1_48486198             0/0
## 575  NW_019211454.1_48491589             0/0
## 576  NW_019211454.1_48525461            <NA>
## 577  NW_019211454.1_49033072             0/0
## 578  NW_019211454.1_49033419             0/0
## 579  NW_019211454.1_49035336            <NA>
## 580  NW_019211454.1_49043289             0/0
## 581  NW_019211454.1_49057108            <NA>
## 582  NW_019211454.1_49378114             0/0
## 583  NW_019211454.1_50191124             0/0
## 584  NW_019211454.1_52386002            <NA>
## 585  NW_019211454.1_53431088            <NA>
## 586  NW_019211454.1_56129168            <NA>
## 587  NW_019211454.1_56754553            <NA>
## 588  NW_019211454.1_56790992             0/0
## 589  NW_019211454.1_58461103             0/0
## 590  NW_019211454.1_60040889             0/0
## 591  NW_019211454.1_60045747             0/0
## 592  NW_019211454.1_60052833             0/0
## 593  NW_019211454.1_60070105             0/0
## 594  NW_019211454.1_60084987            <NA>
## 595  NW_019211454.1_60086106             0/0
## 596  NW_019211454.1_60090982            <NA>
## 597  NW_019211454.1_60091330            <NA>
## 598  NW_019211454.1_60094151             0/0
## 599  NW_019211454.1_60094842             0/0
## 600  NW_019211454.1_60096936            <NA>
## 601  NW_019211454.1_60116457             0/0
## 602  NW_019211454.1_60116528            <NA>
## 603  NW_019211454.1_60119713             0/0
## 604  NW_019211454.1_60120460            <NA>
## 605  NW_019211454.1_60126262            <NA>
## 606  NW_019211454.1_60143176            <NA>
## 607  NW_019211454.1_60150897            <NA>
## 608  NW_019211454.1_60178344             0/0
## 609  NW_019211454.1_60199503             0/0
## 610  NW_019211454.1_60214441            <NA>
## 611  NW_019211454.1_60214456            <NA>
## 612  NW_019211454.1_60227203             0/0
## 613  NW_019211454.1_60233060             0/0
## 614  NW_019211454.1_60294325             0/0
## 615  NW_019211454.1_60295124            <NA>
## 616  NW_019211454.1_60300810            <NA>
## 617  NW_019211454.1_60300845            <NA>
## 618  NW_019211454.1_60302269             0/0
## 619  NW_019211454.1_60309499             0/0
## 620  NW_019211454.1_60314137            <NA>
## 621  NW_019211454.1_60870511            <NA>
## 622  NW_019211454.1_60871712            <NA>
## 623  NW_019211454.1_60886131             0/0
## 624  NW_019211454.1_60892756            <NA>
## 625  NW_019211454.1_60898037            <NA>
## 626  NW_019211454.1_60903426             0/0
## 627  NW_019211454.1_60910004             0/0
## 628  NW_019211454.1_60918838             0/0
## 629  NW_019211454.1_60936792             0/0
## 630  NW_019211454.1_60999522             0/0
## 631  NW_019211454.1_61005443            <NA>
## 632  NW_019211454.1_61010017             0/0
## 633  NW_019211454.1_61036882             0/0
## 634  NW_019211454.1_61059853             0/0
## 635  NW_019211454.1_61060838             0/0
## 636  NW_019211454.1_61060872             0/0
## 637  NW_019211454.1_61062560             0/0
## 638  NW_019211454.1_61064207             0/0
## 639  NW_019211454.1_61066588             0/0
## 640  NW_019211454.1_61071808            <NA>
## 641  NW_019211454.1_61100520            <NA>
## 642  NW_019211454.1_62221894             0/0
## 643  NW_019211454.1_62223198             0/0
## 644  NW_019211454.1_62228815            <NA>
## 645  NW_019211454.1_62253745            <NA>
## 646  NW_019211454.1_62297656             0/0
## 647  NW_019211454.1_62362013             0/0
## 648  NW_019211454.1_62379659            <NA>
## 649  NW_019211454.1_62384719            <NA>
## 650  NW_019211454.1_62742862            <NA>
## 651  NW_019211454.1_62744523             0/0
## 652  NW_019211454.1_62755203             0/0
## 653  NW_019211454.1_62760791            <NA>
## 654  NW_019211454.1_62766087             0/0
## 655  NW_019211454.1_62845696            <NA>
## 656  NW_019211454.1_66357056            <NA>
## 657  NW_019211454.1_67642464            <NA>
## 658  NW_019211454.1_68070552            <NA>
## 659  NW_019211454.1_68081673            <NA>
## 660  NW_019211454.1_68208481             0/0
## 661  NW_019211454.1_68210796             0/0
## 662  NW_019211454.1_68249627            <NA>
## 663  NW_019211454.1_68249776             0/0
## 664  NW_019211454.1_68250522            <NA>
## 665  NW_019211454.1_68252004             0/0
## 666  NW_019211454.1_68252460            <NA>
## 667  NW_019211454.1_68268811            <NA>
## 668  NW_019211454.1_68303585            <NA>
## 669  NW_019211454.1_68310183             0/0
## 670  NW_019211454.1_68334629            <NA>
## 671  NW_019211454.1_68385302            <NA>
## 672  NW_019211454.1_68390801            <NA>
## 673  NW_019211454.1_68392311            <NA>
## 674  NW_019211454.1_68410266             0/0
## 675  NW_019211454.1_68426227             0/0
## 676  NW_019211454.1_68436675             0/0
## 677  NW_019211454.1_68524498             0/0
## 678  NW_019211454.1_68539594             0/0
## 679  NW_019211454.1_68539621             0/0
## 680  NW_019211454.1_68565197             0/0
## 681  NW_019211454.1_69692222             0/0
## 682  NW_019211454.1_69694549            <NA>
## 683  NW_019211454.1_69694715             0/0
## 684  NW_019211454.1_69702475             0/0
## 685  NW_019211454.1_69765615             0/0
## 686  NW_019211454.1_69772933             0/0
## 687  NW_019211454.1_69777852             0/0
## 688  NW_019211454.1_69795954            <NA>
## 689  NW_019211454.1_70186883             0/0
## 690  NW_019211454.1_72789477             0/0
## 691  NW_019211454.1_72789486             0/0
## 692  NW_019211454.1_72789502             0/0
## 693  NW_019211454.1_72789511             0/0
## 694  NW_019211454.1_74126558            <NA>
## 695  NW_019211454.1_74128294            <NA>
## 696  NW_019211454.1_74290774             0/0
## 697  NW_019211454.1_74292728             0/0
## 698  NW_019211454.1_74341144            <NA>
## 699  NW_019211454.1_74352633             0/0
## 700  NW_019211454.1_74375274            <NA>
## 701  NW_019211454.1_74383872            <NA>
## 702  NW_019211454.1_74386880             0/0
## 703  NW_019211454.1_74397070            <NA>
## 704  NW_019211454.1_74668449             0/0
## 705  NW_019211454.1_75034547            <NA>
## 706  NW_019211454.1_75046501            <NA>
## 707  NW_019211454.1_75059444            <NA>
## 708  NW_019211454.1_75090669            <NA>
## 709  NW_019211454.1_75101747             0/0
## 710  NW_019211454.1_75101827             0/0
## 711  NW_019211454.1_75101887             0/0
## 712  NW_019211454.1_75112840             0/0
## 713  NW_019211454.1_75133554            <NA>
## 714  NW_019211454.1_75136801             0/0
## 715  NW_019211454.1_75147647             0/0
## 716  NW_019211454.1_75149807            <NA>
## 717  NW_019211454.1_75158568             0/0
## 718  NW_019211454.1_75163728             0/0
## 719  NW_019211454.1_75194641            <NA>
## 720  NW_019211454.1_75195645            <NA>
## 721  NW_019211454.1_76223546            <NA>
## 722  NW_019211454.1_76223565             0/0
## 723  NW_019211454.1_76234397            <NA>
## 724  NW_019211454.1_76591412            <NA>
## 725    NW_019211455.1_297151            <NA>
## 726   NW_019211455.1_1723309             0/0
## 727   NW_019211455.1_1735588             0/0
## 728   NW_019211455.1_1755178             0/0
## 729   NW_019211455.1_1757742            <NA>
## 730   NW_019211455.1_1759287            <NA>
## 731   NW_019211455.1_1786467            <NA>
## 732   NW_019211455.1_1790136            <NA>
## 733   NW_019211455.1_1796024             0/0
## 734   NW_019211455.1_1799665             0/0
## 735   NW_019211455.1_1817006            <NA>
## 736   NW_019211455.1_1817248            <NA>
## 737   NW_019211455.1_1818976            <NA>
## 738   NW_019211455.1_1827439             0/0
## 739   NW_019211455.1_1844582            <NA>
## 740   NW_019211455.1_1998846             0/0
## 741   NW_019211455.1_2002301             0/0
## 742   NW_019211455.1_2003141            <NA>
## 743   NW_019211455.1_2018035             0/0
## 744   NW_019211455.1_2024149             0/0
## 745   NW_019211455.1_2035966            <NA>
## 746   NW_019211455.1_2037845            <NA>
## 747   NW_019211455.1_2049751            <NA>
## 748   NW_019211455.1_2055201             0/0
## 749   NW_019211455.1_2086517             0/0
## 750   NW_019211455.1_2090967             0/0
## 751   NW_019211455.1_2105994             0/0
## 752   NW_019211455.1_2156648             0/0
## 753   NW_019211455.1_2269993             0/0
## 754   NW_019211455.1_2270186             0/0
## 755   NW_019211455.1_2273866             0/0
## 756   NW_019211455.1_6296845             0/0
## 757   NW_019211455.1_6297935            <NA>
## 758   NW_019211455.1_6310894            <NA>
## 759   NW_019211455.1_6311013            <NA>
## 760   NW_019211455.1_6313858            <NA>
## 761   NW_019211455.1_6328519            <NA>
## 762   NW_019211455.1_6341511            <NA>
## 763   NW_019211455.1_6342371            <NA>
## 764   NW_019211455.1_6388395            <NA>
## 765   NW_019211455.1_6862184            <NA>
## 766   NW_019211455.1_6892496             0/0
## 767   NW_019211455.1_6901533             0/0
## 768   NW_019211455.1_7375407             0/0
## 769   NW_019211455.1_7621250            <NA>
## 770   NW_019211455.1_7760562             0/0
## 771   NW_019211455.1_9623242             0/0
## 772   NW_019211455.1_9624267             0/0
## 773   NW_019211455.1_9986181            <NA>
## 774  NW_019211455.1_10308072             0/0
## 775  NW_019211455.1_10351749             0/0
## 776  NW_019211455.1_10712687             0/0
## 777  NW_019211455.1_10724719            <NA>
## 778  NW_019211455.1_10739248             0/0
## 779  NW_019211455.1_10787412             0/0
## 780  NW_019211455.1_10797398             0/0
## 781  NW_019211455.1_10802538            <NA>
## 782  NW_019211455.1_10815261            <NA>
## 783  NW_019211455.1_10834864            <NA>
## 784  NW_019211455.1_10837041             0/0
## 785  NW_019211455.1_10842021             0/0
## 786  NW_019211455.1_12991126            <NA>
## 787  NW_019211455.1_17149864             0/0
## 788  NW_019211455.1_17185737             0/0
## 789  NW_019211455.1_17252612             0/0
## 790  NW_019211455.1_17563362             0/0
## 791  NW_019211455.1_17585370             0/0
## 792  NW_019211455.1_17595372             0/0
## 793  NW_019211455.1_17653912             0/0
## 794  NW_019211455.1_17673820             0/0
## 795  NW_019211455.1_17678811             0/0
## 796  NW_019211455.1_17740334            <NA>
## 797  NW_019211455.1_17797673            <NA>
## 798  NW_019211455.1_17817512            <NA>
## 799  NW_019211455.1_17879102             0/0
## 800  NW_019211455.1_17882985             0/0
## 801  NW_019211455.1_17896957             0/0
## 802  NW_019211455.1_17906927            <NA>
## 803  NW_019211455.1_17982258             0/0
## 804  NW_019211455.1_18353510             0/0
## 805  NW_019211455.1_22222784            <NA>
## 806  NW_019211455.1_23246093             0/0
## 807  NW_019211455.1_23495473             0/0
## 808  NW_019211455.1_24612052             0/0
## 809  NW_019211455.1_25486809             0/0
## 810  NW_019211455.1_25625408            <NA>
## 811  NW_019211455.1_25882252             0/0
## 812  NW_019211455.1_25905800            <NA>
## 813  NW_019211455.1_25972186            <NA>
## 814  NW_019211455.1_25980384             0/0
## 815  NW_019211455.1_25996836            <NA>
## 816  NW_019211455.1_26000220             0/0
## 817  NW_019211455.1_26031305            <NA>
## 818  NW_019211455.1_26130475             0/0
## 819  NW_019211455.1_26465563            <NA>
## 820  NW_019211455.1_27706793            <NA>
## 821  NW_019211455.1_27775552            <NA>
## 822  NW_019211455.1_27797560            <NA>
## 823  NW_019211455.1_27819122             0/0
## 824  NW_019211455.1_27826549             0/0
## 825  NW_019211455.1_27826921            <NA>
## 826  NW_019211455.1_27835138            <NA>
## 827  NW_019211455.1_27835266            <NA>
## 828  NW_019211455.1_27835677             0/0
## 829  NW_019211455.1_27838161             0/0
## 830  NW_019211455.1_27840881             0/0
## 831  NW_019211455.1_27957912             0/0
## 832  NW_019211455.1_28957975             0/0
## 833  NW_019211455.1_28986332            <NA>
## 834  NW_019211455.1_28986700             0/0
## 835  NW_019211455.1_29111946             0/0
## 836  NW_019211455.1_29120354            <NA>
## 837  NW_019211455.1_29127976            <NA>
## 838  NW_019211455.1_29135396            <NA>
## 839  NW_019211455.1_29140383             0/0
## 840  NW_019211455.1_29142218            <NA>
## 841  NW_019211455.1_29157976             0/0
## 842  NW_019211455.1_29158671             0/0
## 843  NW_019211455.1_29177263            <NA>
## 844  NW_019211455.1_29180050            <NA>
## 845  NW_019211455.1_29185410             0/0
## 846  NW_019211455.1_29186191             0/0
## 847  NW_019211455.1_29200301             0/0
## 848  NW_019211455.1_29202343             0/0
## 849  NW_019211455.1_29260943            <NA>
## 850  NW_019211455.1_29261900             0/0
## 851  NW_019211455.1_29353407            <NA>
## 852  NW_019211455.1_29362379            <NA>
## 853  NW_019211455.1_29387221             0/0
## 854  NW_019211455.1_30201474             0/0
## 855  NW_019211455.1_30227266             0/0
## 856  NW_019211455.1_30231935             0/0
## 857  NW_019211455.1_30738434             0/0
## 858  NW_019211455.1_31207747            <NA>
## 859  NW_019211455.1_32296843             0/0
## 860  NW_019211455.1_32337677             0/0
## 861  NW_019211455.1_32339825             0/0
## 862  NW_019211455.1_32339846             0/0
## 863  NW_019211455.1_32343219            <NA>
## 864  NW_019211455.1_32346622            <NA>
## 865  NW_019211455.1_32356611            <NA>
## 866  NW_019211455.1_32361819            <NA>
## 867  NW_019211455.1_32363967            <NA>
## 868  NW_019211455.1_32369078            <NA>
## 869  NW_019211455.1_32376932            <NA>
## 870  NW_019211455.1_32388092            <NA>
## 871  NW_019211455.1_32390675            <NA>
## 872  NW_019211455.1_32418640             0/0
## 873  NW_019211455.1_32419186            <NA>
## 874  NW_019211455.1_32426852            <NA>
## 875  NW_019211455.1_32430044             0/0
## 876  NW_019211455.1_32434705             0/0
## 877  NW_019211455.1_32441005             0/0
## 878  NW_019211455.1_32445017             0/0
## 879  NW_019211455.1_32468374             0/0
## 880  NW_019211455.1_32469897             0/0
## 881  NW_019211455.1_32470012             0/0
## 882  NW_019211455.1_32474203            <NA>
## 883  NW_019211455.1_32493142            <NA>
## 884  NW_019211455.1_32496172            <NA>
## 885  NW_019211455.1_32497410             0/0
## 886  NW_019211455.1_32507379             0/0
## 887  NW_019211455.1_32512867             0/0
## 888  NW_019211455.1_32515507             0/0
## 889  NW_019211455.1_32522896             0/0
## 890  NW_019211455.1_32523059            <NA>
## 891  NW_019211455.1_32553392             0/0
## 892  NW_019211455.1_32881736             0/0
## 893  NW_019211455.1_32894397            <NA>
## 894  NW_019211455.1_32894669             0/0
## 895  NW_019211455.1_32938799             0/0
## 896  NW_019211455.1_32952007             0/0
## 897  NW_019211455.1_32963757             0/0
## 898  NW_019211455.1_32970140            <NA>
## 899  NW_019211455.1_32971481            <NA>
## 900  NW_019211455.1_32971605             0/0
## 901  NW_019211455.1_32981224            <NA>
## 902  NW_019211455.1_32984952             0/0
## 903  NW_019211455.1_32985410            <NA>
## 904  NW_019211455.1_32994239            <NA>
## 905  NW_019211455.1_33027539             0/0
## 906  NW_019211455.1_33054097             0/0
## 907  NW_019211455.1_33057384             0/0
## 908  NW_019211455.1_33069852            <NA>
## 909  NW_019211455.1_33077797             0/0
## 910  NW_019211455.1_33081788            <NA>
## 911  NW_019211455.1_33081937            <NA>
## 912  NW_019211455.1_33084526            <NA>
## 913  NW_019211455.1_33085106             0/0
## 914  NW_019211455.1_33087714            <NA>
## 915  NW_019211455.1_33092936            <NA>
## 916  NW_019211455.1_33102611            <NA>
## 917  NW_019211455.1_33114831             0/0
## 918  NW_019211455.1_33123530             0/0
## 919  NW_019211455.1_33126806             0/0
## 920  NW_019211455.1_33130026            <NA>
## 921  NW_019211455.1_33200929             0/0
## 922  NW_019211455.1_33226353             0/0
## 923  NW_019211455.1_33262612             0/0
## 924  NW_019211455.1_33292229             0/0
## 925  NW_019211455.1_33293808             0/0
## 926  NW_019211455.1_33299254            <NA>
## 927  NW_019211455.1_33300015             0/0
## 928  NW_019211455.1_33302561             0/0
## 929  NW_019211455.1_33313586             0/0
## 930  NW_019211455.1_33319912            <NA>
## 931  NW_019211455.1_33320213             0/0
## 932  NW_019211455.1_33321615            <NA>
## 933  NW_019211455.1_33350023             0/0
## 934  NW_019211455.1_33352594             0/0
## 935  NW_019211455.1_33359991            <NA>
## 936  NW_019211455.1_33362669            <NA>
## 937  NW_019211455.1_33374791             0/0
## 938  NW_019211455.1_33377643            <NA>
## 939  NW_019211455.1_33403086            <NA>
## 940  NW_019211455.1_33423204             0/0
## 941  NW_019211455.1_33438484            <NA>
## 942  NW_019211455.1_33439183             0/0
## 943  NW_019211455.1_33466919             0/0
## 944  NW_019211455.1_33476515             0/0
## 945  NW_019211455.1_33486164            <NA>
## 946  NW_019211455.1_33532935             0/0
## 947  NW_019211455.1_33556976             0/0
## 948  NW_019211455.1_33603645             0/0
## 949  NW_019211455.1_33630706             0/0
## 950  NW_019211455.1_33635799             0/0
## 951  NW_019211455.1_33745725            <NA>
## 952  NW_019211455.1_33751908             0/0
## 953  NW_019211455.1_33752085             0/0
## 954  NW_019211455.1_33843457             0/0
## 955  NW_019211455.1_33885647            <NA>
## 956  NW_019211455.1_33954535             0/0
## 957  NW_019211455.1_33967120             0/0
## 958  NW_019211455.1_33970776            <NA>
## 959  NW_019211455.1_43992022             0/0
## 960  NW_019211455.1_44010288             0/0
## 961  NW_019211455.1_44018236             0/0
## 962  NW_019211455.1_44080482            <NA>
## 963  NW_019211455.1_44094732            <NA>
## 964  NW_019211455.1_44115058             0/0
## 965  NW_019211455.1_44126474            <NA>
## 966  NW_019211455.1_44150586            <NA>
## 967  NW_019211455.1_44150866            <NA>
## 968  NW_019211455.1_44168728             0/0
## 969  NW_019211455.1_44169811             0/0
## 970  NW_019211455.1_44273021             0/0
## 971  NW_019211455.1_44345826             0/0
## 972  NW_019211455.1_44379285             0/0
## 973  NW_019211455.1_44425219            <NA>
## 974  NW_019211455.1_45942338             0/0
## 975  NW_019211455.1_45947753             0/0
## 976  NW_019211455.1_45952601            <NA>
## 977  NW_019211455.1_45999743             0/0
## 978  NW_019211455.1_46023005            <NA>
## 979  NW_019211455.1_46056725             0/0
## 980  NW_019211455.1_46143785             0/0
## 981  NW_019211455.1_46173048             0/0
## 982  NW_019211455.1_46212761             0/0
## 983  NW_019211455.1_46230879            <NA>
## 984  NW_019211455.1_46232773            <NA>
## 985  NW_019211455.1_46280279             0/0
## 986  NW_019211455.1_46352577             0/0
## 987  NW_019211455.1_46358307             0/0
## 988  NW_019211455.1_46398917             0/0
## 989  NW_019211455.1_46402800            <NA>
## 990  NW_019211455.1_46403999            <NA>
## 991  NW_019211455.1_46422811             0/1
## 992  NW_019211455.1_46422983             0/0
## 993  NW_019211455.1_46438010             0/0
## 994  NW_019211455.1_46463185            <NA>
## 995  NW_019211455.1_46482603             0/0
## 996  NW_019211455.1_46490910            <NA>
## 997  NW_019211455.1_46540183            <NA>
## 998  NW_019211455.1_46569030            <NA>
## 999  NW_019211455.1_46628280            <NA>
## 1000 NW_019211455.1_46635397             0/0
## 1001 NW_019211455.1_46645131             0/0
## 1002 NW_019211455.1_46648663             0/0
## 1003 NW_019211455.1_46657450             0/0
## 1004 NW_019211455.1_46715628             0/0
## 1005 NW_019211455.1_46758223             0/0
## 1006 NW_019211455.1_46758853            <NA>
## 1007 NW_019211455.1_46835426             0/0
## 1008 NW_019211455.1_46840769            <NA>
## 1009 NW_019211455.1_46856747             0/0
## 1010 NW_019211455.1_46856806             0/0
## 1011 NW_019211455.1_46884511             0/0
## 1012 NW_019211455.1_46898007             0/0
## 1013 NW_019211455.1_46914680            <NA>
## 1014 NW_019211455.1_46918514            <NA>
## 1015 NW_019211455.1_46953584             0/0
## 1016 NW_019211455.1_46984811             0/0
## 1017 NW_019211455.1_47022470             0/0
## 1018 NW_019211455.1_47042771             0/0
## 1019 NW_019211455.1_47077891            <NA>
## 1020 NW_019211455.1_47144157            <NA>
## 1021 NW_019211455.1_47179282            <NA>
## 1022 NW_019211455.1_47179293            <NA>
## 1023 NW_019211455.1_47179374            <NA>
## 1024 NW_019211455.1_47193481             0/0
## 1025 NW_019211455.1_47202010             0/0
## 1026 NW_019211455.1_47202462             0/0
## 1027 NW_019211455.1_47212866            <NA>
## 1028 NW_019211455.1_47237571             0/0
## 1029 NW_019211455.1_47248373            <NA>
## 1030 NW_019211455.1_47401099             0/0
## 1031 NW_019211455.1_47425962            <NA>
## 1032 NW_019211455.1_47452310             0/0
## 1033 NW_019211455.1_47701640            <NA>
## 1034 NW_019211455.1_48198008             0/0
## 1035 NW_019211455.1_48208542             0/0
## 1036 NW_019211455.1_48245617            <NA>
## 1037 NW_019211455.1_48258947            <NA>
## 1038 NW_019211455.1_48258992            <NA>
## 1039 NW_019211455.1_51129038            <NA>
## 1040 NW_019211455.1_51307592             0/0
## 1041 NW_019211455.1_52825989             0/0
## 1042 NW_019211455.1_52846975            <NA>
## 1043 NW_019211455.1_52879066             0/0
## 1044 NW_019211455.1_52882124            <NA>
## 1045 NW_019211455.1_52886915             0/0
## 1046 NW_019211455.1_52890367             0/0
## 1047 NW_019211455.1_52893316            <NA>
## 1048 NW_019211455.1_52911389             0/0
## 1049 NW_019211455.1_52937516             0/0
## 1050 NW_019211455.1_52943625            <NA>
## 1051 NW_019211455.1_52943738             0/0
## 1052 NW_019211455.1_52947377            <NA>
## 1053 NW_019211455.1_52995750             0/0
## 1054 NW_019211455.1_53540635             0/0
## 1055 NW_019211455.1_54380106            <NA>
## 1056 NW_019211455.1_54381574            <NA>
## 1057 NW_019211455.1_54386280             0/0
## 1058 NW_019211455.1_54389650            <NA>
## 1059 NW_019211455.1_54395668            <NA>
## 1060 NW_019211455.1_54402567             0/0
## 1061 NW_019211455.1_54411340            <NA>
## 1062 NW_019211455.1_54421016            <NA>
## 1063 NW_019211455.1_54491786             0/0
## 1064 NW_019211455.1_54491868            <NA>
## 1065 NW_019211455.1_54493357            <NA>
## 1066 NW_019211455.1_54498426            <NA>
## 1067 NW_019211455.1_54527577             0/0
## 1068 NW_019211455.1_54553702             0/0
## 1069 NW_019211455.1_54576375            <NA>
## 1070 NW_019211455.1_54607362            <NA>
## 1071 NW_019211455.1_54624652             0/0
## 1072 NW_019211455.1_54644131             0/0
## 1073 NW_019211455.1_54651895             0/0
## 1074 NW_019211455.1_54657681             0/0
## 1075 NW_019211455.1_54669062            <NA>
## 1076 NW_019211455.1_54669579             0/0
## 1077 NW_019211455.1_54687546             0/0
## 1078 NW_019211455.1_54717899             0/0
## 1079 NW_019211455.1_54724953             0/0
## 1080 NW_019211455.1_54732453             0/0
## 1081 NW_019211455.1_54737739             0/0
## 1082 NW_019211455.1_54737924             0/0
## 1083 NW_019211455.1_54758026            <NA>
## 1084 NW_019211455.1_54759583            <NA>
## 1085 NW_019211455.1_54762916            <NA>
## 1086 NW_019211455.1_54782102             0/0
## 1087 NW_019211455.1_54783615             0/0
## 1088 NW_019211455.1_54863497             0/0
## 1089 NW_019211455.1_54879180             0/0
## 1090 NW_019211455.1_54883036             0/0
## 1091 NW_019211455.1_54965397             0/0
## 1092 NW_019211455.1_54979264             0/0
## 1093 NW_019211455.1_54983705             0/0
## 1094 NW_019211455.1_55648597             0/0
## 1095 NW_019211455.1_55689498             0/0
## 1096 NW_019211455.1_55705842             0/0
## 1097 NW_019211455.1_55776197            <NA>
## 1098 NW_019211455.1_55802770            <NA>
## 1099 NW_019211455.1_55855535             0/0
## 1100 NW_019211455.1_55857271             0/0
## 1101 NW_019211455.1_55862495            <NA>
## 1102 NW_019211455.1_55873586             0/0
## 1103 NW_019211455.1_55895887            <NA>
## 1104 NW_019211455.1_55899517            <NA>
## 1105 NW_019211455.1_55904098            <NA>
## 1106 NW_019211455.1_55941737            <NA>
## 1107 NW_019211455.1_55958958             0/0
## 1108 NW_019211455.1_55998586             0/0
## 1109 NW_019211455.1_56308420            <NA>
## 1110 NW_019211455.1_56962613            <NA>
## 1111   NW_019211456.1_852213            <NA>
## 1112   NW_019211456.1_862934             0/0
## 1113  NW_019211456.1_1514091             0/0
## 1114  NW_019211456.1_1557845             0/0
## 1115  NW_019211456.1_1557913             0/0
## 1116  NW_019211456.1_1575938            <NA>
## 1117  NW_019211456.1_1576213             0/0
## 1118  NW_019211456.1_1581618            <NA>
## 1119  NW_019211456.1_1583683             0/0
## 1120  NW_019211456.1_1634695             0/0
## 1121  NW_019211456.1_1657692            <NA>
## 1122  NW_019211456.1_1662405             0/0
## 1123  NW_019211456.1_1671159            <NA>
## 1124  NW_019211456.1_1682880             0/0
## 1125  NW_019211456.1_3857748             0/0
## 1126  NW_019211456.1_3886580             0/0
## 1127  NW_019211456.1_3906070            <NA>
## 1128  NW_019211456.1_3932173             0/0
## 1129  NW_019211456.1_3938341             0/0
## 1130  NW_019211456.1_3941809             0/0
## 1131  NW_019211456.1_3960878             0/0
## 1132  NW_019211456.1_3983269             0/0
## 1133  NW_019211456.1_3996604            <NA>
## 1134  NW_019211456.1_4007274             0/0
## 1135  NW_019211456.1_4010857             0/0
## 1136  NW_019211456.1_4035995             0/0
## 1137  NW_019211456.1_4049358            <NA>
## 1138  NW_019211456.1_4061699             0/0
## 1139  NW_019211456.1_4067032             0/0
## 1140  NW_019211456.1_4092318             0/0
## 1141  NW_019211456.1_4097747             0/0
## 1142  NW_019211456.1_4098080             0/0
## 1143  NW_019211456.1_4123215             0/0
## 1144  NW_019211456.1_4138181            <NA>
## 1145  NW_019211456.1_4303393             0/0
## 1146  NW_019211456.1_5959685            <NA>
## 1147  NW_019211456.1_8208869             0/0
## 1148  NW_019211456.1_9135241            <NA>
## 1149  NW_019211456.1_9568584            <NA>
## 1150  NW_019211456.1_9614543            <NA>
## 1151  NW_019211456.1_9615304             0/0
## 1152 NW_019211456.1_10707372             0/1
## 1153 NW_019211456.1_11003117             0/0
## 1154 NW_019211456.1_11028134             0/0
## 1155 NW_019211456.1_12647375            <NA>
## 1156 NW_019211456.1_12878824             0/0
## 1157 NW_019211456.1_12880280             0/1
## 1158 NW_019211456.1_12892547             0/0
## 1159 NW_019211456.1_12893809             0/0
## 1160 NW_019211456.1_13032008             0/1
## 1161 NW_019211456.1_13066662            <NA>
## 1162 NW_019211456.1_13659623            <NA>
## 1163 NW_019211456.1_13681742             0/0
## 1164 NW_019211456.1_15021677             0/0
## 1165 NW_019211456.1_15035896            <NA>
## 1166 NW_019211456.1_15048085            <NA>
## 1167 NW_019211456.1_15048146            <NA>
## 1168 NW_019211456.1_15053283            <NA>
## 1169 NW_019211456.1_15054416             0/0
## 1170 NW_019211456.1_15055335             0/0
## 1171 NW_019211456.1_15058545            <NA>
## 1172 NW_019211456.1_15061215            <NA>
## 1173 NW_019211456.1_15093465             0/0
## 1174 NW_019211456.1_15103042            <NA>
## 1175 NW_019211456.1_15105080            <NA>
## 1176 NW_019211456.1_15106684             0/0
## 1177 NW_019211456.1_15139477            <NA>
## 1178 NW_019211456.1_15141612            <NA>
## 1179 NW_019211456.1_15165693             0/0
## 1180 NW_019211456.1_15201102            <NA>
## 1181 NW_019211456.1_15207348            <NA>
## 1182 NW_019211456.1_15224287             0/1
## 1183 NW_019211456.1_15237881             0/0
## 1184 NW_019211456.1_15243567             0/0
## 1185 NW_019211456.1_15243642            <NA>
## 1186 NW_019211456.1_15261402             0/0
## 1187 NW_019211456.1_16320717            <NA>
## 1188 NW_019211456.1_16605371             0/0
## 1189 NW_019211456.1_16625728             0/0
## 1190 NW_019211456.1_17066297             0/0
## 1191 NW_019211456.1_17069158            <NA>
## 1192 NW_019211456.1_17085080            <NA>
## 1193 NW_019211456.1_17089494             0/0
## 1194 NW_019211456.1_17093759            <NA>
## 1195 NW_019211456.1_17094459             0/0
## 1196 NW_019211456.1_17107652            <NA>
## 1197 NW_019211456.1_17196021             0/0
## 1198 NW_019211456.1_17206444             0/0
## 1199 NW_019211456.1_17206854             0/0
## 1200 NW_019211456.1_17216707            <NA>
## 1201 NW_019211456.1_17216735            <NA>
## 1202 NW_019211456.1_17230296             0/0
## 1203 NW_019211456.1_17238107             0/0
## 1204 NW_019211456.1_17271127            <NA>
## 1205 NW_019211456.1_17272472             0/0
## 1206 NW_019211456.1_17303572             0/0
## 1207 NW_019211456.1_17303779             0/0
## 1208 NW_019211456.1_18084169             0/0
## 1209 NW_019211456.1_18088405             0/0
## 1210 NW_019211456.1_18091807             0/0
## 1211 NW_019211456.1_18094347            <NA>
## 1212 NW_019211456.1_18097440             0/0
## 1213 NW_019211456.1_18106795            <NA>
## 1214 NW_019211456.1_18109673             0/0
## 1215 NW_019211456.1_18110081            <NA>
## 1216 NW_019211456.1_18111667            <NA>
## 1217 NW_019211456.1_18117542            <NA>
## 1218 NW_019211456.1_18129541             0/0
## 1219 NW_019211456.1_18141664             0/0
## 1220 NW_019211456.1_18142270            <NA>
## 1221 NW_019211456.1_18168105             0/0
## 1222 NW_019211456.1_18168755            <NA>
## 1223 NW_019211456.1_18177475             0/0
## 1224 NW_019211456.1_18192142             0/0
## 1225 NW_019211456.1_18207954             0/0
## 1226 NW_019211456.1_18242958             0/0
## 1227 NW_019211456.1_18254612             0/0
## 1228 NW_019211456.1_18257798            <NA>
## 1229 NW_019211456.1_18268466             0/0
## 1230 NW_019211456.1_18284274            <NA>
## 1231 NW_019211456.1_18288521             0/0
## 1232 NW_019211456.1_18348868            <NA>
## 1233 NW_019211456.1_18349330            <NA>
## 1234 NW_019211456.1_18352718             0/0
## 1235 NW_019211456.1_18357740             0/0
## 1236 NW_019211456.1_18360835            <NA>
## 1237 NW_019211456.1_18377185            <NA>
## 1238 NW_019211456.1_18391231             0/0
## 1239 NW_019211456.1_18405490             0/0
## 1240 NW_019211456.1_18406228             0/0
## 1241 NW_019211456.1_18406565             0/0
## 1242 NW_019211456.1_18409154             0/0
## 1243 NW_019211456.1_18435890            <NA>
## 1244 NW_019211456.1_18452545            <NA>
## 1245 NW_019211456.1_18456498             0/0
## 1246 NW_019211456.1_18456594             0/0
## 1247 NW_019211456.1_18474923            <NA>
## 1248 NW_019211456.1_18476816            <NA>
## 1249 NW_019211456.1_18481702             0/0
## 1250 NW_019211456.1_18481713             0/0
## 1251 NW_019211456.1_18485302             0/0
## 1252 NW_019211456.1_18502236             0/0
## 1253 NW_019211456.1_18508056             0/0
## 1254 NW_019211456.1_18524257             0/0
## 1255 NW_019211456.1_18526248             0/0
## 1256 NW_019211456.1_18527486            <NA>
## 1257 NW_019211456.1_18541153            <NA>
## 1258 NW_019211456.1_18545531            <NA>
## 1259 NW_019211456.1_18546450            <NA>
## 1260 NW_019211456.1_18577461             0/1
## 1261 NW_019211456.1_18577756             0/0
## 1262 NW_019211456.1_18577759             0/0
## 1263 NW_019211456.1_19439475             0/0
## 1264 NW_019211456.1_20181164             0/0
## 1265 NW_019211456.1_20221719            <NA>
## 1266 NW_019211456.1_20225265             0/0
## 1267 NW_019211456.1_20245567             0/0
## 1268 NW_019211456.1_20262180            <NA>
## 1269 NW_019211456.1_20635312             0/0
## 1270 NW_019211456.1_20723040             0/0
## 1271 NW_019211456.1_20731173             0/0
## 1272 NW_019211456.1_20738096            <NA>
## 1273 NW_019211456.1_20738622             0/0
## 1274 NW_019211456.1_20751785             0/0
## 1275 NW_019211456.1_20907700             0/0
## 1276 NW_019211456.1_20925775             0/0
## 1277 NW_019211456.1_20931498             0/0
## 1278 NW_019211456.1_21507157             0/0
## 1279 NW_019211456.1_21539353             0/0
## 1280 NW_019211456.1_21544242             0/0
## 1281 NW_019211456.1_21548005             0/0
## 1282 NW_019211456.1_21578834            <NA>
## 1283 NW_019211456.1_21582779            <NA>
## 1284 NW_019211456.1_21590017            <NA>
## 1285 NW_019211456.1_21623040            <NA>
## 1286 NW_019211456.1_21623240            <NA>
## 1287 NW_019211456.1_22371006            <NA>
## 1288 NW_019211456.1_22845656            <NA>
## 1289 NW_019211456.1_22921703             0/0
## 1290 NW_019211456.1_22928635            <NA>
## 1291 NW_019211456.1_22943032             0/0
## 1292 NW_019211456.1_22962417             0/0
## 1293 NW_019211456.1_22980371            <NA>
## 1294 NW_019211456.1_23001515             0/0
## 1295 NW_019211456.1_23034338             0/0
## 1296 NW_019211456.1_23043887             0/0
## 1297 NW_019211456.1_23062534            <NA>
## 1298 NW_019211456.1_23089635             0/1
## 1299 NW_019211456.1_23116200            <NA>
## 1300 NW_019211456.1_23126151            <NA>
## 1301 NW_019211456.1_23140667            <NA>
## 1302 NW_019211456.1_23185651             0/0
## 1303 NW_019211456.1_23231210             0/0
## 1304 NW_019211456.1_23240943            <NA>
## 1305 NW_019211456.1_23247083             0/0
## 1306 NW_019211456.1_23283676            <NA>
## 1307 NW_019211456.1_23283737            <NA>
## 1308 NW_019211456.1_23284592             0/0
## 1309 NW_019211456.1_23310357             0/0
## 1310 NW_019211456.1_23387640             0/0
## 1311 NW_019211456.1_23388875            <NA>
## 1312 NW_019211456.1_23389214             0/0
## 1313 NW_019211456.1_23392562             0/0
## 1314 NW_019211456.1_23392934             0/0
## 1315 NW_019211456.1_23392953             0/0
## 1316 NW_019211456.1_23393223             0/0
## 1317 NW_019211456.1_23393805             0/0
## 1318 NW_019211456.1_23397535            <NA>
## 1319 NW_019211456.1_23398890            <NA>
## 1320 NW_019211456.1_23413176            <NA>
## 1321 NW_019211456.1_23413427             0/0
## 1322 NW_019211456.1_23413501             0/0
## 1323 NW_019211456.1_23432825             0/0
## 1324 NW_019211456.1_23447346             0/0
## 1325 NW_019211456.1_23503305            <NA>
## 1326 NW_019211456.1_23548980             0/0
## 1327 NW_019211456.1_23556990            <NA>
## 1328 NW_019211456.1_23569363             0/0
## 1329 NW_019211456.1_24801804             0/0
## 1330 NW_019211456.1_24825156            <NA>
## 1331 NW_019211456.1_24828112             0/0
## 1332 NW_019211456.1_24853255             0/0
## 1333 NW_019211456.1_24862822             0/0
## 1334 NW_019211456.1_24870581            <NA>
## 1335 NW_019211456.1_24876402            <NA>
## 1336 NW_019211456.1_24882138             0/0
## 1337 NW_019211456.1_24883429             0/0
## 1338 NW_019211456.1_24886256             0/0
## 1339 NW_019211456.1_24886443             0/0
## 1340 NW_019211456.1_24888644             0/0
## 1341 NW_019211456.1_24889364             0/0
## 1342 NW_019211456.1_24892386             0/0
## 1343 NW_019211456.1_24896884             0/0
## 1344 NW_019211456.1_24896930             0/0
## 1345 NW_019211456.1_24900777            <NA>
## 1346 NW_019211456.1_25978208            <NA>
## 1347 NW_019211456.1_25978248            <NA>
## 1348 NW_019211456.1_27009527             0/1
## 1349 NW_019211456.1_27186402             0/0
## 1350 NW_019211456.1_27219976            <NA>
## 1351 NW_019211456.1_28608466             0/0
## 1352 NW_019211456.1_28608768             0/0
## 1353 NW_019211456.1_28608973             0/0
## 1354 NW_019211456.1_28623307             0/0
## 1355 NW_019211456.1_28638661            <NA>
## 1356 NW_019211456.1_28679024            <NA>
## 1357 NW_019211456.1_28681580            <NA>
## 1358 NW_019211456.1_28684214            <NA>
## 1359 NW_019211456.1_28698184             0/0
## 1360 NW_019211456.1_28702632            <NA>
## 1361 NW_019211456.1_28705307             0/0
## 1362 NW_019211456.1_28709211            <NA>
## 1363 NW_019211456.1_28709934             0/0
## 1364 NW_019211456.1_28713209            <NA>
## 1365 NW_019211456.1_28738435            <NA>
## 1366 NW_019211456.1_28746482            <NA>
## 1367 NW_019211456.1_28756949            <NA>
## 1368 NW_019211456.1_28760675            <NA>
## 1369 NW_019211456.1_28771652            <NA>
## 1370 NW_019211456.1_28774982            <NA>
## 1371 NW_019211456.1_28775099            <NA>
## 1372 NW_019211456.1_28775116            <NA>
## 1373 NW_019211456.1_28776037            <NA>
## 1374 NW_019211456.1_28776666             0/0
## 1375 NW_019211456.1_28783428             0/0
## 1376 NW_019211456.1_28810341             0/0
## 1377 NW_019211456.1_28831177            <NA>
## 1378 NW_019211456.1_28835793            <NA>
## 1379 NW_019211456.1_29148731             0/0
## 1380 NW_019211456.1_29185663            <NA>
## 1381 NW_019211456.1_29851529             0/0
## 1382 NW_019211456.1_32158684            <NA>
## 1383 NW_019211456.1_40655890             0/0
## 1384 NW_019211456.1_41613285             0/0
## 1385 NW_019211456.1_42335654             0/0
## 1386 NW_019211456.1_42343263            <NA>
## 1387 NW_019211456.1_42347209            <NA>
## 1388 NW_019211456.1_42347330             0/0
## 1389 NW_019211456.1_42349826             0/0
## 1390 NW_019211456.1_42349837             0/0
## 1391 NW_019211456.1_42350229             0/0
## 1392 NW_019211456.1_42367184             0/0
## 1393 NW_019211456.1_42377338            <NA>
## 1394 NW_019211456.1_42380344             0/0
## 1395 NW_019211456.1_42382732            <NA>
## 1396 NW_019211456.1_42392011            <NA>
## 1397 NW_019211456.1_42400142             0/0
## 1398 NW_019211456.1_42412002             0/0
## 1399 NW_019211456.1_42412302            <NA>
## 1400 NW_019211456.1_42412364            <NA>
## 1401 NW_019211456.1_42412490             0/0
## 1402 NW_019211456.1_42413110            <NA>
## 1403 NW_019211456.1_42414609            <NA>
## 1404 NW_019211456.1_42415812             0/0
## 1405 NW_019211456.1_42422591             0/0
## 1406 NW_019211456.1_42476727             0/0
## 1407 NW_019211456.1_42482262             0/0
## 1408 NW_019211456.1_42485133             0/0
## 1409 NW_019211456.1_42485611             0/0
## 1410 NW_019211456.1_42515426            <NA>
## 1411 NW_019211456.1_42518980             0/0
## 1412 NW_019211456.1_42520527            <NA>
## 1413 NW_019211456.1_42522460            <NA>
## 1414 NW_019211456.1_42526371            <NA>
## 1415 NW_019211456.1_42526580            <NA>
## 1416 NW_019211456.1_42528120             0/0
## 1417 NW_019211456.1_42537823            <NA>
## 1418 NW_019211456.1_42541979            <NA>
## 1419 NW_019211456.1_42993189            <NA>
## 1420 NW_019211456.1_42993307            <NA>
## 1421 NW_019211456.1_45175731             0/0
## 1422 NW_019211456.1_45184752            <NA>
## 1423 NW_019211456.1_45185303            <NA>
## 1424 NW_019211456.1_45185562             0/0
## 1425 NW_019211456.1_45192221             0/0
## 1426 NW_019211456.1_45202115             0/0
## 1427 NW_019211456.1_45208993             0/0
## 1428 NW_019211456.1_45209008             0/0
## 1429 NW_019211456.1_45491709            <NA>
## 1430 NW_019211456.1_47124660             0/0
## 1431 NW_019211456.1_47125448            <NA>
## 1432 NW_019211456.1_47141944             0/0
## 1433 NW_019211456.1_47155456             0/0
## 1434 NW_019211456.1_47165647             0/0
## 1435 NW_019211456.1_47186039             0/0
## 1436 NW_019211456.1_47263170             0/0
## 1437 NW_019211456.1_47282197             0/0
## 1438 NW_019211456.1_47294803             0/0
## 1439 NW_019211456.1_47301651             0/0
## 1440 NW_019211456.1_47328092             0/0
## 1441 NW_019211456.1_47340652             0/0
## 1442 NW_019211456.1_47343892            <NA>
## 1443 NW_019211456.1_47346605            <NA>
## 1444 NW_019211456.1_47349062            <NA>
## 1445 NW_019211456.1_47370081            <NA>
## 1446 NW_019211456.1_47375562            <NA>
## 1447 NW_019211456.1_47409520             0/0
## 1448 NW_019211456.1_47421842             0/0
## 1449 NW_019211456.1_47469223             0/0
## 1450 NW_019211456.1_47475278             0/0
## 1451 NW_019211456.1_47483686            <NA>
## 1452 NW_019211456.1_47483847             0/0
## 1453 NW_019211456.1_47488539             0/0
## 1454 NW_019211456.1_47508413             0/0
## 1455 NW_019211456.1_47514026            <NA>
## 1456 NW_019211456.1_47514172             0/0
## 1457 NW_019211456.1_47518355             0/0
## 1458 NW_019211456.1_47519263             0/0
## 1459 NW_019211456.1_47520838            <NA>
## 1460 NW_019211456.1_47535679             0/0
## 1461 NW_019211456.1_47545095             0/0
## 1462 NW_019211456.1_47552379             0/0
## 1463 NW_019211456.1_47574407             0/1
## 1464 NW_019211456.1_47588360             0/0
## 1465 NW_019211456.1_47591931             0/0
## 1466 NW_019211456.1_47604128             0/0
## 1467 NW_019211456.1_47627929             0/0
## 1468 NW_019211456.1_47629543             0/0
## 1469 NW_019211456.1_47638249             0/0
## 1470 NW_019211456.1_47644566             0/0
## 1471 NW_019211456.1_47651302             0/0
## 1472 NW_019211456.1_47651579             0/0
## 1473 NW_019211456.1_47651788             0/0
## 1474 NW_019211456.1_47657124             0/0
## 1475 NW_019211456.1_47659056            <NA>
## 1476 NW_019211456.1_47663936             0/0
## 1477 NW_019211456.1_47666195            <NA>
## 1478 NW_019211456.1_47692273             0/0
## 1479 NW_019211456.1_47731162             0/0
## 1480 NW_019211456.1_47754687            <NA>
## 1481 NW_019211456.1_47761819             0/0
## 1482 NW_019211456.1_47801183             0/0
## 1483 NW_019211456.1_47804982            <NA>
## 1484 NW_019211456.1_47805594            <NA>
## 1485 NW_019211456.1_47833970            <NA>
## 1486 NW_019211456.1_47833974            <NA>
## 1487 NW_019211456.1_47861589             0/0
## 1488 NW_019211456.1_47869361             0/0
## 1489 NW_019211456.1_47874421             0/0
## 1490 NW_019211456.1_47875404            <NA>
## 1491 NW_019211456.1_47877435            <NA>
## 1492 NW_019211456.1_47886430             0/0
## 1493 NW_019211456.1_47893392            <NA>
## 1494 NW_019211456.1_47895141            <NA>
## 1495 NW_019211456.1_47896588            <NA>
## 1496 NW_019211456.1_47900582            <NA>
## 1497 NW_019211456.1_47902566             0/0
## 1498 NW_019211456.1_47908621            <NA>
## 1499 NW_019211456.1_47919263             0/0
## 1500 NW_019211456.1_47920455            <NA>
## 1501 NW_019211456.1_47920705            <NA>
## 1502 NW_019211456.1_47927215            <NA>
## 1503 NW_019211456.1_48371353             0/0
## 1504 NW_019211456.1_48387305             0/0
## 1505 NW_019211456.1_48387575            <NA>
## 1506 NW_019211456.1_49269730            <NA>
## 1507 NW_019211456.1_49279645             0/0
## 1508 NW_019211456.1_49282315            <NA>
## 1509 NW_019211456.1_49286349            <NA>
## 1510 NW_019211456.1_49288167            <NA>
## 1511 NW_019211456.1_49418838             0/0
## 1512 NW_019211456.1_49448825            <NA>
## 1513 NW_019211456.1_49449249            <NA>
## 1514 NW_019211456.1_49505240             0/0
## 1515 NW_019211456.1_49701065            <NA>
## 1516 NW_019211456.1_49758609             0/0
## 1517 NW_019211456.1_49831647             0/0
## 1518 NW_019211456.1_50325913             0/0
## 1519 NW_019211456.1_51947413            <NA>
## 1520 NW_019211456.1_54333703            <NA>
## 1521   NW_019211457.1_133204             0/0
## 1522  NW_019211457.1_1075332             0/0
## 1523  NW_019211457.1_1087768            <NA>
## 1524  NW_019211457.1_1094772             0/0
## 1525  NW_019211457.1_1123705            <NA>
## 1526  NW_019211457.1_3268466            <NA>
## 1527  NW_019211457.1_5390569            <NA>
## 1528  NW_019211457.1_5424203             0/0
## 1529  NW_019211457.1_5435145            <NA>
## 1530  NW_019211457.1_5435387             0/0
## 1531  NW_019211457.1_5439526            <NA>
## 1532  NW_019211457.1_5444660            <NA>
## 1533  NW_019211457.1_5446232             0/0
## 1534  NW_019211457.1_5447333             0/0
## 1535  NW_019211457.1_5461792            <NA>
## 1536  NW_019211457.1_5487400             0/0
## 1537  NW_019211457.1_5529960            <NA>
## 1538  NW_019211457.1_5539256             0/0
## 1539  NW_019211457.1_5540649             0/0
## 1540  NW_019211457.1_5547972             0/0
## 1541  NW_019211457.1_5563980             0/0
## 1542  NW_019211457.1_5605089            <NA>
## 1543  NW_019211457.1_5631514             0/0
## 1544  NW_019211457.1_5633509            <NA>
## 1545  NW_019211457.1_5638034             0/0
## 1546  NW_019211457.1_5690429            <NA>
## 1547  NW_019211457.1_5699077            <NA>
## 1548  NW_019211457.1_5700088             0/0
## 1549  NW_019211457.1_5703761            <NA>
## 1550  NW_019211457.1_5704004            <NA>
## 1551  NW_019211457.1_5704332            <NA>
## 1552  NW_019211457.1_7355710            <NA>
## 1553  NW_019211457.1_7379365             0/0
## 1554  NW_019211457.1_7395096             0/0
## 1555  NW_019211457.1_7422731            <NA>
## 1556  NW_019211457.1_7423898            <NA>
## 1557  NW_019211457.1_7430860            <NA>
## 1558  NW_019211457.1_7439215             0/0
## 1559  NW_019211457.1_7455656            <NA>
## 1560  NW_019211457.1_7455936             0/0
## 1561  NW_019211457.1_7462242            <NA>
## 1562  NW_019211457.1_7462527             0/0
## 1563  NW_019211457.1_7471753             0/0
## 1564  NW_019211457.1_7483869            <NA>
## 1565  NW_019211457.1_7494182            <NA>
## 1566  NW_019211457.1_7494351            <NA>
## 1567  NW_019211457.1_7509547             0/0
## 1568  NW_019211457.1_7547221            <NA>
## 1569  NW_019211457.1_7548136             0/0
## 1570  NW_019211457.1_7567575            <NA>
## 1571  NW_019211457.1_7568075             0/0
## 1572  NW_019211457.1_7570857            <NA>
## 1573  NW_019211457.1_7571443            <NA>
## 1574  NW_019211457.1_7578642             0/0
## 1575  NW_019211457.1_7580661            <NA>
## 1576  NW_019211457.1_7580666            <NA>
## 1577  NW_019211457.1_9628162             0/0
## 1578 NW_019211457.1_11930151            <NA>
## 1579 NW_019211457.1_12787763             0/0
## 1580 NW_019211457.1_13880645            <NA>
## 1581 NW_019211457.1_16383889            <NA>
## 1582 NW_019211457.1_16386363             0/0
## 1583 NW_019211457.1_16387199            <NA>
## 1584 NW_019211457.1_16387259             0/0
## 1585 NW_019211457.1_16394617             0/0
## 1586 NW_019211457.1_16402138             0/0
## 1587 NW_019211457.1_16410263             0/0
## 1588 NW_019211457.1_16419533            <NA>
## 1589 NW_019211457.1_16429444             0/0
## 1590 NW_019211457.1_16432326            <NA>
## 1591 NW_019211457.1_16458017            <NA>
## 1592 NW_019211457.1_16477347            <NA>
## 1593 NW_019211457.1_16484039             0/0
## 1594 NW_019211457.1_16497602             0/0
## 1595 NW_019211457.1_16501696            <NA>
## 1596 NW_019211457.1_16514827            <NA>
## 1597 NW_019211457.1_16701481            <NA>
## 1598 NW_019211457.1_16866355             0/0
## 1599 NW_019211457.1_17203852             0/0
## 1600 NW_019211457.1_17424004            <NA>
## 1601 NW_019211457.1_17424511            <NA>
## 1602 NW_019211457.1_17429169             0/0
## 1603 NW_019211457.1_17436733            <NA>
## 1604 NW_019211457.1_18658856             0/0
## 1605 NW_019211457.1_18688856             0/0
## 1606 NW_019211457.1_18710827            <NA>
## 1607 NW_019211457.1_18732148             0/0
## 1608 NW_019211457.1_18745006             0/0
## 1609 NW_019211457.1_18759511             0/0
## 1610 NW_019211457.1_18777241            <NA>
## 1611 NW_019211457.1_18820754             0/0
## 1612 NW_019211457.1_18834672            <NA>
## 1613 NW_019211457.1_18856501            <NA>
## 1614 NW_019211457.1_18916077             0/0
## 1615 NW_019211457.1_18942418            <NA>
## 1616 NW_019211457.1_18955864             0/0
## 1617 NW_019211457.1_18967159             0/0
## 1618 NW_019211457.1_18976389             0/0
## 1619 NW_019211457.1_19394120             0/0
## 1620 NW_019211457.1_19455436             0/0
## 1621 NW_019211457.1_19481462             0/0
## 1622 NW_019211457.1_19494042             0/0
## 1623 NW_019211457.1_19494694             0/0
## 1624 NW_019211457.1_19497539             0/0
## 1625 NW_019211457.1_19503800            <NA>
## 1626 NW_019211457.1_21646715            <NA>
## 1627 NW_019211457.1_22464706            <NA>
## 1628 NW_019211457.1_23454754             0/0
## 1629 NW_019211457.1_24968168             0/0
## 1630 NW_019211457.1_25588770             0/0
## 1631 NW_019211457.1_25606378             0/0
## 1632 NW_019211457.1_25608425             0/0
## 1633 NW_019211457.1_27909898            <NA>
## 1634 NW_019211457.1_30253563             0/0
## 1635 NW_019211457.1_30257390            <NA>
## 1636 NW_019211457.1_30267423            <NA>
## 1637 NW_019211457.1_30311075             0/0
## 1638 NW_019211457.1_30380288            <NA>
## 1639 NW_019211457.1_30389811            <NA>
## 1640 NW_019211457.1_30432550             0/0
## 1641 NW_019211457.1_30435284             0/0
## 1642 NW_019211457.1_30439161            <NA>
## 1643 NW_019211457.1_30442378            <NA>
## 1644 NW_019211457.1_30444685             0/0
## 1645 NW_019211457.1_34163229            <NA>
## 1646 NW_019211457.1_34519920             0/0
## 1647 NW_019211457.1_34548550             0/0
## 1648 NW_019211457.1_34555605             0/0
## 1649 NW_019211457.1_34573863            <NA>
## 1650 NW_019211457.1_34618145             0/0
## 1651 NW_019211457.1_34683117             0/0
## 1652 NW_019211457.1_34775151             0/0
## 1653 NW_019211457.1_34815741            <NA>
## 1654 NW_019211457.1_34817873             0/0
## 1655 NW_019211457.1_34823714             0/0
## 1656 NW_019211457.1_34828332             0/0
## 1657 NW_019211457.1_34845237             0/0
## 1658 NW_019211457.1_34856510            <NA>
## 1659 NW_019211457.1_34857001            <NA>
## 1660 NW_019211457.1_34862112             0/0
## 1661 NW_019211457.1_34863195            <NA>
## 1662 NW_019211457.1_34864568             0/0
## 1663 NW_019211457.1_37055534            <NA>
## 1664 NW_019211457.1_37061654             0/0
## 1665 NW_019211457.1_37073286            <NA>
## 1666 NW_019211457.1_37089049            <NA>
## 1667 NW_019211457.1_37091343            <NA>
## 1668 NW_019211457.1_38131790             0/0
## 1669 NW_019211457.1_38148303             0/0
## 1670 NW_019211457.1_38164723             0/0
## 1671 NW_019211457.1_38171473             0/0
## 1672 NW_019211457.1_38171511             0/0
## 1673 NW_019211457.1_38365755             0/0
## 1674 NW_019211457.1_38381478             0/0
## 1675 NW_019211457.1_38382033            <NA>
## 1676 NW_019211457.1_38399758             0/0
## 1677 NW_019211457.1_38399807             0/0
## 1678 NW_019211457.1_38410188             0/0
## 1679 NW_019211457.1_38443887             0/0
## 1680 NW_019211457.1_38463182             0/0
## 1681 NW_019211457.1_38654930            <NA>
## 1682 NW_019211457.1_39381290             0/0
## 1683 NW_019211457.1_39616341             0/0
## 1684 NW_019211457.1_39636521             0/0
## 1685 NW_019211457.1_39651693             0/0
## 1686 NW_019211457.1_39696486            <NA>
## 1687 NW_019211457.1_39958419            <NA>
## 1688 NW_019211457.1_39958571            <NA>
## 1689 NW_019211457.1_39996121             0/0
## 1690 NW_019211457.1_40102951             0/0
## 1691 NW_019211457.1_40260277            <NA>
## 1692 NW_019211457.1_40363203             0/0
## 1693 NW_019211457.1_40368698             0/0
## 1694 NW_019211457.1_40430968            <NA>
## 1695 NW_019211457.1_40452705             0/1
## 1696 NW_019211457.1_40475997             0/0
## 1697 NW_019211457.1_40548622             0/0
## 1698 NW_019211457.1_40565372            <NA>
## 1699 NW_019211457.1_40595641             0/0
## 1700 NW_019211457.1_40596634            <NA>
## 1701 NW_019211457.1_40672507             0/0
## 1702 NW_019211457.1_40680237            <NA>
## 1703 NW_019211457.1_40687798             0/0
## 1704 NW_019211457.1_40690400             0/0
## 1705 NW_019211457.1_40706014             0/0
## 1706 NW_019211457.1_40771663            <NA>
## 1707 NW_019211457.1_43018820            <NA>
## 1708 NW_019211457.1_43026455            <NA>
## 1709 NW_019211457.1_43027748             0/0
## 1710 NW_019211457.1_43056390             0/0
## 1711 NW_019211457.1_43062236             0/0
## 1712 NW_019211457.1_43064173            <NA>
## 1713 NW_019211457.1_43075503            <NA>
## 1714 NW_019211457.1_43082306            <NA>
## 1715 NW_019211457.1_43230758            <NA>
## 1716 NW_019211457.1_43404179            <NA>
## 1717 NW_019211457.1_43412111             0/0
## 1718 NW_019211457.1_43419956             0/0
## 1719 NW_019211457.1_43423239            <NA>
## 1720 NW_019211457.1_43423554             0/0
## 1721 NW_019211457.1_43434782            <NA>
## 1722 NW_019211457.1_43441999             0/0
## 1723 NW_019211457.1_43448690             0/0
## 1724 NW_019211457.1_43460214             0/0
## 1725 NW_019211457.1_43477460             0/0
## 1726 NW_019211457.1_43491491            <NA>
## 1727 NW_019211457.1_43494064             0/0
## 1728 NW_019211457.1_43494184             0/0
## 1729 NW_019211457.1_43500273            <NA>
## 1730 NW_019211457.1_43506468            <NA>
## 1731 NW_019211457.1_43513947             0/0
## 1732 NW_019211457.1_43522475             0/0
## 1733 NW_019211457.1_43524965             0/0
## 1734 NW_019211457.1_43560657             0/0
## 1735 NW_019211457.1_43567633             0/0
## 1736 NW_019211457.1_43570115            <NA>
## 1737 NW_019211457.1_43579126            <NA>
## 1738 NW_019211457.1_43599283             0/0
## 1739 NW_019211457.1_43624623            <NA>
## 1740 NW_019211457.1_43628278            <NA>
## 1741 NW_019211457.1_43646659             0/0
## 1742 NW_019211457.1_43664111            <NA>
## 1743 NW_019211457.1_43668917             0/0
## 1744 NW_019211457.1_43672968            <NA>
## 1745 NW_019211457.1_43688371             0/0
## 1746 NW_019211457.1_43691038             0/0
## 1747 NW_019211457.1_43692692            <NA>
## 1748 NW_019211457.1_43699438             0/0
## 1749 NW_019211457.1_43715387             0/0
## 1750 NW_019211457.1_43716304            <NA>
## 1751 NW_019211457.1_43719012             0/0
## 1752 NW_019211457.1_43734279            <NA>
## 1753 NW_019211457.1_43762754            <NA>
## 1754 NW_019211457.1_43765573             0/0
## 1755 NW_019211457.1_43769507             0/0
## 1756 NW_019211457.1_43770686            <NA>
## 1757 NW_019211457.1_43775412             0/0
## 1758 NW_019211457.1_43781953            <NA>
## 1759 NW_019211457.1_43784201             0/0
## 1760 NW_019211457.1_43798996             0/0
## 1761 NW_019211457.1_43799878            <NA>
## 1762 NW_019211457.1_43819800            <NA>
## 1763 NW_019211457.1_43829987             0/0
## 1764 NW_019211457.1_43833419            <NA>
## 1765 NW_019211457.1_43841047             0/0
## 1766 NW_019211457.1_43865442             0/0
## 1767 NW_019211457.1_43866548             0/0
## 1768 NW_019211457.1_43890042            <NA>
## 1769 NW_019211457.1_43899806            <NA>
## 1770 NW_019211457.1_43911156             0/0
## 1771 NW_019211457.1_43932093            <NA>
## 1772 NW_019211457.1_43932155            <NA>
## 1773 NW_019211457.1_43937584            <NA>
## 1774 NW_019211457.1_43946568            <NA>
## 1775 NW_019211457.1_43955715            <NA>
## 1776 NW_019211457.1_43961122             0/0
## 1777 NW_019211457.1_43988546             0/0
## 1778 NW_019211457.1_44025576            <NA>
## 1779 NW_019211457.1_44064629             0/0
## 1780 NW_019211457.1_44095649             0/0
## 1781 NW_019211457.1_44111224             0/0
## 1782 NW_019211457.1_44130335            <NA>
## 1783 NW_019211457.1_44134277            <NA>
## 1784 NW_019211457.1_44141355             0/0
## 1785 NW_019211457.1_44142885            <NA>
## 1786 NW_019211457.1_44167928             0/0
## 1787 NW_019211457.1_44178254             0/0
## 1788 NW_019211457.1_44200136             0/0
## 1789 NW_019211457.1_44235626             0/0
## 1790 NW_019211457.1_44240251            <NA>
## 1791 NW_019211457.1_44255839             0/0
## 1792 NW_019211457.1_44275899             0/0
## 1793 NW_019211457.1_44295414            <NA>
## 1794 NW_019211457.1_44295502            <NA>
## 1795 NW_019211457.1_45023343            <NA>
## 1796 NW_019211457.1_45034118             0/0
## 1797 NW_019211457.1_45035398            <NA>
## 1798 NW_019211457.1_45042444             0/0
## 1799 NW_019211457.1_45048100             0/0
## 1800 NW_019211457.1_45063030             0/0
## 1801 NW_019211457.1_45071546            <NA>
## 1802 NW_019211457.1_45073516            <NA>
## 1803 NW_019211457.1_50318127             0/0
## 1804 NW_019211457.1_50962362            <NA>
## 1805 NW_019211457.1_51443167             0/0
## 1806 NW_019211457.1_51478363             0/0
## 1807 NW_019211457.1_51513257             0/0
## 1808 NW_019211457.1_51546372            <NA>
## 1809 NW_019211457.1_51581192            <NA>
## 1810 NW_019211457.1_51584412            <NA>
## 1811 NW_019211457.1_51588475             0/0
## 1812 NW_019211457.1_51683102             0/0
## 1813 NW_019211457.1_51705140            <NA>
## 1814 NW_019211457.1_51710311             0/0
## 1815 NW_019211457.1_51761327             0/0
## 1816 NW_019211457.1_51763547             0/0
## 1817 NW_019211457.1_51769826            <NA>
## 1818 NW_019211457.1_51775782             0/0
## 1819 NW_019211457.1_51779339             0/0
## 1820 NW_019211457.1_51873098             0/0
## 1821 NW_019211457.1_51944703             0/0
## 1822 NW_019211457.1_51961139            <NA>
## 1823 NW_019211457.1_51961441            <NA>
## 1824 NW_019211457.1_51976089             0/0
## 1825 NW_019211457.1_52045948             0/0
## 1826 NW_019211457.1_52069480             0/0
## 1827 NW_019211457.1_52104315            <NA>
## 1828 NW_019211457.1_52183824            <NA>
## 1829 NW_019211457.1_52184827             0/0
## 1830 NW_019211457.1_52196943            <NA>
## 1831 NW_019211457.1_52249527            <NA>
## 1832 NW_019211457.1_52258187             0/0
## 1833 NW_019211457.1_52302222            <NA>
## 1834 NW_019211457.1_52321663             0/0
## 1835 NW_019211457.1_52353812             0/0
## 1836 NW_019211457.1_52363162            <NA>
## 1837 NW_019211457.1_52365285             0/0
## 1838 NW_019211457.1_52365614             0/0
## 1839 NW_019211457.1_52387359             0/0
## 1840 NW_019211457.1_52402281            <NA>
## 1841 NW_019211457.1_52404236             0/0
## 1842 NW_019211457.1_52417174            <NA>
## 1843 NW_019211457.1_52417275             0/0
## 1844 NW_019211457.1_52417534             0/0
## 1845 NW_019211457.1_52443043            <NA>
## 1846 NW_019211457.1_52454244            <NA>
## 1847 NW_019211457.1_52469489             0/0
## 1848 NW_019211457.1_52478758             0/0
## 1849 NW_019211457.1_52478780             0/0
## 1850 NW_019211457.1_52479560             0/0
## 1851 NW_019211457.1_52479570             0/0
## 1852 NW_019211457.1_52482140             0/0
## 1853 NW_019211457.1_52661003             0/0
## 1854  NW_019211458.1_1298382            <NA>
## 1855  NW_019211458.1_1583182            <NA>
## 1856  NW_019211458.1_2084304             0/0
## 1857  NW_019211458.1_2563007            <NA>
## 1858  NW_019211458.1_2912793            <NA>
## 1859  NW_019211458.1_3486581             0/0
## 1860  NW_019211458.1_3501170             0/0
## 1861  NW_019211458.1_3517571             0/0
## 1862  NW_019211458.1_3538441             0/0
## 1863  NW_019211458.1_3563440             0/0
## 1864  NW_019211458.1_3594149            <NA>
## 1865  NW_019211458.1_3603498            <NA>
## 1866  NW_019211458.1_3615957            <NA>
## 1867  NW_019211458.1_3632231             0/0
## 1868  NW_019211458.1_3634791             0/0
## 1869  NW_019211458.1_3637222            <NA>
## 1870  NW_019211458.1_3638327            <NA>
## 1871  NW_019211458.1_3638631            <NA>
## 1872  NW_019211458.1_3639626            <NA>
## 1873  NW_019211458.1_3640124            <NA>
## 1874  NW_019211458.1_3660510             0/0
## 1875  NW_019211458.1_3674364            <NA>
## 1876  NW_019211458.1_3677885            <NA>
## 1877  NW_019211458.1_3682080             0/0
## 1878  NW_019211458.1_3697615             0/0
## 1879  NW_019211458.1_3705441             0/0
## 1880  NW_019211458.1_3717645            <NA>
## 1881  NW_019211458.1_3761906            <NA>
## 1882  NW_019211458.1_3904036             0/0
## 1883  NW_019211458.1_3915474            <NA>
## 1884  NW_019211458.1_3932901             0/0
## 1885  NW_019211458.1_3958477            <NA>
## 1886  NW_019211458.1_3972721            <NA>
## 1887  NW_019211458.1_4030552            <NA>
## 1888  NW_019211458.1_4073158            <NA>
## 1889  NW_019211458.1_4084498            <NA>
## 1890  NW_019211458.1_4147248             0/0
## 1891  NW_019211458.1_4161593             0/0
## 1892  NW_019211458.1_4204532            <NA>
## 1893  NW_019211458.1_4237489            <NA>
## 1894  NW_019211458.1_4259639             0/0
## 1895  NW_019211458.1_4265158            <NA>
## 1896  NW_019211458.1_4321607             0/0
## 1897  NW_019211458.1_4328941             0/0
## 1898  NW_019211458.1_4394285            <NA>
## 1899  NW_019211458.1_4503041            <NA>
## 1900  NW_019211458.1_4537007             0/0
## 1901  NW_019211458.1_4561542            <NA>
## 1902  NW_019211458.1_4562049            <NA>
## 1903  NW_019211458.1_4567805             0/0
## 1904  NW_019211458.1_4599456            <NA>
## 1905  NW_019211458.1_4601710             0/0
## 1906  NW_019211458.1_4718955            <NA>
## 1907  NW_019211458.1_4731554            <NA>
## 1908  NW_019211458.1_4736124            <NA>
## 1909  NW_019211458.1_4767724            <NA>
## 1910  NW_019211458.1_4791848            <NA>
## 1911  NW_019211458.1_4807639            <NA>
## 1912  NW_019211458.1_4853777            <NA>
## 1913  NW_019211458.1_4859954             0/0
## 1914  NW_019211458.1_4864440             0/0
## 1915  NW_019211458.1_4873032             0/0
## 1916  NW_019211458.1_4925166             0/0
## 1917  NW_019211458.1_5039431            <NA>
## 1918  NW_019211458.1_5044338             0/0
## 1919  NW_019211458.1_5063604             1/1
## 1920  NW_019211458.1_5086921            <NA>
## 1921  NW_019211458.1_5110485             0/0
## 1922  NW_019211458.1_5118361            <NA>
## 1923  NW_019211458.1_5128990            <NA>
## 1924  NW_019211458.1_5173802             0/0
## 1925  NW_019211458.1_5192566             0/0
## 1926  NW_019211458.1_5192913            <NA>
## 1927  NW_019211458.1_5199497             0/0
## 1928  NW_019211458.1_5269581             0/0
## 1929  NW_019211458.1_5299536             0/0
## 1930  NW_019211458.1_5299674            <NA>
## 1931  NW_019211458.1_5307804            <NA>
## 1932  NW_019211458.1_5374857             0/0
## 1933  NW_019211458.1_5479501             0/0
## 1934  NW_019211458.1_5494423            <NA>
## 1935  NW_019211458.1_5497967            <NA>
## 1936  NW_019211458.1_5503263             0/0
## 1937  NW_019211458.1_5504618             0/0
## 1938  NW_019211458.1_5511453             0/0
## 1939  NW_019211458.1_5523299            <NA>
## 1940  NW_019211458.1_5531127            <NA>
## 1941  NW_019211458.1_5535808             0/0
## 1942  NW_019211458.1_5546803             0/0
## 1943  NW_019211458.1_5547623             0/0
## 1944  NW_019211458.1_5557778             0/0
## 1945  NW_019211458.1_5574279             0/0
## 1946  NW_019211458.1_5577966            <NA>
## 1947  NW_019211458.1_5595097            <NA>
## 1948  NW_019211458.1_5600497             0/0
## 1949  NW_019211458.1_5600613            <NA>
## 1950  NW_019211458.1_5613882             0/0
## 1951  NW_019211458.1_5631738            <NA>
## 1952  NW_019211458.1_5680791             0/0
## 1953  NW_019211458.1_5694042            <NA>
## 1954  NW_019211458.1_5722833             0/0
## 1955  NW_019211458.1_5741658             0/0
## 1956  NW_019211458.1_5741703            <NA>
## 1957  NW_019211458.1_5804801            <NA>
## 1958  NW_019211458.1_5804897            <NA>
## 1959  NW_019211458.1_5823319            <NA>
## 1960  NW_019211458.1_5847489            <NA>
## 1961  NW_019211458.1_5903748             0/0
## 1962  NW_019211458.1_5912650            <NA>
## 1963  NW_019211458.1_5914918             0/0
## 1964  NW_019211458.1_5917146             0/0
## 1965  NW_019211458.1_5924336             0/0
## 1966  NW_019211458.1_5929745             0/0
## 1967  NW_019211458.1_5934807            <NA>
## 1968  NW_019211458.1_5935652             0/0
## 1969  NW_019211458.1_5937197             0/0
## 1970  NW_019211458.1_5938445            <NA>
## 1971  NW_019211458.1_5938601             0/0
## 1972  NW_019211458.1_5965772             0/0
## 1973  NW_019211458.1_5979869            <NA>
## 1974  NW_019211458.1_5992447             0/0
## 1975  NW_019211458.1_6009470            <NA>
## 1976  NW_019211458.1_6011168             0/0
## 1977  NW_019211458.1_6030794             0/0
## 1978  NW_019211458.1_6042570            <NA>
## 1979  NW_019211458.1_6057177             0/0
## 1980  NW_019211458.1_6058524             0/0
## 1981  NW_019211458.1_6058691            <NA>
## 1982  NW_019211458.1_6061677             0/0
## 1983  NW_019211458.1_6061843             0/0
## 1984  NW_019211458.1_6069829             0/0
## 1985  NW_019211458.1_6069908             0/0
## 1986  NW_019211458.1_6069982             0/0
## 1987  NW_019211458.1_6071819            <NA>
## 1988  NW_019211458.1_6144848            <NA>
## 1989  NW_019211458.1_6151080            <NA>
## 1990  NW_019211458.1_6191019             0/0
## 1991  NW_019211458.1_6191325             0/0
## 1992  NW_019211458.1_6202189            <NA>
## 1993  NW_019211458.1_6223063            <NA>
## 1994  NW_019211458.1_6448227             0/0
## 1995  NW_019211458.1_6565784            <NA>
## 1996  NW_019211458.1_6568869            <NA>
## 1997  NW_019211458.1_6572907             0/0
## 1998  NW_019211458.1_6576809             0/0
## 1999  NW_019211458.1_6598025             0/0
## 2000  NW_019211458.1_6600418            <NA>
## 2001  NW_019211458.1_6619620            <NA>
## 2002  NW_019211458.1_6625453             0/0
## 2003  NW_019211458.1_6640507            <NA>
## 2004  NW_019211458.1_6683613             0/0
## 2005  NW_019211458.1_6751037             0/0
## 2006  NW_019211458.1_6754083             0/0
## 2007  NW_019211458.1_6764061             0/0
## 2008  NW_019211458.1_6764939            <NA>
## 2009  NW_019211458.1_6778749            <NA>
## 2010  NW_019211458.1_6782718             0/0
## 2011  NW_019211458.1_6785629            <NA>
## 2012  NW_019211458.1_6795715            <NA>
## 2013  NW_019211458.1_6803452            <NA>
## 2014  NW_019211458.1_6811442            <NA>
## 2015  NW_019211458.1_6827923            <NA>
## 2016  NW_019211458.1_6842486             0/0
## 2017  NW_019211458.1_6885232             0/0
## 2018  NW_019211458.1_6933656            <NA>
## 2019  NW_019211458.1_6934656             0/0
## 2020  NW_019211458.1_6936668             0/0
## 2021  NW_019211458.1_6961890            <NA>
## 2022  NW_019211458.1_6961981            <NA>
## 2023  NW_019211458.1_6967444            <NA>
## 2024  NW_019211458.1_6980583             0/0
## 2025  NW_019211458.1_6981187             0/0
## 2026  NW_019211458.1_7208511             0/0
## 2027  NW_019211458.1_7235294            <NA>
## 2028  NW_019211458.1_7264863            <NA>
## 2029  NW_019211458.1_7265347             0/0
## 2030  NW_019211458.1_7631653            <NA>
## 2031  NW_019211458.1_7640452            <NA>
## 2032  NW_019211458.1_7774133             0/0
## 2033  NW_019211458.1_7786800            <NA>
## 2034  NW_019211458.1_7787879             0/0
## 2035  NW_019211458.1_7791296            <NA>
## 2036  NW_019211458.1_7799492             0/0
## 2037  NW_019211458.1_7866350             0/0
## 2038  NW_019211458.1_7888320             0/0
## 2039  NW_019211458.1_7889638             0/0
## 2040  NW_019211458.1_7954921            <NA>
## 2041  NW_019211458.1_7974357             0/0
## 2042  NW_019211458.1_7984851             0/0
## 2043  NW_019211458.1_8117696             0/0
## 2044  NW_019211458.1_8237922            <NA>
## 2045  NW_019211458.1_8239891             0/0
## 2046  NW_019211458.1_8251037            <NA>
## 2047  NW_019211458.1_8336887            <NA>
## 2048  NW_019211458.1_8369585            <NA>
## 2049  NW_019211458.1_8409414             0/0
## 2050  NW_019211458.1_8436682             0/0
## 2051  NW_019211458.1_8451711             0/0
## 2052  NW_019211458.1_8466673             0/0
## 2053  NW_019211458.1_8496190            <NA>
## 2054  NW_019211458.1_8504608            <NA>
## 2055  NW_019211458.1_8523779             0/0
## 2056  NW_019211458.1_8541409            <NA>
## 2057  NW_019211458.1_8568317            <NA>
## 2058  NW_019211458.1_8753223             0/0
## 2059  NW_019211458.1_8772588             0/0
## 2060  NW_019211458.1_8801693             0/0
## 2061  NW_019211458.1_8824576             0/0
## 2062  NW_019211458.1_8842277             0/0
## 2063  NW_019211458.1_8845053            <NA>
## 2064  NW_019211458.1_8877896             0/0
## 2065  NW_019211458.1_8887963             0/0
## 2066  NW_019211458.1_8896964             0/0
## 2067  NW_019211458.1_8898831            <NA>
## 2068  NW_019211458.1_8904369             0/0
## 2069  NW_019211458.1_8908098            <NA>
## 2070  NW_019211458.1_8908442             0/0
## 2071  NW_019211458.1_8915801            <NA>
## 2072  NW_019211458.1_8925533            <NA>
## 2073  NW_019211458.1_8929328            <NA>
## 2074  NW_019211458.1_8929373            <NA>
## 2075  NW_019211458.1_8941054             0/0
## 2076  NW_019211458.1_8954838            <NA>
## 2077 NW_019211458.1_11432208             0/0
## 2078 NW_019211458.1_11481565            <NA>
## 2079 NW_019211458.1_11481676            <NA>
## 2080 NW_019211458.1_11482583             0/0
## 2081 NW_019211458.1_11505070             0/0
## 2082 NW_019211458.1_11518442             0/0
## 2083 NW_019211458.1_11519731             0/0
## 2084 NW_019211458.1_11570514             0/0
## 2085 NW_019211458.1_11601340             0/0
## 2086 NW_019211458.1_11713555            <NA>
## 2087 NW_019211458.1_11713592            <NA>
## 2088 NW_019211458.1_11815375             0/1
## 2089 NW_019211458.1_12922078            <NA>
## 2090 NW_019211458.1_14116016            <NA>
## 2091 NW_019211458.1_15139477             0/0
## 2092 NW_019211458.1_15824535             0/0
## 2093 NW_019211458.1_15827758             0/0
## 2094 NW_019211458.1_15836665            <NA>
## 2095 NW_019211458.1_15859437             0/0
## 2096 NW_019211458.1_15878169            <NA>
## 2097 NW_019211458.1_15900717            <NA>
## 2098 NW_019211458.1_15905441            <NA>
## 2099 NW_019211458.1_15915235             0/0
## 2100 NW_019211458.1_15953469            <NA>
## 2101 NW_019211458.1_15991683             0/0
## 2102 NW_019211458.1_16001815             0/0
## 2103 NW_019211458.1_16007541             0/0
## 2104 NW_019211458.1_16023747             0/0
## 2105 NW_019211458.1_16023907             0/0
## 2106 NW_019211458.1_16024023             0/0
## 2107 NW_019211458.1_16066803            <NA>
## 2108 NW_019211458.1_16069923             0/0
## 2109 NW_019211458.1_16070177             0/0
## 2110 NW_019211458.1_16088760             0/0
## 2111 NW_019211458.1_16095165             0/0
## 2112 NW_019211458.1_16095833            <NA>
## 2113 NW_019211458.1_19013418            <NA>
## 2114 NW_019211458.1_19018068             0/0
## 2115 NW_019211458.1_19027189            <NA>
## 2116 NW_019211458.1_19028926             0/0
## 2117 NW_019211458.1_19031120             0/0
## 2118 NW_019211458.1_19050065             0/0
## 2119 NW_019211458.1_19050171            <NA>
## 2120 NW_019211458.1_19055947            <NA>
## 2121 NW_019211458.1_19066857            <NA>
## 2122 NW_019211458.1_19068543             0/0
## 2123 NW_019211458.1_19068845             0/0
## 2124 NW_019211458.1_19078895             0/0
## 2125 NW_019211458.1_19078983             0/0
## 2126 NW_019211458.1_19079095             0/0
## 2127 NW_019211458.1_19084806             0/0
## 2128 NW_019211458.1_19084873             0/0
## 2129 NW_019211458.1_19091764            <NA>
## 2130 NW_019211458.1_19092961             0/0
## 2131 NW_019211458.1_19095912            <NA>
## 2132 NW_019211458.1_19109324            <NA>
## 2133 NW_019211458.1_19113505             0/0
## 2134 NW_019211458.1_19119879            <NA>
## 2135 NW_019211458.1_19119961            <NA>
## 2136 NW_019211458.1_19124025             0/0
## 2137 NW_019211458.1_19181055             0/0
## 2138 NW_019211458.1_19181271             0/0
## 2139 NW_019211458.1_19196850            <NA>
## 2140 NW_019211458.1_19201768            <NA>
## 2141 NW_019211458.1_19215862             0/0
## 2142 NW_019211458.1_19232339             0/0
## 2143 NW_019211458.1_19253053             0/0
## 2144 NW_019211458.1_19257350            <NA>
## 2145 NW_019211458.1_19267663             0/0
## 2146 NW_019211458.1_19286737            <NA>
## 2147 NW_019211458.1_19294933             0/0
## 2148 NW_019211458.1_19338531             0/0
## 2149 NW_019211458.1_19354473             0/0
## 2150 NW_019211458.1_19354796            <NA>
## 2151 NW_019211458.1_19365119             0/0
## 2152 NW_019211458.1_19371573             0/0
## 2153 NW_019211458.1_19396211            <NA>
## 2154 NW_019211458.1_19397292             0/0
## 2155 NW_019211458.1_19400732             0/0
## 2156 NW_019211458.1_19405277             0/0
## 2157 NW_019211458.1_19405888             0/0
## 2158 NW_019211458.1_19596512            <NA>
## 2159 NW_019211458.1_19601614             0/0
## 2160 NW_019211458.1_19602033             0/0
## 2161 NW_019211458.1_19602247             0/0
## 2162 NW_019211458.1_19604009             0/0
## 2163 NW_019211458.1_19605186            <NA>
## 2164 NW_019211458.1_19616226             0/0
## 2165 NW_019211458.1_19624725            <NA>
## 2166 NW_019211458.1_19624973            <NA>
## 2167 NW_019211458.1_19626291            <NA>
## 2168 NW_019211458.1_19638190             0/0
## 2169 NW_019211458.1_19650797             0/0
## 2170 NW_019211458.1_19659026            <NA>
## 2171 NW_019211458.1_19660483             0/0
## 2172 NW_019211458.1_19660856            <NA>
## 2173 NW_019211458.1_19678845            <NA>
## 2174 NW_019211458.1_19725910             0/0
## 2175 NW_019211458.1_19730714             0/0
## 2176 NW_019211458.1_19735907             0/0
## 2177 NW_019211458.1_19737396            <NA>
## 2178 NW_019211458.1_19740580            <NA>
## 2179 NW_019211458.1_19836095             0/0
## 2180 NW_019211458.1_19837623            <NA>
## 2181 NW_019211458.1_19841141            <NA>
## 2182 NW_019211458.1_19842654            <NA>
## 2183 NW_019211458.1_19873331             0/0
## 2184 NW_019211458.1_19880995            <NA>
## 2185 NW_019211458.1_19881862            <NA>
## 2186 NW_019211458.1_19887566            <NA>
## 2187 NW_019211458.1_19901298            <NA>
## 2188 NW_019211458.1_19908020            <NA>
## 2189 NW_019211458.1_19908744            <NA>
## 2190 NW_019211458.1_19908783            <NA>
## 2191 NW_019211458.1_19920008             0/0
## 2192 NW_019211458.1_19992993             0/1
## 2193 NW_019211458.1_20838753            <NA>
## 2194 NW_019211458.1_21027260             0/0
## 2195 NW_019211458.1_27062739            <NA>
## 2196 NW_019211458.1_27189781             0/0
## 2197 NW_019211458.1_27206565             0/0
## 2198 NW_019211458.1_27211355            <NA>
## 2199 NW_019211458.1_27216651             0/0
## 2200 NW_019211458.1_27219710             0/0
## 2201 NW_019211458.1_27223068            <NA>
## 2202 NW_019211458.1_27240750             0/0
## 2203 NW_019211458.1_27258002             0/0
## 2204 NW_019211458.1_27265899             0/0
## 2205 NW_019211458.1_27270567            <NA>
## 2206 NW_019211458.1_27273649             0/0
## 2207 NW_019211458.1_27319003             0/0
## 2208 NW_019211458.1_27339853            <NA>
## 2209 NW_019211458.1_27350350            <NA>
## 2210 NW_019211458.1_27355400             0/0
## 2211 NW_019211458.1_27802083             0/0
## 2212 NW_019211458.1_27803679            <NA>
## 2213 NW_019211458.1_27913829             0/0
## 2214 NW_019211458.1_28179038             0/0
## 2215 NW_019211458.1_28183832             0/0
## 2216 NW_019211458.1_28187161            <NA>
## 2217 NW_019211458.1_28190473             0/0
## 2218 NW_019211458.1_28219707             0/0
## 2219 NW_019211458.1_28232864             0/0
## 2220 NW_019211458.1_29252933            <NA>
## 2221 NW_019211458.1_32332413            <NA>
## 2222 NW_019211458.1_32361992            <NA>
## 2223 NW_019211458.1_32398801             0/0
## 2224 NW_019211458.1_32400417            <NA>
## 2225 NW_019211458.1_32434604             0/0
## 2226 NW_019211458.1_32723983             0/0
## 2227 NW_019211458.1_32761400             0/0
## 2228 NW_019211458.1_32801232            <NA>
## 2229 NW_019211458.1_32819244             0/0
## 2230 NW_019211458.1_32824333            <NA>
## 2231 NW_019211458.1_32841501             0/0
## 2232 NW_019211458.1_32858811             0/0
## 2233 NW_019211458.1_32860345             0/0
## 2234 NW_019211458.1_32875388            <NA>
## 2235 NW_019211458.1_32879957            <NA>
## 2236 NW_019211458.1_32882940             0/0
## 2237 NW_019211458.1_33022961             0/0
## 2238 NW_019211458.1_33031453            <NA>
## 2239 NW_019211458.1_33063753             0/0
## 2240 NW_019211458.1_33064558            <NA>
## 2241 NW_019211458.1_33107893            <NA>
## 2242 NW_019211458.1_33109485            <NA>
## 2243 NW_019211458.1_33233555             0/1
## 2244 NW_019211458.1_33522104            <NA>
## 2245 NW_019211458.1_33804852             0/0
## 2246 NW_019211458.1_34302052            <NA>
## 2247 NW_019211458.1_34308469            <NA>
## 2248 NW_019211458.1_34352594            <NA>
## 2249 NW_019211458.1_34359656            <NA>
## 2250 NW_019211458.1_34365462            <NA>
## 2251 NW_019211458.1_34463012             0/0
## 2252 NW_019211458.1_35512875             0/0
## 2253 NW_019211458.1_35528359             0/0
## 2254 NW_019211458.1_35533449            <NA>
## 2255 NW_019211458.1_35540304            <NA>
## 2256 NW_019211458.1_35582307             0/0
## 2257 NW_019211458.1_35596997            <NA>
## 2258 NW_019211458.1_35598175             0/0
## 2259 NW_019211458.1_35610287             0/0
## 2260 NW_019211458.1_35761965             0/0
## 2261 NW_019211458.1_35828090            <NA>
## 2262 NW_019211458.1_36241182            <NA>
## 2263 NW_019211458.1_36242813            <NA>
## 2264 NW_019211458.1_36257868             0/0
## 2265 NW_019211458.1_36273185             0/0
## 2266 NW_019211458.1_36274829            <NA>
## 2267 NW_019211458.1_36299576             0/0
## 2268 NW_019211458.1_36323521            <NA>
## 2269 NW_019211458.1_36339994            <NA>
## 2270 NW_019211458.1_36345371             0/0
## 2271 NW_019211458.1_36346267             0/0
## 2272 NW_019211458.1_36353126             0/0
## 2273 NW_019211458.1_36354283            <NA>
## 2274 NW_019211458.1_36865715             0/0
## 2275 NW_019211458.1_36876584             0/0
## 2276 NW_019211458.1_36876619            <NA>
## 2277 NW_019211458.1_36878418             0/0
## 2278 NW_019211458.1_36913417            <NA>
## 2279 NW_019211458.1_36914408             0/0
## 2280 NW_019211458.1_36914923            <NA>
## 2281 NW_019211458.1_36915002             0/0
## 2282 NW_019211458.1_36959919            <NA>
## 2283 NW_019211458.1_36999868            <NA>
## 2284 NW_019211458.1_37022712             0/0
## 2285 NW_019211458.1_37041696            <NA>
## 2286 NW_019211458.1_37052376            <NA>
## 2287 NW_019211458.1_37073289             0/0
## 2288 NW_019211458.1_37120474             0/0
## 2289 NW_019211458.1_37217687            <NA>
## 2290 NW_019211458.1_37227564             0/0
## 2291 NW_019211458.1_37229032             0/0
## 2292 NW_019211458.1_37230357            <NA>
## 2293 NW_019211458.1_37243315             0/0
## 2294 NW_019211458.1_37283307             0/0
## 2295 NW_019211458.1_37283314            <NA>
## 2296 NW_019211458.1_37320842             0/0
## 2297 NW_019211458.1_37351519             0/0
## 2298 NW_019211458.1_37399338             0/0
## 2299 NW_019211458.1_37408537             0/0
## 2300 NW_019211458.1_37416973            <NA>
## 2301 NW_019211458.1_37423183             0/0
## 2302 NW_019211458.1_37451577            <NA>
## 2303 NW_019211458.1_37492743            <NA>
## 2304 NW_019211458.1_37557391             0/0
## 2305 NW_019211458.1_37621585            <NA>
## 2306 NW_019211458.1_37621641            <NA>
## 2307 NW_019211458.1_37679835            <NA>
## 2308 NW_019211458.1_37680170            <NA>
## 2309 NW_019211458.1_37682511             0/0
## 2310 NW_019211458.1_37684867             0/0
## 2311 NW_019211458.1_37714114            <NA>
## 2312 NW_019211458.1_37782487             0/0
## 2313 NW_019211458.1_37855379             0/0
## 2314 NW_019211458.1_37869307            <NA>
## 2315 NW_019211458.1_37903630            <NA>
## 2316 NW_019211458.1_37908386            <NA>
## 2317 NW_019211458.1_37950222             0/0
## 2318 NW_019211458.1_37953562             0/0
## 2319 NW_019211458.1_37960489             0/0
## 2320 NW_019211458.1_37978174            <NA>
## 2321 NW_019211458.1_38021076             0/0
## 2322 NW_019211458.1_38023027             0/0
## 2323 NW_019211458.1_38023475             0/0
## 2324 NW_019211458.1_38024944            <NA>
## 2325 NW_019211458.1_38025836             0/0
## 2326 NW_019211458.1_38027113            <NA>
## 2327 NW_019211458.1_38520712            <NA>
## 2328 NW_019211458.1_38538361             0/0
## 2329 NW_019211458.1_38548438            <NA>
## 2330 NW_019211458.1_38559480             0/0
## 2331 NW_019211458.1_38563467             0/0
## 2332 NW_019211458.1_38570601             0/0
## 2333 NW_019211458.1_38572714             0/0
## 2334 NW_019211458.1_38582213            <NA>
## 2335 NW_019211458.1_38590371             0/0
## 2336 NW_019211458.1_38636482            <NA>
## 2337 NW_019211458.1_38661992            <NA>
## 2338 NW_019211458.1_38699410             0/0
## 2339 NW_019211458.1_38719109             0/0
## 2340 NW_019211458.1_38932488             0/0
## 2341 NW_019211458.1_38932501             0/0
## 2342 NW_019211458.1_38932509             0/0
## 2343 NW_019211458.1_40638176             0/0
## 2344 NW_019211458.1_41788723            <NA>
## 2345 NW_019211458.1_41788885             0/0
## 2346 NW_019211458.1_41791574             0/0
## 2347 NW_019211458.1_41796965             0/0
## 2348 NW_019211458.1_41804698             0/0
## 2349 NW_019211458.1_41807060             0/0
## 2350 NW_019211458.1_41811488             0/0
## 2351 NW_019211458.1_41821532             0/0
## 2352 NW_019211458.1_41826848             0/0
## 2353 NW_019211458.1_41830093             0/0
## 2354 NW_019211458.1_41839896            <NA>
## 2355 NW_019211458.1_41887673            <NA>
## 2356   NW_019211459.1_510323             0/0
## 2357   NW_019211459.1_510350             0/0
## 2358   NW_019211459.1_532411            <NA>
## 2359   NW_019211459.1_624725            <NA>
## 2360   NW_019211459.1_632570             0/0
## 2361  NW_019211459.1_1347466            <NA>
## 2362  NW_019211459.1_1367453            <NA>
## 2363  NW_019211459.1_1700114            <NA>
## 2364  NW_019211459.1_1939843            <NA>
## 2365  NW_019211459.1_1941683            <NA>
## 2366  NW_019211459.1_2519543             0/0
## 2367  NW_019211459.1_2524265            <NA>
## 2368  NW_019211459.1_2542932             0/0
## 2369  NW_019211459.1_2574745             0/0
## 2370  NW_019211459.1_2598829             0/0
## 2371  NW_019211459.1_2610146            <NA>
## 2372  NW_019211459.1_2653117             0/0
## 2373  NW_019211459.1_2686546            <NA>
## 2374  NW_019211459.1_2693033             0/0
## 2375  NW_019211459.1_2696228             0/0
## 2376  NW_019211459.1_2748114             0/0
## 2377  NW_019211459.1_2795335            <NA>
## 2378  NW_019211459.1_2795797             0/0
## 2379  NW_019211459.1_2797543             0/0
## 2380  NW_019211459.1_2801120            <NA>
## 2381  NW_019211459.1_2804298            <NA>
## 2382  NW_019211459.1_2832639             0/0
## 2383  NW_019211459.1_2837319             0/0
## 2384  NW_019211459.1_2851154             0/0
## 2385  NW_019211459.1_2852395            <NA>
## 2386  NW_019211459.1_2859310            <NA>
## 2387  NW_019211459.1_2871438             0/0
## 2388  NW_019211459.1_2878083             0/0
## 2389  NW_019211459.1_2878105             0/0
## 2390  NW_019211459.1_2884416             0/0
## 2391  NW_019211459.1_2888854             0/0
## 2392  NW_019211459.1_2924150             0/0
## 2393  NW_019211459.1_2935854            <NA>
## 2394  NW_019211459.1_2947695             0/0
## 2395  NW_019211459.1_2978438             0/0
## 2396  NW_019211459.1_2979143            <NA>
## 2397  NW_019211459.1_2979163            <NA>
## 2398  NW_019211459.1_2981924             0/0
## 2399  NW_019211459.1_2998454             0/0
## 2400  NW_019211459.1_3006130            <NA>
## 2401  NW_019211459.1_3019585            <NA>
## 2402  NW_019211459.1_3032266             0/0
## 2403  NW_019211459.1_3048648            <NA>
## 2404  NW_019211459.1_3060683             0/0
## 2405  NW_019211459.1_3078000             0/0
## 2406  NW_019211459.1_3091163            <NA>
## 2407  NW_019211459.1_3105630             0/0
## 2408  NW_019211459.1_3140777             0/0
## 2409  NW_019211459.1_3162304             0/0
## 2410  NW_019211459.1_3175769            <NA>
## 2411  NW_019211459.1_3197656             0/0
## 2412  NW_019211459.1_3207550            <NA>
## 2413  NW_019211459.1_3218094            <NA>
## 2414  NW_019211459.1_3231117             0/0
## 2415  NW_019211459.1_3232596            <NA>
## 2416  NW_019211459.1_3247637             0/0
## 2417  NW_019211459.1_3251357             0/0
## 2418  NW_019211459.1_3251413             0/0
## 2419  NW_019211459.1_3263437            <NA>
## 2420  NW_019211459.1_3270237             0/0
## 2421  NW_019211459.1_3274126            <NA>
## 2422  NW_019211459.1_3277218             0/0
## 2423  NW_019211459.1_3283014             0/0
## 2424  NW_019211459.1_3285341             0/0
## 2425  NW_019211459.1_3289139            <NA>
## 2426  NW_019211459.1_3289174             0/0
## 2427  NW_019211459.1_3291738            <NA>
## 2428  NW_019211459.1_4021907             0/0
## 2429  NW_019211459.1_4228622             0/0
## 2430  NW_019211459.1_4259114             0/0
## 2431  NW_019211459.1_4262291            <NA>
## 2432  NW_019211459.1_4297197             0/1
## 2433  NW_019211459.1_4310364            <NA>
## 2434  NW_019211459.1_4320315             0/0
## 2435  NW_019211459.1_4397084            <NA>
## 2436  NW_019211459.1_4402585            <NA>
## 2437  NW_019211459.1_4422127            <NA>
## 2438  NW_019211459.1_4434550             0/0
## 2439  NW_019211459.1_4452017            <NA>
## 2440  NW_019211459.1_4457498            <NA>
## 2441  NW_019211459.1_4460795             0/0
## 2442  NW_019211459.1_4464236            <NA>
## 2443  NW_019211459.1_4473615            <NA>
## 2444  NW_019211459.1_4511023             0/0
## 2445  NW_019211459.1_4533525             0/0
## 2446  NW_019211459.1_4543765             0/0
## 2447  NW_019211459.1_4544044             0/0
## 2448  NW_019211459.1_4552610            <NA>
## 2449  NW_019211459.1_4563263             0/0
## 2450  NW_019211459.1_4568845             0/0
## 2451  NW_019211459.1_4599001            <NA>
## 2452  NW_019211459.1_4604999             0/0
## 2453  NW_019211459.1_4623715             0/0
## 2454  NW_019211459.1_4628790            <NA>
## 2455  NW_019211459.1_4629483             0/0
## 2456  NW_019211459.1_4629910             0/0
## 2457  NW_019211459.1_4640567             0/0
## 2458  NW_019211459.1_4648926            <NA>
## 2459  NW_019211459.1_4666027            <NA>
## 2460  NW_019211459.1_4673162             0/0
## 2461  NW_019211459.1_4679058             0/1
## 2462  NW_019211459.1_4730980            <NA>
## 2463  NW_019211459.1_4739012            <NA>
## 2464  NW_019211459.1_4758633             0/0
## 2465  NW_019211459.1_4875501             0/0
## 2466  NW_019211459.1_4955607             0/0
## 2467  NW_019211459.1_4962353            <NA>
## 2468  NW_019211459.1_5000908             0/0
## 2469  NW_019211459.1_5022473            <NA>
## 2470  NW_019211459.1_5033023            <NA>
## 2471  NW_019211459.1_5033099            <NA>
## 2472  NW_019211459.1_5036602             0/0
## 2473  NW_019211459.1_5083040             0/0
## 2474  NW_019211459.1_5128239             0/1
## 2475  NW_019211459.1_5128514             0/0
## 2476  NW_019211459.1_5167203            <NA>
## 2477  NW_019211459.1_5214712            <NA>
## 2478  NW_019211459.1_5275437            <NA>
## 2479  NW_019211459.1_6992579             0/0
## 2480  NW_019211459.1_7023185             0/0
## 2481  NW_019211459.1_8300019            <NA>
## 2482  NW_019211459.1_8899022            <NA>
## 2483  NW_019211459.1_8933702            <NA>
## 2484  NW_019211459.1_8935178             0/0
## 2485  NW_019211459.1_8935884             0/0
## 2486  NW_019211459.1_8937580             0/0
## 2487  NW_019211459.1_8946785            <NA>
## 2488  NW_019211459.1_8958209             0/0
## 2489  NW_019211459.1_8973609             0/0
## 2490  NW_019211459.1_8982284            <NA>
## 2491  NW_019211459.1_9008014             0/0
## 2492  NW_019211459.1_9059347             0/0
## 2493  NW_019211459.1_9241092             0/0
## 2494 NW_019211459.1_10087026            <NA>
## 2495 NW_019211459.1_10107747             0/1
## 2496 NW_019211459.1_10107794             0/1
## 2497 NW_019211459.1_10108541             0/0
## 2498 NW_019211459.1_10112298            <NA>
## 2499 NW_019211459.1_10121231            <NA>
## 2500 NW_019211459.1_10135711            <NA>
## 2501 NW_019211459.1_10146712            <NA>
## 2502 NW_019211459.1_10153761            <NA>
## 2503 NW_019211459.1_10154127            <NA>
## 2504 NW_019211459.1_10160385             0/0
## 2505 NW_019211459.1_10160600             0/0
## 2506 NW_019211459.1_10160707             0/0
## 2507 NW_019211459.1_10169116            <NA>
## 2508 NW_019211459.1_10182778             0/0
## 2509 NW_019211459.1_10198062            <NA>
## 2510 NW_019211459.1_10208896             0/0
## 2511 NW_019211459.1_10241638             0/0
## 2512 NW_019211459.1_10261575             0/0
## 2513 NW_019211459.1_10262865            <NA>
## 2514 NW_019211459.1_10265595            <NA>
## 2515 NW_019211459.1_10267105            <NA>
## 2516 NW_019211459.1_10281952            <NA>
## 2517 NW_019211459.1_10330557             0/0
## 2518 NW_019211459.1_10344060            <NA>
## 2519 NW_019211459.1_10350581             0/0
## 2520 NW_019211459.1_10364949             0/0
## 2521 NW_019211459.1_10366758            <NA>
## 2522 NW_019211459.1_10379197            <NA>
## 2523 NW_019211459.1_10380616             0/0
## 2524 NW_019211459.1_10380950            <NA>
## 2525 NW_019211459.1_10382268            <NA>
## 2526 NW_019211459.1_10382577             0/0
## 2527 NW_019211459.1_10391111             0/0
## 2528 NW_019211459.1_10399252             0/0
## 2529 NW_019211459.1_10433176            <NA>
## 2530 NW_019211459.1_10500630            <NA>
## 2531 NW_019211459.1_10550951            <NA>
## 2532 NW_019211459.1_10589130            <NA>
## 2533 NW_019211459.1_10592546             0/0
## 2534 NW_019211459.1_10606334            <NA>
## 2535 NW_019211459.1_10634012             0/0
## 2536 NW_019211459.1_10640340             0/0
## 2537 NW_019211459.1_10642636            <NA>
## 2538 NW_019211459.1_10646579            <NA>
## 2539 NW_019211459.1_10650724             0/0
## 2540 NW_019211459.1_10655757            <NA>
## 2541 NW_019211459.1_10672227             0/0
## 2542 NW_019211459.1_10674343            <NA>
## 2543 NW_019211459.1_10694592             0/0
## 2544 NW_019211459.1_10723436             0/0
## 2545 NW_019211459.1_10726444            <NA>
## 2546 NW_019211459.1_10727591            <NA>
## 2547 NW_019211459.1_10728287            <NA>
## 2548 NW_019211459.1_10736383            <NA>
## 2549 NW_019211459.1_10749078             0/0
## 2550 NW_019211459.1_10753960             0/0
## 2551 NW_019211459.1_10757714             0/0
## 2552 NW_019211459.1_10797461             0/0
## 2553 NW_019211459.1_10809098             0/0
## 2554 NW_019211459.1_10816645             0/0
## 2555 NW_019211459.1_10828995            <NA>
## 2556 NW_019211459.1_10836412             0/0
## 2557 NW_019211459.1_10839862             0/0
## 2558 NW_019211459.1_10840342            <NA>
## 2559 NW_019211459.1_10844485            <NA>
## 2560 NW_019211459.1_10852954            <NA>
## 2561 NW_019211459.1_10870232            <NA>
## 2562 NW_019211459.1_10873062            <NA>
## 2563 NW_019211459.1_10893246            <NA>
## 2564 NW_019211459.1_10903958            <NA>
## 2565 NW_019211459.1_10913313            <NA>
## 2566 NW_019211459.1_10919810             0/0
## 2567 NW_019211459.1_10920443             0/0
## 2568 NW_019211459.1_10928329             0/0
## 2569 NW_019211459.1_10938695            <NA>
## 2570 NW_019211459.1_10981677             0/0
## 2571 NW_019211459.1_11003838            <NA>
## 2572 NW_019211459.1_11039657            <NA>
## 2573 NW_019211459.1_11048749             0/0
## 2574 NW_019211459.1_11052886            <NA>
## 2575 NW_019211459.1_11061956             0/0
## 2576 NW_019211459.1_11069914            <NA>
## 2577 NW_019211459.1_11071223             0/0
## 2578 NW_019211459.1_11073068            <NA>
## 2579 NW_019211459.1_11086706             0/0
## 2580 NW_019211459.1_11099083            <NA>
## 2581 NW_019211459.1_11126677             0/0
## 2582 NW_019211459.1_11154036            <NA>
## 2583 NW_019211459.1_11157686             0/0
## 2584 NW_019211459.1_11223902             0/0
## 2585 NW_019211459.1_11225067             0/0
## 2586 NW_019211459.1_11316868            <NA>
## 2587 NW_019211459.1_11353270            <NA>
## 2588 NW_019211459.1_11354712            <NA>
## 2589 NW_019211459.1_11357258            <NA>
## 2590 NW_019211459.1_11361710             0/0
## 2591 NW_019211459.1_11374565            <NA>
## 2592 NW_019211459.1_11412125             0/0
## 2593 NW_019211459.1_11414731             0/0
## 2594 NW_019211459.1_11415983             0/0
## 2595 NW_019211459.1_11416038             0/0
## 2596 NW_019211459.1_11420923             0/0
## 2597 NW_019211459.1_11432749             0/0
## 2598 NW_019211459.1_11432839             0/0
## 2599 NW_019211459.1_11452783            <NA>
## 2600 NW_019211459.1_11458297            <NA>
## 2601 NW_019211459.1_11458978            <NA>
## 2602 NW_019211459.1_11472264            <NA>
## 2603 NW_019211459.1_11475340            <NA>
## 2604 NW_019211459.1_11475505             0/0
## 2605 NW_019211459.1_11504195            <NA>
## 2606 NW_019211459.1_12521869            <NA>
## 2607 NW_019211459.1_12563651             0/0
## 2608 NW_019211459.1_12567074             0/0
## 2609 NW_019211459.1_12578233            <NA>
## 2610 NW_019211459.1_14820804            <NA>
## 2611 NW_019211459.1_14831426             0/0
## 2612 NW_019211459.1_14872279             0/0
## 2613 NW_019211459.1_14947379            <NA>
## 2614 NW_019211459.1_14950937            <NA>
## 2615 NW_019211459.1_14951266            <NA>
## 2616 NW_019211459.1_14952178            <NA>
## 2617 NW_019211459.1_14956802            <NA>
## 2618 NW_019211459.1_14962368            <NA>
## 2619 NW_019211459.1_14980326            <NA>
## 2620 NW_019211459.1_15000266             0/0
## 2621 NW_019211459.1_15000613            <NA>
## 2622 NW_019211459.1_15003096             0/0
## 2623 NW_019211459.1_15003526            <NA>
## 2624 NW_019211459.1_15003533            <NA>
## 2625 NW_019211459.1_15008193            <NA>
## 2626 NW_019211459.1_15011154            <NA>
## 2627 NW_019211459.1_15024701             0/0
## 2628 NW_019211459.1_15036510             0/0
## 2629 NW_019211459.1_15077907             0/0
## 2630 NW_019211459.1_15124384            <NA>
## 2631 NW_019211459.1_15155897             0/0
## 2632 NW_019211459.1_15158179            <NA>
## 2633 NW_019211459.1_15159187            <NA>
## 2634 NW_019211459.1_15177454            <NA>
## 2635 NW_019211459.1_15193479            <NA>
## 2636 NW_019211459.1_15198257            <NA>
## 2637 NW_019211459.1_15198999             0/0
## 2638 NW_019211459.1_15200811            <NA>
## 2639 NW_019211459.1_15216901             0/0
## 2640 NW_019211459.1_15216918             0/0
## 2641 NW_019211459.1_15217012             0/0
## 2642 NW_019211459.1_15220950             0/0
## 2643 NW_019211459.1_15224199             0/0
## 2644 NW_019211459.1_15240515             0/0
## 2645 NW_019211459.1_15240547             0/0
## 2646 NW_019211459.1_15249165             0/0
## 2647 NW_019211459.1_15258611            <NA>
## 2648 NW_019211459.1_15259791             0/0
## 2649 NW_019211459.1_15264885             0/0
## 2650 NW_019211459.1_15266513             0/0
## 2651 NW_019211459.1_15298609            <NA>
## 2652 NW_019211459.1_15308432             0/0
## 2653 NW_019211459.1_15322608             0/0
## 2654 NW_019211459.1_15326372             0/0
## 2655 NW_019211459.1_15349553             0/0
## 2656 NW_019211459.1_15377453             0/0
## 2657 NW_019211459.1_15377886            <NA>
## 2658 NW_019211459.1_15378537            <NA>
## 2659 NW_019211459.1_15395119             0/0
## 2660 NW_019211459.1_15417953             0/0
## 2661 NW_019211459.1_15419511            <NA>
## 2662 NW_019211459.1_15449351            <NA>
## 2663 NW_019211459.1_15451862             0/0
## 2664 NW_019211459.1_15478785             0/0
## 2665 NW_019211459.1_16507456             0/0
## 2666 NW_019211459.1_16510496             0/0
## 2667 NW_019211459.1_16510686             0/0
## 2668 NW_019211459.1_16546534             0/0
## 2669 NW_019211459.1_16574418            <NA>
## 2670 NW_019211459.1_16579165            <NA>
## 2671 NW_019211459.1_16579796            <NA>
## 2672 NW_019211459.1_16597717            <NA>
## 2673 NW_019211459.1_16649658             0/0
## 2674 NW_019211459.1_16666296            <NA>
## 2675 NW_019211459.1_16668755            <NA>
## 2676 NW_019211459.1_16705278             0/0
## 2677 NW_019211459.1_16757856             0/0
## 2678 NW_019211459.1_16821859            <NA>
## 2679 NW_019211459.1_16849069             0/0
## 2680 NW_019211459.1_16906948            <NA>
## 2681 NW_019211459.1_17455485             0/0
## 2682 NW_019211459.1_17473661            <NA>
## 2683 NW_019211459.1_17501696             0/0
## 2684 NW_019211459.1_17502670             0/0
## 2685 NW_019211459.1_17504758             0/0
## 2686 NW_019211459.1_17506325            <NA>
## 2687 NW_019211459.1_17550842             0/0
## 2688 NW_019211459.1_17550866             0/0
## 2689 NW_019211459.1_17551060             0/0
## 2690 NW_019211459.1_17704709             0/0
## 2691 NW_019211459.1_17747882             0/0
## 2692 NW_019211459.1_17754516            <NA>
## 2693 NW_019211459.1_17757452            <NA>
## 2694 NW_019211459.1_17757936             0/0
## 2695 NW_019211459.1_17762475             0/0
## 2696 NW_019211459.1_17783200            <NA>
## 2697 NW_019211459.1_17783312            <NA>
## 2698 NW_019211459.1_17853495             0/0
## 2699 NW_019211459.1_17853976            <NA>
## 2700 NW_019211459.1_17868550             0/0
## 2701 NW_019211459.1_17873236            <NA>
## 2702 NW_019211459.1_18278233             0/0
## 2703 NW_019211459.1_18322494             0/0
## 2704 NW_019211459.1_18396482            <NA>
## 2705 NW_019211459.1_18513643            <NA>
## 2706 NW_019211459.1_18518165            <NA>
## 2707 NW_019211459.1_19392838             0/0
## 2708 NW_019211459.1_19657853            <NA>
## 2709 NW_019211459.1_19685137            <NA>
## 2710 NW_019211459.1_19711917            <NA>
## 2711 NW_019211459.1_19736608            <NA>
## 2712 NW_019211459.1_19758847            <NA>
## 2713 NW_019211459.1_19778873             0/0
## 2714 NW_019211459.1_19838151            <NA>
## 2715 NW_019211459.1_19866491            <NA>
## 2716 NW_019211459.1_19866551            <NA>
## 2717 NW_019211459.1_19866562            <NA>
## 2718 NW_019211459.1_21101757            <NA>
## 2719 NW_019211459.1_21104520             0/0
## 2720 NW_019211459.1_21105358            <NA>
## 2721 NW_019211459.1_21108346            <NA>
## 2722 NW_019211459.1_21110045            <NA>
## 2723 NW_019211459.1_21127498            <NA>
## 2724 NW_019211459.1_21138929            <NA>
## 2725 NW_019211459.1_21145296            <NA>
## 2726 NW_019211459.1_21145698            <NA>
## 2727 NW_019211459.1_21147876            <NA>
## 2728 NW_019211459.1_21147886            <NA>
## 2729 NW_019211459.1_21147979            <NA>
## 2730 NW_019211459.1_21172405            <NA>
## 2731 NW_019211459.1_21177339            <NA>
## 2732 NW_019211459.1_21223719            <NA>
## 2733 NW_019211459.1_21225338             0/0
## 2734 NW_019211459.1_21229882            <NA>
## 2735 NW_019211459.1_21230817             0/0
## 2736 NW_019211459.1_21230879             0/0
## 2737 NW_019211459.1_21233465             0/0
## 2738 NW_019211459.1_21277299            <NA>
## 2739 NW_019211459.1_21280179            <NA>
## 2740 NW_019211459.1_21291378            <NA>
## 2741 NW_019211459.1_21291829             0/0
## 2742 NW_019211459.1_21363440             0/0
## 2743 NW_019211459.1_21409066             0/0
## 2744 NW_019211459.1_21560113            <NA>
## 2745 NW_019211459.1_21567705            <NA>
## 2746 NW_019211459.1_22189151            <NA>
## 2747 NW_019211459.1_23806425            <NA>
## 2748 NW_019211459.1_23806429            <NA>
## 2749 NW_019211459.1_24412142             0/0
## 2750 NW_019211459.1_24439494            <NA>
## 2751 NW_019211459.1_24450964             0/0
## 2752 NW_019211459.1_24454464            <NA>
## 2753 NW_019211459.1_24542874            <NA>
## 2754 NW_019211459.1_24590387             0/0
## 2755 NW_019211459.1_24604843             0/0
## 2756 NW_019211459.1_24817465            <NA>
## 2757 NW_019211459.1_24875829            <NA>
## 2758 NW_019211459.1_24877313            <NA>
## 2759 NW_019211459.1_24879194            <NA>
## 2760 NW_019211459.1_24887150            <NA>
## 2761 NW_019211459.1_24902074             0/0
## 2762 NW_019211459.1_24910069            <NA>
## 2763 NW_019211459.1_24919601            <NA>
## 2764 NW_019211459.1_24920560            <NA>
## 2765 NW_019211459.1_24935585             0/0
## 2766 NW_019211459.1_24937340             0/0
## 2767 NW_019211459.1_24943479             0/0
## 2768 NW_019211459.1_24947785             0/0
## 2769 NW_019211459.1_24980552             0/0
## 2770 NW_019211459.1_25026504            <NA>
## 2771 NW_019211459.1_25145986             0/0
## 2772 NW_019211459.1_25146371             0/0
## 2773 NW_019211459.1_25302867            <NA>
## 2774 NW_019211459.1_25307360             0/0
## 2775 NW_019211459.1_25308155             0/0
## 2776 NW_019211459.1_25314277             0/0
## 2777 NW_019211459.1_25327841            <NA>
## 2778 NW_019211459.1_25339976            <NA>
## 2779 NW_019211459.1_25349143             0/0
## 2780 NW_019211459.1_25368088             0/0
## 2781 NW_019211459.1_25402092            <NA>
## 2782 NW_019211459.1_25404651             0/0
## 2783 NW_019211459.1_25421345            <NA>
## 2784 NW_019211459.1_25423717            <NA>
## 2785 NW_019211459.1_25590948             0/0
## 2786 NW_019211459.1_25695764             0/0
## 2787 NW_019211459.1_25697756            <NA>
## 2788 NW_019211459.1_26994249             0/0
## 2789 NW_019211459.1_27024200            <NA>
## 2790 NW_019211459.1_27069604            <NA>
## 2791 NW_019211459.1_27078644            <NA>
## 2792 NW_019211459.1_27088479            <NA>
## 2793 NW_019211459.1_27095147             0/0
## 2794 NW_019211459.1_27101407            <NA>
## 2795 NW_019211459.1_27130898             0/0
## 2796 NW_019211459.1_27141459             0/0
## 2797 NW_019211459.1_27171958            <NA>
## 2798 NW_019211459.1_27176129             0/0
## 2799 NW_019211459.1_27191312            <NA>
## 2800 NW_019211459.1_27195799            <NA>
## 2801 NW_019211459.1_27209794             0/0
## 2802 NW_019211459.1_27313697            <NA>
## 2803 NW_019211459.1_27337349            <NA>
## 2804 NW_019211459.1_27366708            <NA>
## 2805 NW_019211459.1_27397853             0/0
## 2806 NW_019211459.1_28727429             0/0
## 2807 NW_019211459.1_28729763             0/0
## 2808 NW_019211459.1_28729772             0/0
## 2809 NW_019211459.1_28729914            <NA>
## 2810 NW_019211459.1_28731221            <NA>
## 2811 NW_019211459.1_28731819            <NA>
## 2812 NW_019211459.1_28764056             0/0
## 2813 NW_019211459.1_28787309             0/0
## 2814 NW_019211459.1_28823791            <NA>
## 2815 NW_019211459.1_28824457             0/0
## 2816 NW_019211459.1_28840633            <NA>
## 2817 NW_019211459.1_28848882            <NA>
## 2818 NW_019211459.1_28868041             0/0
## 2819 NW_019211459.1_28890187             0/0
## 2820 NW_019211459.1_28890903             0/0
## 2821 NW_019211459.1_28918063             0/0
## 2822 NW_019211459.1_28924314            <NA>
## 2823 NW_019211459.1_28924554            <NA>
## 2824 NW_019211459.1_28935370             0/0
## 2825 NW_019211459.1_28941202             0/0
## 2826 NW_019211459.1_28946902             0/0
## 2827 NW_019211459.1_28947280            <NA>
## 2828 NW_019211459.1_28954103             0/0
## 2829 NW_019211459.1_28960351             0/0
## 2830 NW_019211459.1_28985557             0/0
## 2831 NW_019211459.1_28986755            <NA>
## 2832 NW_019211459.1_28988095            <NA>
## 2833 NW_019211459.1_28989525             0/0
## 2834 NW_019211459.1_29018094            <NA>
## 2835 NW_019211459.1_29030427             0/0
## 2836 NW_019211459.1_29429319            <NA>
## 2837 NW_019211459.1_31582812            <NA>
## 2838 NW_019211460.1_12931634             0/0
## 2839 NW_019211460.1_15044225            <NA>
## 2840 NW_019211460.1_15062783            <NA>
## 2841 NW_019211460.1_15074460            <NA>
## 2842 NW_019211460.1_15080077            <NA>
## 2843 NW_019211460.1_15093868             0/0
## 2844 NW_019211460.1_15099849             0/0
## 2845 NW_019211460.1_15105019             0/0
## 2846 NW_019211460.1_15106952             0/0
## 2847 NW_019211460.1_15111841             0/0
## 2848 NW_019211460.1_15112998            <NA>
## 2849 NW_019211460.1_15137622             0/0
## 2850 NW_019211460.1_15137703            <NA>
## 2851 NW_019211460.1_15142577            <NA>
## 2852 NW_019211460.1_15143847            <NA>
## 2853 NW_019211460.1_15156152             0/0
## 2854 NW_019211460.1_15160358            <NA>
## 2855 NW_019211460.1_15166292            <NA>
## 2856 NW_019211460.1_15166629             0/0
## 2857 NW_019211460.1_15170564             0/0
## 2858 NW_019211460.1_15178081            <NA>
## 2859 NW_019211460.1_15664598            <NA>
## 2860 NW_019211460.1_16839024             0/0
## 2861 NW_019211460.1_16863056            <NA>
## 2862 NW_019211460.1_16866004             0/0
## 2863 NW_019211460.1_16866343            <NA>
## 2864 NW_019211460.1_16866368            <NA>
## 2865 NW_019211460.1_16877968             0/0
## 2866 NW_019211460.1_16885811             0/0
## 2867 NW_019211460.1_16899481            <NA>
## 2868 NW_019211460.1_16904981            <NA>
## 2869 NW_019211460.1_16908001             0/0
## 2870 NW_019211460.1_16913047             0/0
## 2871 NW_019211460.1_16913688             0/0
## 2872 NW_019211460.1_16920377             0/0
## 2873 NW_019211460.1_16922744             0/0
## 2874 NW_019211460.1_16938607             0/0
## 2875 NW_019211460.1_16941948            <NA>
## 2876 NW_019211460.1_16941957            <NA>
## 2877 NW_019211460.1_16941983             0/0
## 2878 NW_019211460.1_16951890             0/0
## 2879 NW_019211460.1_16960747             0/0
## 2880 NW_019211460.1_16962543            <NA>
## 2881 NW_019211460.1_16969978             0/0
## 2882 NW_019211460.1_16969985             0/0
## 2883 NW_019211460.1_16989117             0/0
## 2884 NW_019211460.1_16989571            <NA>
## 2885 NW_019211460.1_18071833            <NA>
## 2886 NW_019211460.1_18105555            <NA>
## 2887 NW_019211460.1_18177915            <NA>
## 2888 NW_019211460.1_18188131            <NA>
## 2889 NW_019211460.1_18206694            <NA>
## 2890 NW_019211460.1_18207348            <NA>
## 2891 NW_019211460.1_19105952             0/0
## 2892 NW_019211460.1_19106353            <NA>
## 2893 NW_019211460.1_19108135             0/0
## 2894 NW_019211460.1_19123354            <NA>
## 2895 NW_019211460.1_19124448            <NA>
## 2896 NW_019211460.1_19131285             0/0
## 2897 NW_019211460.1_19137831            <NA>
## 2898 NW_019211460.1_19152321            <NA>
## 2899 NW_019211460.1_19153419             0/0
## 2900 NW_019211460.1_19165257             0/0
## 2901 NW_019211460.1_22806463             0/0
## 2902 NW_019211460.1_23247293             0/0
## 2903 NW_019211460.1_23344582             0/0
## 2904 NW_019211460.1_23802699            <NA>
## 2905 NW_019211460.1_24002852             0/0
## 2906 NW_019211460.1_24018316             0/0
## 2907 NW_019211460.1_24027395            <NA>
## 2908 NW_019211460.1_24042122            <NA>
## 2909 NW_019211460.1_24053598             0/0
## 2910 NW_019211460.1_24724896            <NA>
## 2911 NW_019211460.1_24825526            <NA>
## 2912 NW_019211460.1_24841003            <NA>
## 2913 NW_019211460.1_24841986             0/0
## 2914 NW_019211460.1_24878441             0/0
## 2915 NW_019211460.1_24885118            <NA>
## 2916 NW_019211460.1_24887517            <NA>
## 2917 NW_019211460.1_24915912             0/0
## 2918 NW_019211460.1_25003419             0/0
## 2919 NW_019211460.1_25014341            <NA>
## 2920 NW_019211460.1_25018135            <NA>
## 2921 NW_019211460.1_25022889             0/0
## 2922 NW_019211460.1_25030682             0/0
## 2923 NW_019211460.1_25045908             0/0
## 2924 NW_019211460.1_25127481            <NA>
## 2925 NW_019211460.1_25128531            <NA>
## 2926 NW_019211460.1_25161208            <NA>
## 2927 NW_019211460.1_25178534             0/0
## 2928 NW_019211460.1_25182753             0/0
## 2929 NW_019211460.1_25195380             0/0
## 2930 NW_019211460.1_25195395             0/0
## 2931 NW_019211460.1_25198752            <NA>
## 2932 NW_019211460.1_25198828            <NA>
## 2933 NW_019211460.1_25209591            <NA>
## 2934 NW_019211460.1_25232531             0/0
## 2935 NW_019211460.1_25247339             0/0
## 2936 NW_019211460.1_25249614            <NA>
## 2937 NW_019211460.1_25256932             0/0
## 2938 NW_019211460.1_25262705            <NA>
## 2939 NW_019211460.1_25276446            <NA>
## 2940 NW_019211460.1_25276553            <NA>
## 2941 NW_019211460.1_25276669             0/0
## 2942 NW_019211460.1_25278561             0/0
## 2943 NW_019211460.1_25287719             0/0
## 2944 NW_019211460.1_25288753            <NA>
## 2945 NW_019211460.1_25295444            <NA>
## 2946 NW_019211460.1_25307870             0/0
## 2947 NW_019211460.1_25320490             0/0
## 2948 NW_019211460.1_25323653             0/0
## 2949 NW_019211460.1_25339648             0/0
## 2950 NW_019211460.1_25350152            <NA>
## 2951 NW_019211460.1_25361016            <NA>
## 2952 NW_019211460.1_25368991             0/0
## 2953 NW_019211460.1_25379747            <NA>
## 2954 NW_019211460.1_25384542             0/0
## 2955 NW_019211460.1_25396698            <NA>
## 2956 NW_019211460.1_25416104             0/0
## 2957 NW_019211460.1_25424159             0/0
## 2958 NW_019211460.1_25433021             0/0
## 2959 NW_019211460.1_25439608            <NA>
## 2960 NW_019211460.1_25439872             0/0
## 2961 NW_019211460.1_25573640            <NA>
## 2962 NW_019211460.1_25574326            <NA>
## 2963 NW_019211460.1_25576696            <NA>
## 2964 NW_019211460.1_25577312             0/0
## 2965 NW_019211460.1_25577476            <NA>
## 2966 NW_019211460.1_25579474            <NA>
## 2967 NW_019211460.1_26171391             0/0
## 2968 NW_019211460.1_26171400             0/0
## 2969 NW_019211460.1_26182516             0/0
## 2970 NW_019211460.1_26187482            <NA>
## 2971 NW_019211460.1_26196068             0/0
## 2972 NW_019211460.1_26197331             0/0
## 2973 NW_019211460.1_26198191             0/0
## 2974 NW_019211460.1_26232280             0/0
## 2975 NW_019211460.1_26232919             0/0
## 2976 NW_019211460.1_26235727             0/0
## 2977 NW_019211460.1_26243136             0/0
## 2978 NW_019211460.1_29665252            <NA>
## 2979 NW_019211460.1_29716761            <NA>
## 2980 NW_019211460.1_31235175            <NA>
## 2981 NW_019211460.1_31244709            <NA>
## 2982 NW_019211460.1_31269312             0/0
## 2983 NW_019211460.1_31272915            <NA>
## 2984 NW_019211460.1_31320858             0/0
## 2985 NW_019211460.1_31376306             0/0
## 2986 NW_019211460.1_31380075            <NA>
## 2987 NW_019211460.1_31387222             0/0
## 2988 NW_019211460.1_31387393            <NA>
## 2989 NW_019211460.1_31389335             0/0
## 2990 NW_019211460.1_31393856             0/0
## 2991 NW_019211460.1_31398725             0/0
## 2992 NW_019211460.1_31399415             0/0
## 2993 NW_019211460.1_31410941             0/0
## 2994 NW_019211460.1_31414455            <NA>
## 2995 NW_019211460.1_31419681             0/0
## 2996 NW_019211460.1_31419694             0/0
## 2997 NW_019211460.1_31425021            <NA>
## 2998 NW_019211460.1_31427558             0/0
## 2999 NW_019211460.1_31429025             0/0
## 3000 NW_019211460.1_31431478             0/0
## 3001 NW_019211460.1_31432913             0/0
## 3002 NW_019211460.1_31436979            <NA>
## 3003 NW_019211460.1_31438471             0/0
## 3004 NW_019211460.1_31442110            <NA>
## 3005 NW_019211460.1_31488643            <NA>
## 3006 NW_019211460.1_31489789             0/0
## 3007 NW_019211460.1_31491343            <NA>
## 3008 NW_019211460.1_31492247             0/1
## 3009 NW_019211460.1_31493843            <NA>
## 3010 NW_019211460.1_31495429             0/0
## 3011 NW_019211460.1_31503669             0/0
## 3012 NW_019211460.1_31503834             0/0
## 3013 NW_019211460.1_31513600            <NA>
## 3014 NW_019211460.1_31597890            <NA>
## 3015 NW_019211460.1_31599116             0/0
## 3016 NW_019211460.1_31601492             0/0
## 3017 NW_019211460.1_31610769            <NA>
## 3018 NW_019211460.1_31648212            <NA>
## 3019 NW_019211460.1_31648968             0/0
## 3020 NW_019211460.1_31652326            <NA>
## 3021 NW_019211460.1_31674644            <NA>
## 3022 NW_019211460.1_31675965             0/0
## 3023 NW_019211460.1_31676774             0/0
## 3024 NW_019211460.1_31696355             0/0
## 3025 NW_019211460.1_31703516            <NA>
## 3026 NW_019211460.1_31713732             0/0
## 3027 NW_019211460.1_31722091            <NA>
## 3028 NW_019211460.1_31722242             0/0
## 3029 NW_019211460.1_31728935            <NA>
## 3030 NW_019211460.1_31741624             0/0
## 3031 NW_019211460.1_32089125            <NA>
## 3032 NW_019211460.1_32760140             0/0
## 3033 NW_019211460.1_32783801             0/0
## 3034 NW_019211460.1_32785599            <NA>
## 3035 NW_019211460.1_32791825            <NA>
## 3036 NW_019211460.1_32791840            <NA>
## 3037 NW_019211460.1_32997494             0/0
## 3038 NW_019211460.1_33010381             0/0
## 3039 NW_019211460.1_33032400             0/0
## 3040 NW_019211460.1_33053490            <NA>
## 3041 NW_019211460.1_33057515            <NA>
## 3042 NW_019211460.1_33060141             0/0
## 3043 NW_019211460.1_33062371             0/0
## 3044 NW_019211460.1_33062504             0/0
## 3045 NW_019211460.1_33062519             0/0
## 3046 NW_019211460.1_33066045            <NA>
## 3047 NW_019211460.1_33076181            <NA>
## 3048 NW_019211460.1_33097115             0/0
## 3049 NW_019211460.1_33287487             0/0
## 3050 NW_019211460.1_33287839            <NA>
## 3051 NW_019211460.1_33304817             0/0
## 3052 NW_019211460.1_33309273             0/0
## 3053 NW_019211460.1_33310506             0/0
## 3054 NW_019211460.1_33312790             0/0
## 3055 NW_019211460.1_33313369            <NA>
## 3056 NW_019211460.1_33359694            <NA>
## 3057 NW_019211460.1_33366977             0/0
## 3058 NW_019211460.1_33379971             0/0
## 3059 NW_019211460.1_33404378             0/0
## 3060 NW_019211460.1_33415732            <NA>
## 3061 NW_019211460.1_33428353            <NA>
## 3062 NW_019211460.1_33429143            <NA>
## 3063 NW_019211460.1_33433042            <NA>
## 3064 NW_019211460.1_33495253             0/0
## 3065 NW_019211460.1_33522404             0/0
## 3066 NW_019211460.1_33522600             0/0
## 3067 NW_019211460.1_33527995            <NA>
## 3068 NW_019211460.1_33529980            <NA>
## 3069 NW_019211460.1_33530202             0/0
## 3070 NW_019211460.1_33534021            <NA>
## 3071 NW_019211460.1_33535617            <NA>
## 3072 NW_019211460.1_33917439            <NA>
## 3073 NW_019211460.1_34014687             0/0
## 3074 NW_019211460.1_34041959            <NA>
## 3075 NW_019211460.1_34043457             0/0
## 3076 NW_019211460.1_34058147            <NA>
## 3077 NW_019211460.1_34058870             0/0
## 3078 NW_019211460.1_34062816            <NA>
## 3079 NW_019211460.1_34070680             0/0
## 3080 NW_019211460.1_34126479            <NA>
## 3081 NW_019211460.1_34137993             0/0
## 3082 NW_019211460.1_34138303             0/0
## 3083 NW_019211460.1_34142318             0/0
## 3084 NW_019211460.1_34154512            <NA>
## 3085 NW_019211460.1_34165787             0/0
## 3086 NW_019211460.1_34205199             0/1
## 3087 NW_019211460.1_34217942             0/0
## 3088 NW_019211460.1_34218229             0/0
## 3089 NW_019211460.1_34232725            <NA>
## 3090 NW_019211460.1_34245101            <NA>
## 3091 NW_019211460.1_34261994             0/0
## 3092 NW_019211460.1_34262025             0/0
## 3093 NW_019211460.1_34268652            <NA>
## 3094 NW_019211460.1_34291766            <NA>
## 3095 NW_019211460.1_34292705            <NA>
## 3096 NW_019211460.1_34294018             0/0
## 3097 NW_019211460.1_34300158             0/0
## 3098 NW_019211460.1_34300503             0/0
## 3099 NW_019211460.1_34301950            <NA>
## 3100 NW_019211460.1_34308067             0/0
## 3101 NW_019211460.1_34309708             0/0
## 3102 NW_019211460.1_34318077             0/0
## 3103 NW_019211460.1_34324796            <NA>
## 3104 NW_019211460.1_34384750             0/0
## 3105 NW_019211460.1_34586356            <NA>
## 3106 NW_019211460.1_36084650             0/0
## 3107 NW_019211460.1_36089631             0/0
## 3108 NW_019211460.1_36106639            <NA>
## 3109 NW_019211460.1_36107351            <NA>
## 3110 NW_019211460.1_36116166            <NA>
## 3111 NW_019211460.1_36143041             0/0
## 3112 NW_019211460.1_36151630             0/0
## 3113 NW_019211460.1_36156137            <NA>
## 3114 NW_019211460.1_36170384             0/0
## 3115 NW_019211460.1_36183796             0/0
## 3116 NW_019211460.1_36185700             0/0
## 3117 NW_019211460.1_36186688             0/0
## 3118 NW_019211460.1_36189677             0/0
## 3119 NW_019211460.1_36415971            <NA>
## 3120 NW_019211460.1_36427509             0/0
## 3121 NW_019211460.1_36440396             0/0
## 3122 NW_019211460.1_36495785             0/0
## 3123 NW_019211460.1_36495846             0/0
## 3124 NW_019211460.1_36496401             0/0
## 3125 NW_019211460.1_36505988            <NA>
## 3126 NW_019211460.1_36506890             0/0
## 3127 NW_019211460.1_36506925             0/0
## 3128 NW_019211460.1_36506986             0/0
## 3129 NW_019211460.1_36530456             0/0
## 3130 NW_019211460.1_36537102             0/0
## 3131 NW_019211460.1_36540823             0/0
## 3132 NW_019211460.1_36574100             0/0
## 3133 NW_019211460.1_36575932            <NA>
## 3134 NW_019211460.1_36583026            <NA>
## 3135 NW_019211460.1_36611771             0/0
## 3136 NW_019211460.1_36620185             0/0
## 3137 NW_019211460.1_36659362            <NA>
## 3138 NW_019211460.1_36682242             0/0
## 3139 NW_019211460.1_36701049             0/0
## 3140 NW_019211460.1_36719644            <NA>
## 3141 NW_019211460.1_36734940            <NA>
## 3142 NW_019211460.1_36735078             0/0
## 3143 NW_019211460.1_36741319            <NA>
## 3144 NW_019211460.1_36742951             0/0
## 3145 NW_019211460.1_36772888             0/0
## 3146 NW_019211460.1_36805687            <NA>
## 3147 NW_019211460.1_36826968             0/0
## 3148 NW_019211460.1_36827037             0/0
## 3149 NW_019211460.1_36831567            <NA>
## 3150 NW_019211460.1_36840509             0/0
## 3151 NW_019211460.1_36840579            <NA>
## 3152 NW_019211460.1_36845653             0/0
## 3153 NW_019211460.1_36850173            <NA>
## 3154 NW_019211460.1_36857838             0/0
## 3155 NW_019211460.1_36865521            <NA>
## 3156 NW_019211460.1_36871446            <NA>
## 3157 NW_019211460.1_36890857            <NA>
## 3158 NW_019211460.1_36932253             0/0
## 3159 NW_019211460.1_36960392            <NA>
## 3160 NW_019211460.1_36961695            <NA>
## 3161 NW_019211460.1_36963691            <NA>
## 3162 NW_019211460.1_36983052             0/0
## 3163 NW_019211460.1_37022696             0/0
## 3164 NW_019211460.1_37025392             0/0
## 3165 NW_019211460.1_37027355            <NA>
## 3166 NW_019211460.1_37030423            <NA>
## 3167 NW_019211460.1_37033340             0/0
## 3168 NW_019211460.1_37033695             0/0
## 3169 NW_019211460.1_38389758             0/0
## 3170 NW_019211460.1_38423641             0/0
## 3171 NW_019211460.1_38425397            <NA>
## 3172 NW_019211460.1_38501125            <NA>
## 3173 NW_019211460.1_38502540            <NA>
## 3174 NW_019211460.1_38534624             0/0
## 3175 NW_019211460.1_38581841            <NA>
## 3176 NW_019211460.1_38670263             0/0
## 3177 NW_019211460.1_38670283             0/0
## 3178 NW_019211460.1_38670376             0/0
## 3179 NW_019211460.1_38706950            <NA>
## 3180 NW_019211460.1_38721194            <NA>
## 3181 NW_019211460.1_38721858             0/0
## 3182 NW_019211460.1_38730978            <NA>
## 3183 NW_019211460.1_38757613             0/0
## 3184 NW_019211460.1_38777645            <NA>
## 3185 NW_019211460.1_38787839            <NA>
## 3186 NW_019211460.1_38858664             0/0
## 3187 NW_019211460.1_38860985            <NA>
## 3188 NW_019211460.1_38862382             0/1
## 3189 NW_019211460.1_39010532            <NA>
## 3190 NW_019211460.1_39017021            <NA>
## 3191 NW_019211460.1_39017101            <NA>
## 3192 NW_019211460.1_39026115            <NA>
## 3193 NW_019211460.1_39044960             0/0
## 3194 NW_019211460.1_39045225            <NA>
## 3195 NW_019211460.1_39076063             0/0
## 3196 NW_019211460.1_39150362            <NA>
## 3197 NW_019211460.1_39155077             0/0
## 3198 NW_019211460.1_39185495             0/0
## 3199 NW_019211460.1_39299544            <NA>
## 3200 NW_019211460.1_39301652             0/0
## 3201 NW_019211460.1_39380485            <NA>
## 3202 NW_019211460.1_39395023            <NA>
# join all grandkids of 26 families together
grnson0.0_0.0 <- reduce(list_grnson, full_join, by = "site") %>% view()

grndat0.0_0.0 <- reduce(list_grndat, full_join, by = "site") %>% view()

# it looks like the 'false sites' in one family, are not necessarily the false sites in other families. should i take them out?